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Abstract 

This project develops comprehensive modeling and simulation tools for analysis of variable rotor 
speed helicopter propulsion system dynamics. The Comprehensive Variable-Speed Rotorcraft Propulsion 
Modeling (CVSRPM) tool developed in this research is used to investigate coupled rotor/engine/fuel 
control/gearbox/shaft/clutch/fhght control system dynamic interactions for several variable rotor speed 
mission scenarios. In this investigation, a prototypical two-speed Dual-Clutch Transmission (DCT) is 
proposed and designed to achieve 50 percent rotor speed variation. The comprehensive modeling tool 
developed in this study is utilized to analyze the two-speed shift response of both a conventional single 
rotor helicopter and a tiltrotor drive system. In the tiltrotor system, both a Parallel Shift Control (PSC) 
strategy and a Sequential Shift Control (SSC) strategy for constant and variable forward speed mission 
profiles are analyzed. Under the PSC strategy, selecting clutch shift-rate results in a design tradeoff 
between transient engine surge margins and clutch frictional power dissipation. In the case of SSC, clutch 
power dissipation is drastically reduced in exchange for the necessity to disengage one engine at a time 
which requires a multi-DCT drive system topology. In addition to comprehensive simulations, several 
sections are dedicated to detailed analysis of driveline subsystem components under variable speed 
operation. In particular an aeroelastic simulation of a stiff in-plane rotor using nonlinear quasi-steady 
blade element theory was conducted to investigate variable speed rotor dynamics. It was found that 2/rev 
and 4/rev flap and lag vibrations were significant during resonance crossings with 4/rev lagwise loads 
being directly transferred into drive-system torque disturbances. To capture the clutch engagement 
dynamics, a nonlinear stick-slip clutch torque model is developed. Also, a transient gas-turbine engine 
model based on first principles mean-line compressor and turbine approximations is developed. Finally an 
analysis of high frequency gear dynamics including the effect of tooth mesh stiffness variation under 
variable speed operation is conducted including experimental validation. Through exploring the 
interactions between the various subsystems, this investigation provides important insights into the 
continuing development of variable-speed rotorcraft propulsion systems. 

1.0 Introduction 

1.1 Introduction and Motivation 

According to a recent NASA-Army-Industry-University investigation, significant benefits to 
rotorcraft operational performance, effectiveness and acoustic signature could be gained through the 
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(a) 


Figure 1.1. — Some current rotorcraft employing a variable speed rotor 
(achieved with engine speed variations); (a) Bell-BoeingV-22, (b) 

Eurocopter EC130, (c) A160 Humming Bird UAV. 

ability to adjust main rotor speed to accommodate various flight conditions (Ref. 1). In particular, variable 
speed rotor technology is critical to the slowed-rotor compound configuration concepts and would offer 
significant benefits to future Heavy-Lift helicopter and tiltrotor configurations as well as variable 
diameter rotor concepts (Ref. 2). Current rotorcraft propulsion systems are designed around a fixed-ratio 
transmission without the capability to vary rotor speed except by engine speed adjustment. Figure 1.1 
shows several examples of rotorcraft which are required to adjust their rotor speeds over a certain portion 
of their operating envelop. In these cases, the rotor RPM variations are archived through engine speed 
adjustments alone. 

Since the specific fuel consumption of modern gas-turbine engines is optimum within a relatively 
narrow speed range, the ability to achieve a variable speed rotor through engine speed adjustment is 
limited (Ref. 3). Therefore, to fully benefit from variable rotor speed designs, some form of variable ratio 
transmission becomes necessary. 

1.2 V ariable Ratio T r ansmissions 

Recently, there have been many investigations concerning variable ratio transmissions for rotorcraft 
(Refs. 2 to 4) and other applications, such as wind-turbines and automotive drivetrains (Refs. 5 to 10). See 
to Reference 3 for an extensive survey of variable speed transmission concepts. For rotorcraft 
applications, both two-speed transmission and traction-drive Continuously Variable Transmission (CVT) 
concepts based on Split Torque Differential Planetary arraignments have been explored in References 3 to 
5. Furthermore, others have explored, nontraction based, Pericyclic CVT’s (P-CVT) with very promising 
results (Refs. 2, 6, and 7). In particular, Reference 2 explored a P-CVT design to replace a conventional 
fixed-ratio planetary gear transmission. Here, a P-CVT achieved a maximum single stage reduction of 
50:1 and could be continuously varied to approximately 25:1. Furthermore, the P-CVT was favorable in 
terms of power density, part count, and load-sharing (high contact ratios) compared with a conventional 
fixed-ratio planetary gear transmission design. The success of these investigations adds further motivation 
to develop and implement variable speed rotor rotorcraft concepts. 

There have been numerous studies concerning various, ideally driven, propulsion system components. 
To note a few, References 1 1 to 16 explored gear mesh induced vibration in gear trains, References 17 to 
22 explored rotor blade dynamics and aeromechanics, References 23 to 30 studied flexible driveshaft 
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vibration and stability issues, and References 3 1 to 34 explored closed-loop fuel control and modeling 
techniques for gas turbine engines. Also, many coupled rotorcraft engine/drivetrain/rotor system analyses 
have been conducted for constant rotor speed propulsion systems, see References 35 to 38. 

Since many rotordynamic characteristics in the propulsion system, such as rotor aerodynamic 
damping and centrifugal stiffness (Ref. 17), cross-shafting driveline misalignment parametric variations 
(Refs. 25 to 27), and gear mesh stiffness effects (Refs. 1 1 and 12), are RPM dependent, variable speed 
operation and speed shifts will give rise to nonlinear and parametric interaction mechanisms not captured 
by the previous propulsion system investigations. Furthermore, since all variable speed transmission 
arrangements are fundamentally split-path power devices, the effects of power circulation within the 
system must be included. In some cases, the circulating power can be higher than the transmission input 
power which leads to excessive loads and gear train vibrations (Refs. 3 and 4). Thus, both the variable 
RPM effects and the power circulation induced vibration phenomena must be investigated and clearly 
understood in the context of the overall engine/transmission/driveline/rotor system. 

1.3 Research Objectives and Significance 

The overall goal of this program is to develop first principles-based, comprehensive rotorcraft 
propulsion system modeling and analysis tools to account for dynamical interactions between engine, 
rotor, cross-shaft, clutch and gearbox subsystems. 

The variable speed rotorcraft propulsion system components and their interactions which are 
considered in this project are summarized in Figure 1.2. Through the comprehensive analysis, another 
goal is to gain insight into the complex system-level and component-level transient responses under a 
variety of variable speed operating conditions and gear shifting scenarios. 

To fully benefit from the variable speed rotor concept and ensure efficient, reliable, highly loaded 
propulsion system designs, a complete understating of the effects of: (a) variable speed operation and (b) 
variable ratio transmissions designs on the overall propulsion system dynamics must be obtained. To 
address these issues, the objective of the proposed research is to analytically develop and experimentally 
validate a comprehensive dynamics model of a variable speed rotorcraft propulsion system including 
engine/transmission/cross-shafting/rotor interactions. Incorporating recent developments such as, low- 
order turbine engine throttle response models, multisegment flexible driveshaft models, lumped parameter 
multistage gear train models, coupled flap-lag-torsion rotor models, and variable ratio transmission 
concepts into a single integrated analysis will advance the state of the art and yield new knowledge about 
the overall variable speed propulsion system dynamics. 

By accounting for the various propulsion subsystem dynamic interactions, the proposed Comprehensive 
Variable Speed Rotorcraft Propulsion System Model (CVSRPM) will allow better performance prediction, 
give new design insights and enable system level optimization of (a) conventional constant speed helicopter 
propulsion systems, (b) novel variable speed rotorcraft propulsion systems, and (c) various multipath power 
flow configurations (e.g., quad tiltrotors, tandem rotor/pusher configurations, etc.). Furthermore, this 
research will create analytical and computational tools for analyzing the resulting nonlinear, time -varying, 
comprehensive propulsion system model. Specifically, tools for predicting stability limits, vibration 
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Figure 1.2. — Comprehensive variable speed rotorcraft propulsion system elements. 
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amplitudes, shaft and gear tooth stresses, heat generation and parasitic losses will be developed. Finally, 
given the current trend toward tandem, tilt-rotor, multirotor, and co-axial rotor/pusher-prop configurations 
with numerous cross-shafting and multipath power flow arrangements (e.g., Boeing CH-47 and V-22, 
Sikorsky X-2 High Speed Lifter and Heavy Lift quad tiltrotor concepts), the CVSRPM code will be 
utilized to investigate and compare these different configurations for fixed and variable speed design 
cases under a variety of operating conditions. 
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2.0 Variable RPM Rotor 

2.1 Summary 

An aeroelastic simulation of a stiff in-plane rotor in forward fight was conducted to investigate the 
dynamic characteristics of a variable speed rotor during resonance crossing. A finite element analysis 
based on a moderate deflection beam model was employed to capture the coupled flap, lag and torsion 
deflections of rotor blade. The nonlinear quasi-steady blade element theory with table look-up of airfoil 
aerodynamics was utilized to calculate the blade aerodynamic loads. By using Hamilton's principle, 
system equations of motion were derived based on the generalized force formulation. An implicit 
Newmark integration scheme was used to calculate the steady and transient responses. Transient 
aeroelastic responses of a four-blade stiff in-plane rotor are calculated to analyze the blade lagwise root 
bending moment and rotor torque. Rotor systems with identical and dissimilar blades were investigated. 
During the 2/rev resonance crossing, for identical rotors, the transient lagwise root bending moment is 
amplified significantly. The variation of rotor torque is substantially small. Flap motion has vital 
contribution to the steady and transient lagwise loads. The faster the blade crosses the resonance area, the 
smaller the transient lagwise loads and the higher the rotor torque. For dissimilar rotors, 5 percent 
reduction of one blade mass at 0.6 to 0.7R can cause a sharp rise of 2/rev rotor torque. Increasing blade 
lag critical damping from 1 to 5 percent can reduce the peak-peak lagwise root bending moment by 
64.9 percent, and the rotor torque is reduced to the level without dissimilarity. 4/rev lagwise root bending 
moment is transferred to the rotor shaft during the 4/rev resonance crossing. 

Transient aeroelastic response of a stiff in-plane rotor system undergoing variable speed operation in 
forward flight is simulated using a finite element model. During crossing of the fundamental lag mode 
near 2/rev, high transient lag bending moments are observed. Flapping amplitude and duration of the 
resonance crossing event have a strong influence on the peak lagwise root bending moments. Embedded 
chordwise fluidlastic dampers are introduced to control the transient lagwise loads of the variable speed 
rotor during resonance crossing. The design of the fluidlastic damper is based on the analysis of a two 
degree-of-freedom blade-damper system. Results indicate more than 6 percent critical damping can be 
provided to the blade around the resonance rotor speed. Results indicate that approximately 65 percent 
peak-peak moment reduction can be achieved with reasonable damper devices (i.e., less than 5 percent 
blade mass). The stroke of the damper is limited to less than 2.5 percent blade chord length in the worst 
case scenarios (i.e., high flapping). Parametric studies show that tuning port area ratios, loss factors, and 
device mass can be utilized to enhance the performance of the damper, and control the stroke. Damper 
performance is shown to be relatively independent of rotor thrust and advance ratio. 

Rigid blade modeling with coupled flap and lag motions is derived for the comprehensive 
rotor/transmission/engine modeling. The degree of the variation of rotor speed is considered. A nonlinear 
quasi-steady aerodynamic model is utilized, and the lift, drag, and moment coefficients of the aerofoil are 
calculated by a two-dimensional table-look-up method. The Pitt-Peters dynamic inflow model is utilized 
to capture the induced velocity over the rotor during steady and transient states. 

2.2 Lagwise Dynamic Analysis of a Variable Speed Rotor 
2.2.1 Introduction 

The reduction of rotor speed in forward fight provides a feasible and effective means to achieve better 
efficiency and performance compared to constant rotor speed helicopters (Ref. 39), especially for long 
endurance or long range acquirement. The reduction of rotor speed can also reduce rotor noise, and 
improve the life of transmission systems, gears and engines. XV- 15 attempted to adopt a two-speed rotor 
(589 to 517 RPM, 100 to 88 percent) (Ref. 40), which caused rotor dynamic problems at other rotor speed — 
high vibration and loads. The V-22 Osprey also had two-speed rotor: 412 RPM (100 percent) used for 
helicopter mode and for conversion to airplane mode; 333 RPM (81 percent) used for propeller mode in 
forward flight. The rotor speed of the A 160 unmanned helicopter can be slowed down to about 40 percent 


NAS A/CR— 20 13-21 6502 


5 



r 


e 



Figure 2.1. — Rigid blade single degree-of-freedom lag model. 

of its maximum value. Light and stiff blades are utilized in A 160 rotor system to avoid vibration 
problems. To attain the technology goal of NASA Vehicle Systems Program, Johnson et al. examined in 
depth three rotorcraft configurations for the large civil transport concept (Ref. 1). These configurations all 
employed variable speed rotors. The rotor speed can vary in the range of 70.0 to 37.7 RPM (100 to 53.8 
percent), 80.9 to 25.5 RPM (100 to 31.5 percent) and 68.6 to 26.9 RPM (100 to 39.2 percent) for LCTR, 
LCTC, and LABC respectively. The study of variable speed rotors indicates that increasing the range of 
the variation of the rotor speed is beneficial and advantageous. 

With the increase of the variation of rotor speed, some problems associated with the rotor system will 
emerge, such as vibration, loads, stability and so on. For example, the fundamental natural frequency of the 
lag motion for uniform rotor blade with lag hinge offset and hinge spring shown in Figure 2. 1 is (Ref. 41) 

9 3 e k£ , . 

Vq _ 21^7 + /,,Q 2 (l-£>) 

Typical values of e in normal conditions are small. If rotor speed decreases by 50 percent, the 
fundamental lag frequency will increase by about 100 percent. For typical soft in-plane rotors ( v^= 0.65- 
0.80/rev in full rotor speed), the rotor blades go through the 1/rev resonance area. For typical stiff in-plane 
rotors (Vf= 1.4-1.6/rev in full rotor speed), they go through the 2/rev resonance area. 

Normally, the nCl(n = 0, 1,2...) loads from blades will be transferred to the transmission system and 
fuselage in steady state. When a rotor goes through a resonance area, the vibration and loads usually 
increase sharply or dramatically, and a large rise in amplitude will occur. During this transient process, if 
these sharply amplified loads are transferred to the transmission system, they will seriously affect the 
working condition of transmission system. This transfer might cause the damage to transmission shafts, 
gears or engines, especially for dissimilar rotors. 

Present research focuses on the analysis of lagwise root bending moment and rotor torque during 
resonance crossing. System modeling of helicopter rotor, taking into account the variation of rotor speed, 
will be presented first. Hence, a hypothetical stiff in-plane rotor is proposed to investigate the dynamic 
characteristics during resonance crossing. Transient aeroelastic responses of blades and rotor torque will 
be calculated to analyze which harmonic component of the amplified transient lagwise root bending 
moment is transferred to transmission system, especially for dissimilar rotors. 

2.2.2 Analytical Model 

The objective of present report is to investigate the dynamic characteristics of variable speed rotor 
during resonance crossing. To describe the transient process, the major difference between the modeling 
of variable speed rotors and the modeling of constant speed rotors is the consideration of the variation of 
rotor speed. That causes the modification of kinetic terms of rotor blade associated with rotor’s rotation 
degree. Since resonance crossing is a transient process, dynamic inflow model needs to be utilized to 
capture the variation of rotor induced velocity. It is assumed that the blade is undergoing moderately large 
deflection and small strains in flap, lag and torsion. The detail modeling procedure can be found in detail 
in (Ref. 42). 
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2.2.2. 1 


Moderate Deflection Beam Model 


To describe the geometrical nonlinearity of advanced helicopter blades, such as hingeless and 
bearingless rotor blades, Hodges and Dowell put forward the moderate deflection beam model (Ref 43). 
The modified moderate deflection beam model is adopted in this paper (Ref. 44). The axial strain and 
shear strains of arbitrary point (x, 77 , g) on the blade are 


a 

V 

1 

, 1^- 
1 

+ w" (/) 

-z 

tt 1 . tt 

-v (p + W 

fi-^y 


l 2 ; 




V 2 )_ 


= «' + -y + -^- + (^0 + | fo 2 + ? 2 V ' 2 - y 


Yx V =(^-?)(^ + vV ) 

Yxg = Wc + v) (/ + v"w ) 

and the variation of the elastic potential energy is 

SU = X " = | = [WJyEs^SSxx + GYxti^Yxtj + Gy xg 8y xg )dA dL 


( 2 ) 

(3) 

(4) 

(5) 


2.2.2.2 Kinetic Energy 

Usually helicopter rotors have complicated kinematics. Even for hingeless rotor, the blades undergo 
elastic deformations, rigid motions introduced by pitch controls and rotation around rotor shaft. To describe 
the nonlinear coupling characteristics between elastic deflections and rigid rotations for rotor blades, the 
computational method for kinetic energy derived by Zheng (Ref. 45) based on the generalized force 
formulation is employed. The rigid rotations of the flap, lag and pitch hinges are introduced as generalized 
coordinates, as shown in Figure 2.2. The hinge sequences can be changed according to actual rotors. When 
the three rigid rotational generalized coordinates are adopted, the modeling methodology can suit other 
types of helicopter rotors. When calculating the kinetic energy, warping effect is usually not taken into 
account. The position vector of an arbitrary point on the blade in a rotor shaft coordinate frame is 
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( 6 ) 




k]h/fe/][7>fc] 


where these transformation matrices are defined in the Appendix. Thus the variation of the kinetic energy 
of a blade is 


ST = X!l, Q" & i> = ZLiI \\ A ~ p *'it dA dlSq 




and its ith generalized force introduced by the kinetic energy is 


(7) 


( 8 ) 
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According to the definition of tangent mass, damping, and stiffness matrices introduced by kinetic energy, 
the expressions of theses matrices are 


mtJZL 

J dq. 


0R 0R „ 

p dA dl 

iJJa Sqj dq, 
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dq i dq, 


dA dl 


Kl = 


dQj 

dq 


- = -f ff J— — + R- 
JiJJa [ dq. Qq. 


5 2 R 




dq,dq 


dA dl 


J J 


( 9 ) 

( 10 ) 

( 11 ) 


From the expressions of tangent mass, damping, and stiffness matrices and the generalized force vector, it 
can be seen that, if theses matrices are calculated, the position vector of an arbitrary point on the blade 
shown in Equation (6) and its derivative with respect to time and partial derivatives with respect to 
generalized coordinates need to be given out. 

2.2.2.3 Aerodynamics 

A nonlinear quasi-steady aerodynamic model is adopted, and the lift, drag, and moment coefficients 
of the airfoil are calculated by a two-dimensional table-look-up method according to the angle of attack 
and the oncoming air flow (Mach number). The direction of blade section velocity is shown in Figure 2.3. 
The velocity of an arbitrary point on the pitch axis with respect to the local airflow is 


U R 

T 

Ax 
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U T 

> = < 

Ay 
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Up 


Az 




V,Al T l T A T 'A T frl T rs\} T 


( 12 ) 


where p u and p 3 are the components of the air velocity in the rotor tip plane. Using a nonlinear quasi- 
steady aerodynamic model, the forces introduced by the aerodynamics can be calculated using the 
velocity expression, Equation (12). The variation of the work done by the aerodynamics is 


dW A = X'l, QfSqi = | ; (F^ •<SR i +M^ -Sa s )dl 


(13) 
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Figure 2.3. — Velocity components and cross-section definitions. 


The generalized force introduced by the aerodynamics is 



da/ 

d<7, , 


dl 


(14) 


It should be noted that the force vector F A , moment vector M A , position vector R s , and angle vector 
are defined in the rotor shaft coordinate system, which is treated as an inertial coordinate system. Thus, 
after the aerodynamic forces and moments in the deformed coordinate system are calculated, these 
vectors must be transformed to the rotor shaft coordinate system. 

2.2.2 A Inflow Modeling 

The three-state dynamic inflow model is used to determine the inflow distribution over the rotor disk 
during the steady and transient states (Ref. 46). The inflow with the variation of axial position and 
azimuth is 


r r 

A = Aq + A c — cos^ + A S — sin^ 
R R 


where the coefficients are determined by 
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(15) 


(16) 


the coefficient terms shown in the upper equation can be found in (Ref. 46). The wind tunnel trim method 
is applied to calculate the pitch controls during steady states. 
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2.2.2.5 Equations of Motion 

By using Hamilton’s principle, the implicit nonlinear dynamic equations based on the generalized 
force formulation include three parts: elastic potential energy, kinetic energy and work done by the 
aerodynamics. The equations of motion are 

ef(q)+ef(q 5 q,q^) + e^(q,q^)=° = ( 17 ) 


an Implicit-Newmark integration method is utilized to calculate the steady and transient responses 
(Ref 47). This unconditionally stable implicit scheme permits the use of large time integration step. This 
step is determined by the considerations of accuracy. 


2.2.2.6 Wind-Tunnel Trim Analysis 

When the steady-state response is calculated, the corresponding pitch controls to the rotor need to be 
provided. In the present report, the wind-tunnel trim method (Ref. 48) is adopted to trim the rotor. The 
control input vector \sx={6 0 0 lc OisY \ and the output or target vector is y={C T fii c PhY - According to the 
Newton-Raphson method, the recursive expression to calculate the pitch controls in the nth step is 


x„ =x„_ 1 +J- 1 (y-y„_ 1 ) 

where the Jacobian matrix is 
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(18) 


(19) 


because the precision of the Jacobian matrix only influences the number of the iterations, simple 
expressions of the partial derivatives are adopted to give approximations. For example, the physical 
expression for the thrust coefficient (Ref. 49) is 
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according to Equation (20), the partial derivatives associated with the thrust coefficient are 
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( 21 ) 


the other partial derivatives can also be derived in the same way. 

2.2.2 .7 Steady Loads Calculation 

Trimmed values of pitch controls need to be found before the calculation of the steady rotor response. 
At first, some prescribed pitch controls are initialized. The wind-tunnel trim method is used to update the 
pitch controls after every circle integration. Then, the steady response can be attained after several circle 
iterations. Because the externally applied forces, including the centrifugal force and inertial force, are 
reacted by the structure, the root bending moments are calculated using the generalized structural forces 
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corresponding to the degrees of freedom. For example, the variation of the elastic potential energy can be 
expressed as 


SU^F&i i = (22) 

where the generalized degrees are independent. If the structural lagwise root bending moment in the blade 
coordinate frame is desired, the moment is the generalized nodal force F v , corresponding to the degree v ’ 
at the blade root. Because of pitch controls, transformation needs to be conducted to transform the 
bending moments from the blade coordinate frame to the hub coordinate frame. At last, the 1 and 2/rev 
lagwise root bending moments from the periodic response of the lagwise root bending moment are 
extracted. 

2.2.3 Results and Discussions 

2.2.3. 1 Parameters of Baseline Rotor 

For numerical studies, a four-blade stiff in-plane hingeless rotor shown in Figure 2.4 is adopted as a 
hypothetical example. The parameters of the rotor are given in Table 2.1. Every rotor blade is discretized 
by six fifteen degree-of-freedom beam elements. The rotor blade frequencies are shown in Figure 2.5. ‘F’ 
denotes the frequency in flapwise direction; ‘L’ denotes the frequency in lagwise direction and ‘T’ 
denotes the frequency in torsional direction. The total operational rotor speed range is from 150 to 
300 RPM. The first lag frequency goes through the 2/rev resonance area, which is the lowest resonance 
frequency. The variation of the nondimensional fundamental flap and lag frequencies with respect to the 
rotor speed is shown in Figure 2.6, which illustrates that the 2/rev resonance occurs at the rotor speed 
205 RPM (3.42 Hz, 21.5 rad/s). The variation of the lag frequency ratio with respect to the rotor speed 
varies more significantly than that of the flap frequency ratio. Due to the low resonance crossing 
frequency and low damping in lagwise direction, present research concentrates on the lagwise resonance 
crossing analysis and the loads transferred to the transmission system. 



Figure 2.4. — Finite element modeling of four-blade elastic hingeless rotor 
(wind tunnel trim condition). 


NAS A/CR— 20 13-21 6502 


11 


Non-dimensional Fundamental Lag Frequency 


TABLE 2. L— BASELINE ROTOR PROPERTIES 


Parameter Value 

Number of blades 4 

Rotor radius 6.0 m 

Blade mass 60.0 kg 

Blade chord 0.6 m 

Blade coning angle 0° 

Airfoil NACA0012 

Blade linear twist 0° 

Full rotor speed 300 RPM 

Flapwise bending stiffness 1.68><10 5 Nm 2 

Lagwise bending stiffness 1.85><10 6 Nm 2 

T orsional stiffness 8.157xl0 4 Nm 2 



Figure 2.5. — Rotor blade frequencies. 




Figure 2.6. — Fundamental frequency versus rotor speed. 
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2.2.3.2 Steady Loads 

The steady lagwise root bending moments with respect to the rotor speed are shown in Figure 2.7. 

One percent critical structure damping is added to the fundamental lag mode. The wind tunnel trim with 
zero flapping is applied to this analysis. The thrust is constant with respect to the rotor speed, and the 
corresponding thrust coefficient at 150 RPM is 0.010. The forward speed is also constant with respect to 
the rotor speed, and the advance ratio at 150 RPM is 0.30. The 1/rev moment increases significantly when 
the rotor speed approaches 150 or 300 RPM. There is a peak value of the 2/rev lagwise root bending 
moment around 200 RPM, which is close to the 2/rev resonance frequency (205 RPM). The loads near the 
resonance frequency are much higher than the loads at other rotor speeds. The 3/rev loads near 150 RPM 
is much higher, since the first lag frequency ratio at 150 RPM is close to 3.0. 

Due to the Coriolis force, blade flap motion has vital contribution to lagwise loads. With 1° 
longitudinal or lateral flapping, the steady lagwise root bending moments with respect to the rotor speed 
are shown in Figure 2.8 and Figure 2.9. The overall amplitudes increase significantly with the previous 
case. The peak moment near the 2/rev resonance frequency for 1° longitudinal or lateral flap is about 3.03 
or 1.38 times the value with zero flapping. With different fundamental flap frequency for zero flapping, 
the flap wise and lagwise root bending moments with respect to the rotor speed are shown in Figure 2.10 
and Figure 2.1 1. The increase of the blade flap frequency can decrease the lagwise loads significantly, 
and the flapwise loads increase distinctly as the penalty. 



Harmonic 


Figure 2.7. — Steady lagwise bending moment with zero flapping 
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Figure 2.8. — Steady lagwise root bending moment with 1° longitudinal flap. 



Figure 2.9. — Steady lagwise root bending moment with 1° lateral flap. 
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3/rev Load [Nm] 2/rev Load [Nm] 1 /rev Load [Nm] 3/rev Load [Nm] 2/rev Load [Nm] 1 /rev Load [Nm] 



Figure 2.10. — Flapwise loads with different flap stiffness. 



Figure 2.1 1 . — Lagwise loads with different flap frequencies. 
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For rigid blade modeling in hover, the relations between the flap motion and the cyclic pitch controls 
are (Ref. 49) 


e lc = *-{vj- 1)/? 1C + /? 1S (23) 

e ls = B -{vj-i)p ls -p 1 c (24) 

Usually, Vp is very close to 1.0. For different flap frequency ratios, the cyclic pitch controls are shown in 
Figure 2.12. The Lock number 7.58 is adopted in the calculation. The increase of fundamental flap 
frequency ratio causes significant increase of pitch controls. The limitation of pitch controls causes trim 
problems. The insufficient flap motion from the pitch controls makes the helicopter alter its altitude to 
generate enough forces and moments to trim the helicopter. 

From the above analysis, the blades undergo high lagwise loads when operating around the lagwise 
resonance area. For the consideration of blade strength and fatigue, it is of great importance to avoid long 
time working at these rotor speeds. In present research, the rotor speed will go through from 180 to 
240 RPM, or inversely, to prevent rotor blades from high loads. With the consumption of fuel, the 
optimum rotor speed decreases slowly with the helicopter gross weight (Ref. 50). When the rotor speed 
reaches the boundary of the resonance area, it goes through the dangerous region very quickly to avoid 
high loads, as seen in Figure 2.13. The baseline crossing time is 10 sec. 


JL =1°,p. =0° 

r Is 



Rap Frequency Ratio 



Flap Frequency Ratio 


Figure 2.12. — Relations between pitch controls and flap motion in hover. 
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Figure 2.14. — The time history of the rotor speed and its variation. 


2.2.3.3 Identical Rotor Crossing 2/rev Resonance Area 

Baseline results are obtained for a forward flight condition, fi = 0.30 at half full rotor speed, and the 
corresponding thrust coefficient 0.010. When the rotor speed changes, the forward speed and rotor thrust 
are kept constant, but the advance ratio and thrust coefficient will alter. The longitudinal tilt of the rotor 
tip plane can be simply calculated by 


Pic — 


C T pnR 2 (D.R ) 2 



(25) 


flnR 2 typically falls between 0.004 (clean designs) and up to 0.025 (heavy-lift transport or first-generation 
helicopters) (Ref. 49). In present example, 0.004 is adopted, and the corresponding fi\ c is 1.03°. In the 
baseline example, 1° longitudinal flap is applied. The lateral tilt is set to zero. The time histories of the 
rotor speed and the variation of rotor speed (fl) are shown in Figure 2.14. The variation of rotor speed 
increases firstly then decreases. Every 10 RPM, the pitch controls are provided by the steady trim. The 
pitch controls between the discrete values are calculated by the curve-fit method. The pitch controls 
accompanying the variation of rotor speed are shown in Figure 2.15. 


NAS A/CR— 20 13-21 6502 


17 






The time histories of the lagwise root bending moment of blade -2 and the rotor torque are shown in 
Figure 2.16. The lagwise root bending moment is amplified significantly during 10 to 12 sec. The FFT 
(Fast Fourier Transform) analysis illustrates that the resonance frequency (6.84 Hz/410 RPM) is exactly 
twice the rotor speed at that instant, which means the resonance frequency is 2/rev. The corresponding 
time histories of the thrust and the flap wise displacement of the blade tip are shown in Figure 2.17. 
During the resonance crossing, the variation of the rotor thrust and the flapwise tip displacement are 
substantially small. 

The variation of rotor torque is substantially small during the resonance crossing. The shape of the 
rotor torque is the same as the rotor speed. The transient 2/rev lagwise root bending moment has not been 
transferred to the hub. The static rotor torque at the low rotor speed (180 RPM) is much smaller than that 
at high rotor speed (240 RPM), which shows that reducing rotor speed in forward flight can reduce the 
required rotor power. It is economical to decrease rotor speed in forward flight. 
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Figure 2.16. — The lagwise root moment and rotor torque with 1° longitudinal flap. 
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Figure 2.18. — The time histories of the lagwise root bending moment 
and the rotor torque without the variation of rotor speed. 


During the transient process, there is a sharp increase of the rotor torque. Figure 2.18 shows the time 
history of the rotor torque without the variation of rotor speed (fl). The sharp increase of the static rotor 
torque disappears. The inertia of the rotor is 2880 kg/m 2 . At the instant 10.62 sec, the angular acceleration 
is 0.968 rad/sec 2 . The extra torque introduced by the rotor inertia is 2788 Nm. This value is about the 
difference between the transient torque with the variation of rotor speed at that instant and the torque 
without the variation. To reduce the sharp increase of the rotor torque, it is necessary to limit the 
acceleration of rotor rotation (fl). 

2.2.3.3.1 Influence of Flap Motion 

Previous analysis has pointed out that flap motion has significantly influence on steady lagwise loads. 
Figure 2.19 and Figure 2.20 show the time histories of the lagwise root bending moment and rotor torque 
for different longitudinal flapping. The flap motion has significant influence on the transient peak-peak 
moment. The peak-peak lagwise root bending moment with 1°, 2° and 3° longitudinal flapping is 3.13, 
7.77 and 14.01 times that with zero flapping, respectively. With the increase of the flapping, the variation 
of the rotor torque during the resonance crossing also increases significantly. The FFT analysis of the 
rotor torque with 3° longitudinal flapping during the resonance crossing illustrates that the dominative 
frequency of the peak torque is 14.17 Hz. That is four times the frequency of the rotor speed at that instant. 
The 4/rev lagwise load is transferred to the rotor hub, since large flapping causes large higher harmonic 
loads. The level of rotor flapping during the resonance crossing is very important to the peak loading in 
the lead-lag direction. It is important to consider the dynamic trim or proper time histories of cyclic and 
collective pitch controls to control the resonance loads. 
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Figure 2.19. — Lagwise root bending moment for different flapping. 





Figure 2.20. — Rotor torque for different flapping. 


2.2.3.3.2 Influence of the Duration During the Transient Process 

The acceleration of rotor rotation is determined by the duration from the initial instant of the rotor 
speed to the terminal instant. Figure 2.21 shows the responses of the lagwise root bending moment for 
different durations. The lagwise root bending moment increases with the increase of the duration. The 
peak-peak moment increases by 37.2 percent, when the duration increases from 5 to 20 sec. Figure 2.22 
shows the time histories of the rotor torque with different durations. The maximum rotor torque decreases 
significantly with the increase of the duration. The maximum value decreases by 46.3 percent, when the 
period varies from 5 to 20 sec. Long duration causes high lagwise loads, and short duration causes high 
rotor torque due to the large variation of rotor speed. It is necessary to balance the duration to cross the 
resonance area. 
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Figure 2.21 . — Lagwise root bending moment for different rotor 
speed variation durations. 




Figure 2.22. — Rotor torque for different rotor speed variation duration. 


2.2.3.3.3 Influence of Forward Speed 

Helicopter forward speed is limited by the shock wave in advancing blades and stall in retreating 
blades. The reduction of the rotor speed can help to alleviate the shock wave. But it aggravates the stall. 
In present example, the maximum advance ratio at half full rotor speed is 0.4 (136 km/h), and the 
corresponding advance ratio at the full rotor speed is 0.2. The peak-peak lagwise root bending during the 
resonance crossing for different forward speeds are shown in Figure 2.23. The moment increases 
significantly with the forward speed. The influence of the flapping is obvious. For the consideration of 
blade strength, varying the rotor speed at lower forward speed could help to reduce the loads during the 
resonance crossing. 
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Figure 2.23. — The influence of forward speed on the 
peak-peak lagwise root bending moment. 



Figure 2.24. — The influence of rotor thrust on the lagwise 
root bending moment. 


2.2.3.3.4 Influence of Rotor Thrust 

The influence of the rotor thrust on the peak-peak lagwise root bending moment during the resonance 
crossing is shown in Figure 2.24. The forward speed is 102 km/h with wind tunnel trim. The baseline thrust 
is corresponding to the thrust coefficient 0.010 at half rotor speed. Thirty percent increase of the rotor thrust 
causes 50.7 and 21.3 percent increase of the peak-peak moment for 0° and 1° longitudinal flapping 
respectively. Thirty percent decrease of the rotor thrust causes 40.6 and 21.1 percent reduction of the 
moment. It is better to change the rotor speed at low thrust to reduce transient loads during resonance 
crossing. 
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Figure 2.25. — Dissimilarity at 0.6 to 0.7R. 


2.2.3.4 Dissimilar Rotor Crossing 2/rev Resonance Area 

Usually, the manufacturing tolerance among all blades of a rotor is limited to some small value to 
achieve static and dynamic balance. Not only manufacturing tolerance, but also some other factors 
introduce rotor dissimilarities — battle damage to military helicopter rotors, moisture absorption, loss of 
trim mass, a nonoperational lag damper and so on (Refs. 51 to 53). These blade faults usually result in a 
change of mass distribution, bending stiffness, torsional stiffness, lag damping, and so on. In this work, 
the practice in (Ref. 53) is followed to deal with the rotor dissimilarities. The baseline dissimilarity is 
simulated by the modification of blade material property at 0.6 to 0.7R, as seen in Figure 2.25. 

2.2.3.4.1 Unbalanced Mass 

The time history of the lagwise root bending moment and the rotor torque are shown in Figure 2.26, 
when the blade-2 mass distribution is reduced by 5 percent at 0.6 to 0.7R. The wind tunnel trim with 1° 
longitudinal flapping is applied. In the following analysis, this trim is adopted. The variation of the 
lagwise root bending moment is substantially small compared with Figure 2.16. The rotor torque 
increases significantly during the resonance crossing. The difference of mass distribution among different 
blades introduces unbalanced lagwise root bending moments, and the unbalanced moments are transferred 
to the rotor hub, then generate sharp rise in the torque response. The frequency of the major component of 
the transient response is 6.84 Hz. Obviously, this frequency is the 2/rev component at that instant. 

2.2.3.4.2 Influence of Blade Stiffness 

The influence of the flap, lag and torsional stiffness on the transient rotor torque is substantially small, 
when the stiffness is reduced by 5 percent at 0.6 to 0.7R. 

2.2.3.4.3 Influence of Blade Lag Damping 

To investigate the influence of the blade damping on the transient loads, the critical damping of every 
blade is increased from 1 to 5 percent. The time responses of the lagwise root bending moment and the 
rotor torque are shown in Figure 2.27, when the blade-2 mass distribution is reduced by 5 percent at 0.6 to 
0.7R. The increase of the blade lag damping can significantly reduce the transient loads during the 
resonance crossing. The peak-peak lagwise root bending moment is reduced by 64.9 percent. The rotor 
torque is reduced to the level without the dissimilarity. The increase of the blade lag damping is an 
effective means to control the transient loads during the lagwise resonance crossing. However, 
introducing 5 percent critical lag damping into stiff in-plane blades is a difficult task due to the small 
rotation of blade root and the high blade stiffness. 


NAS A/CR— 20 13-21 6502 


24 




Frequency [Hz] 


Figure 2.26. — Lagwise root bending moment and rotor toque with 
5 percent mass dissimilarity. 




Figure 2.27. — Lagwise root bending moment and rotor torque with 
5 percent critical damping. 
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2.2.3A.4 Dissimilarity in Two Blades 

In the previous analysis, one blade property is modified. When the mass distribution of blade-2 and 
blade-3 is reduced by 5 percent simultaneously at 0.6 to 0.7R, the time history of the rotor torque is shown 
in Figure 2.28. The rotor torque is almost the same as the baseline rotor torque shown in Figure 2.16 during 
the 2/rev resonance crossing. When the mass distribution of blade-2 and blade -4 is reduced by 5 percent at 
0.6 to 0.7R, the time history of the rotor torque is also shown in Figure 2.28. It is amplified significantly. 
From the above analysis, it can be concluded that the modifications of two blade mass distributions can 
cause the lagwise root bending moments to cancel each other, or inversely to superimpose together 
depending on their phase relation. 

The principles determining the phenomenon during steady states are applied to explain why the 
lagwise root bending moments superimpose together when the mass distribution of blade-2 and blade -4 is 
modified; on the contrary, they seem to cancel each other when the mass distribution of blade-2 and 
blade-3 is altered. The lagwise root bending moment with baseline blade property can be expressed as 

M — M 0 + I£ = i [M nc cos (nxp m ) + M ns sin(mp m )] (26) 

the lagwise root bending moment with modified blade property can be expressed as 

M' — M' 0 + 22=it M 2c COS + M; s sin(n^ m )] (27) 

since the maximum harmonic component of the lagwise root bending moment is 2/rev during the 2/rev 
resonance crossing, the analysis of 2/rev load is specified. When the mass distribution of blade-2 and 
blade-3 is reduced, the 2/rev rotor torque is 

Q 2 = M 2c cos(2j/> m ) + M 2s sin(2^ m ) + M 2c cos 2 (ip m + f) + M 2s sin 2 (ip m + 

+M 2c cos 2(xp m +n) + M 2s sin 2(xp m + n) +M 2c cos 2 (xp m + ™) + M 2s sin 2(i/> m + = 0 (28) 

if the mass distribution of blade-2 and blade -4 is reduced, the 2/rev rotor torque is 

Q 2 = 2(M 2c - M' c ) cos(2 xp m ) + 2 (M 25 - M' 2s ) sin(2 (29) 

Equation (28) illustrates that the lagwise root bending moments of blade -2 and blade-3 cancel each other 
due to the phase difference 7i/2. The phenomenon illustrated in Equation (29) is on the contrary. 


Blade-2 and Blade-3 



Blade-2 and Blade-4 



Time [sec] 

Figure 2.28. — Rotor torque with two blade dissimilarity. 
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2.2.3.5 Higher Stiffness Rotor During 4/rev Resonance Crossing 

To investigate the dynamics of variable speed rotor during 4/rev resonance crossing, the lag stiffness 
of the blade shown in Table 2.1 is increased to four times its original value. Other parameters of the rotor 
are kept at the same values. The variation of the fundamental lag frequency with respect to the rotor speed 
is shown in Figure 2.29, which illustrates that the rotor goes through the 3/rev (270 RPM) resonance area 
and the 4/rev (202 RPM) resonance area. 

The time histories of the lagwise root bending moment of blade-2 and the rotor torque are shown in 
Figure 2.30. The variation of the lagwise root bending moment is substantially small. There is a sharp rise 
of rotor torque during 8 to 10 sec. The FFT analysis shows that the resonance frequency is 13.19 Hz, 
which is about four times 3.36 Hz. The 4/rev lagwise root bending moment is transferred to the rotor hub. 
The steady response of the rotor torque during the first 5 sec is much larger than that shown in Figure 
2.14. It is obvious that stiffer blades introduce larger variation of rotor torque. 

The time histories of the blade-2 lagwise root bending moment and rotor torque are shown in Figure 
2.3 1, when the mass distribution of blade-2 is reduced by 5 percent at 0.6 to 0.7R using the higher lag 
stiffness. The variation of the lagwise root bending moment is substantially small compared with the 
above response. The rotor torque is similar to that of the identical rotor. 

The critical damping of every blade is increased to 5 times the original value of the identical rotor. 

The time responses of the lagwise root bending moment and the rotor torque are shown in Figure 2.32. 
The increase of blade lag critical damping can significantly reduce the transient loads during the 
resonance crossing. The maximum rotor torque is decreased by 24.6 percent. 
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Figure 2.30. — Baseline lagwise root bending moment when crossing 4/rev 
resonance area. 


2000 



Time [sec] 


10000 

^ B0 00 
z 

g eooo 

o- 

^ 4000 
2000 

0 5 10 15 20 

Time [sec] 

Figure 2.31 . — Response with higher lag stiffness and 5 percent mass 
dissimilarity. 
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Figure 2.32. — Reponses with five times the original damping. 

2.2.4 Summary and Conclusions 

In this research, an aeroelastic simulation is employed to analyze the dynamic characteristics of a 
four-blade stiff in-plane rotor during resonance crossing. Based on the 2/rev resonance crossing analysis, 
the following can be concluded: 

• For identical rotors, 2/rev lagwise root bending moment of single rotor blade is amplified 
significantly during the resonance crossing, and the variation of rotor torque is substantially small. 

• The flap motion contributes significantly to the lagwise loads, whether in steady state or in 
transient process. One degree longitudinal flap causes 213 percent increase of the peak-peak 
lagwise root bending moment compared with the zero flapping case. It is important to reduce the 
flapping during the resonance crossing to control the transient loads. 

• Increasing the duration during the resonance crossing amplifies the peak-peak lagwise root 
bending moment, and decreases the maximum rotor torque. 

• Varying the rotor speed at lower forward speed or lower thrust could help to reduce the transient 
loads. 

• Five percent reduction of blade mass at 0.6 to 0.7R has significant influence on the transient rotor 
torque. Five percent blade critical damping significantly reduces the rotor torque to the level 
without dissimilarity. The peak-peak lagwise root bending moment is reduced by 64.9 percent. 

• The 5 percent reduction of flap lag, or torsional stiffness at 0.6 to 0.7R has substantially small 
influence on the transient lagwise root bending moment and rotor torque. 

• The dissimilarity in two blades can cancel each other, or inversely superimpose together 
depending on their phase relation. 

• Based on the 4/rev resonance crossing analysis, the following can be concluded: 

• For identical rotor, 4/rev lagwise root bending moment is transferred to the rotor shaft, and a 
sharp rise of 4/rev rotor torque occurs. 

• For dissimilar rotor with 5 percent reduction of one blade mass distribution at 0.6 to 0.7R, 4/rev 
lagwise root bending moment is also transferred to the rotor shaft. 

• Increasing blade lag damping significantly reduces the maximum transient rotor torque. 
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From the summary of present study, it can be concluded that attention should be paid to blade 
strength during resonance crossings. It is seriously dangerous for dissimilar rotors to go through 
resonance crossing, which perhaps will cause some damages to transmission shafts, gears and engines. 


2.3 Transient Loads Control of a Variable Speed Rotor during Lagwise Resonance 
Crossing 

2.3.1 Introduction 

For some variable speed rotors, lagwise loads can increase sharply during the lagwise resonance 
crossing (Ref. 54). For example, the fundamental natural frequency of the lag motion for uniform rotor 
blade with lag hinge offset and hinge spring is 


v 2 — 3 e 

^ 2 1—e Ib^ 2 { l-e) 


(30) 


Typical values of e in normal conditions are small. If rotor speed decreases by 50 percent, the 
fundamental lag frequency will increase by about 100 percent. For typical stiff in-plane rotors (v^= 1.4- 
1.6/rev in full rotor speed), they go through the 2/rev resonance area. Rotor blades have to undergo these 
high transient loads. The transfer of these severe loads to transmission systems affects the working 
condition of the transmission systems, which might damage the transmission shafts, gears or engines. It is 
highly necessary to control these high transient resonance loads to protect rotor blades and transmission 
systems. 

For articulated or soft in-plane rotors, large lag dampers are attached to blade root. These lag dampers 
can enhance the stability of the rotor- fuselage system and can also suppress the sudden rise of lagwise 
transient loads. For stiff in-plane rotors, the ground or air resonance problems disappear. Lag dampers are 
not typically required from these rotor systems. For stiff in-plane variable speed rotors during large 
variation of rotor speed, how to manage the high transient loads during resonance crossing can be an 
important and challenging problem to solve. Usually, structural damping is very low, and the lag damping 
introduced by aerodynamics is weak. For stiff in-plane rotors, it is difficult to provide enough lag 
damping through traditional root dampers due to the small deformations at blade root, and large blade lag 
stiffness. 

Zapfe and Lesieutre put forward the concept of distributed inertial dampers using chordwise 
absorbers to provide damping to the beam structure (Ref. 55). Results demonstrated the effectiveness 
using a simply supported beam under tensile load. Hebert and Lesieutre utilized highly distributed tuned 
vibration absorbers to provide lag damping to helicopter blades (Ref. 56). Their studies illustrated that 
method could provide enough damping with a weight penalty equal to only 3 percent blade mass. Kang, 
Smith and Lesieutre carried out a comprehensive analysis of the damping characteristics of a rotor blade 
with an embedded chordwise absorber (Refs. 57 to 59). The potential large stroke of elastomeric dampers 
was a major obstacle for the application of the elastomeric dampers due to the limited space in the blade 
cavity. The concept of fluidlastic dampers was introduced to reduce the stroke. Petrie, Lesieutre and 
Smith utilized embedded fluid elastomer absorber to enhance blade lag damping, and conducted the fluid- 
elastomer design and the lag damping test of a blade with a fluid-elastomer absorber (Ref. 60 and 61). 

Han and Smith introduced embedded chordwise absorbers to reduce blade lagwise loads, and their 
investigations indicated that placing an embedded chordwise absorber at blade tip could significantly 
reduce the 1/rev or 2/rev lagwise load (Ref. 57). The embedded chordwise fluidlastic damper shown in 
Figure 2.33 presents a feasible means to provide enough lag damping to stiff in-plane rotors. 
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Figure 2.33. — Configuration of embedded chordwise damper 


2.3.2 Research Objectives 

This research focuses on the simulation and control of the transient loads of a variable speed rotor 
during resonance crossing. The feasibility of embedded chordwise fluidlastic dampers is explored. A 
device design procedure to select parameters such as mass, loss factor, tuning frequency, and tuning port 
area ratio is developed. The system modeling of helicopter rotor with embedded chordwise dampers is 
presented. A hypothetical stiff in-plane rotor is analyzed during resonance crossing. Transient aeroelastic 
responses are calculated to evaluate the performance of the fluidlastic dampers to reduce the transient 
blade lagwise loads. Comprehensive parametric study including tuning frequency, loss factor, tuning port 
ratio, tuning mass, blade flapping motion, duration during resonance and flight states on the performance 
of embedded fluidlastic dampers are investigated. 


2.3.3 Analytical Model 

2.3.3. 1 Modeling of a Rigid Blade With a Fluidlastic Damper 

The viscoelastic behavior of the elastomeric spring of the fluidlastic damper in frequency domain is 
simulated by 


K = kp(l + irf) 


(31) 


The configuration of the mathematical modeling of a fluidlastic damper is shown in Figure 2.34. The 
inertia of the fluidlastic damper can be amplified by the tuning mass through the leverage effect. Larger G 
can be utilized to reduce the weight and stroke. That is the vital benefit of fluidlastic dampers different 
from elastomeric dampers (Ref. 57). The natural frequency of the fluidlastic damper is 


0) n — 



k v 

+(G-l) 2 m t 


(32) 


To evaluate how much damping is provided to the rotor blade, a two degree-of-freedom model of a 
rigid blade with a fluidlastic damper is employed from Petrie’s work (Ref. 62). The equations of motion 
for the system are 


M 




= F 


(33) 
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Figure 2.34. — Mathematical modeling of fluidlastic damper. 


Where 


M = 


fe m(r - e) 2 dr + (m p + m t )(r - e ) 2 ~[m p - m t (G - l)](r a - e) 


~[m p - m t (G - l)](r a - e) 


m p + m t (G - l) 2 


K = 


r / e R m(r — e)eft 2 dr + + (m p + m t )(r — e)eil 2 ~[ m p ~ m t(G ~ l)]efl 2 

~[m p - m t (G - l)]efi 2 k* a - [m p + m t {G - l) 2 ]fi 2 


F = 


[C F f(. r ~ e ) dr + m p eCl 2 a p0 + m f efi 2 a t0 j 
( m p Sl 2 a p0 — m t (G — l)n 2 a t0 ] 


(34) 

(35) 

(36) 


The eigenvalues and eigenvectors of the coupled system are analyzed to calculate the blade damping. To 
evaluate the stroke of the fluidlastic damper, 1° blade tip motion is given to calculate the stroke of the 
damper. 

Under static state, the spring acts on the prime mass and the fluid has no effect on the static 
displacement of the damper. In this way, the static displacement of the prime mass can be reduced 
significantly due to the large spring. Under rotation, the embedded damper must undergo large centrifugal 
force. The fluid is encapsulated in the damper with physical restriction, which can limit the variation of 
the mass center of the damper. 

To estimate how much the in-plane loads can be reduced during the resonance crossing, the rotor 
response is calculated in the time domain using a comprehensive rotor code. The harmonic excitations in 
the lagwise direction are integral times the rotor speed. During the resonance crossing, the major 
excitation is the resonance frequency. In this report, the baseline variable speed rotor goes through the 
2/rev resonance area, and the 2/rev rotor speed (2Q) is utilized as the excitation frequency (co) to the 
damper. The corresponding steady solution for the single excitation can be expressed by Ae ia *. A is the 
amplitude of the response. For the single degree-of-freedom modeling of the fluidlastic damper, the 
relation for the equivalent of the imaginary part is 


c p o) = r]k v (37) 

From the definition of the critical damping ^ = c p /[2a) n (m p + (G — 1 ) 2 m t )], the relation between the 
critical damping and the excitation frequency is 

= S ( 38 > 

The above equation illustrates that the critical damping decreases with the rotor speed. 
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2.3.3.2 Comprehensive Rotor Modeling With Fluidlastic Dampers 

The comprehensive rotor modeling follows References 54 and 42. A moderate deflection beam model 
is employed to describe the elastic deformations of rotor blades. Rigid rotations including blade hinges 
and rotor rotation degree are introduced based on the generalized force formulation (Ref. 45). Quasi- 
steady aerodynamics with table look-up aerofoil aerodynamics is utilized to describe the blade 
aerodynamics. The rotor induced velocity in the steady and transient states is captured by the Pitt-Peters 
dynamic inflow model. 

For trusted analysis, a time domain model of fluidlastic damper device is required. One degree-of- 
freedom model of the fluidlastic damper is coupled with the rotor modeling. The device is mechanically 
constrained to move only in the blade chordwise direction. The relative displacements of the prime mass 
and the tuning mass of the absorber to the blade are denoted as a p +a p0 and a t +a t0 . The position vector of 
the prime mass can be expressed as 


PX \ 
py t ~ 


PZ j 


l of 




l lp 


(x + u\ 


0 [T rs ] + \ 0 [Tfr][Trs] + \ 0 [T lf ][T fr ][T rs ] + \ v [T pl ][T lf ][T fr ][T rs ] + 


0 


0 


0 


iv 


0 


j dp + UpQ | T ] \T pi 1 1 | | T f r 1 1 T rs ] 

y o 

The relation of the prime mass with the tuning mass is 

a t = -(G - 1 )Op 

Then the relative displacement of the tuning mass to the blade can be expressed as 


(49) 


(50) 



lfl) 


H P \ 


(x + u\ 


oj [r rs ] + |oj MPYJ + j 0 J [Tif][T fr ][T rs ] + ^ v J [T pl ][T lf }[T fr ][T rs \ + 

0 0 ) T 

— (G-l)a p + a t0 [T][T pl ][T lf ][T fr ][T rs ] (51) 

0 ) 


The tangent mass, stiffness, damping matrices and generalized force vector introduced by the kinetic 
energy can be calculated according to the definition in References 42 and 45 using the above position 
vectors. The elastic potential energy is and the work done by the viscous force is \c v c$. The 
fluidlastic damper thus contributes to the generalized force, tangent stiffness and damping matrices. 

Assembling the three components including structure, kinetics and aerodynamics, the equations of 
motion based on the generalized force formulation are achieved. The implicit Newmark integration 
method is utilized to calculate the steady and transient responses in time domain. Every rotor blade is 
discretized by the 15 degree-of-freedom beam elements, as seen in Figure 2.35. 
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Figure 2.35. — Rotor-damper system. 


The variation of the rotor speed is prescribed during the resonance crossing in forward flight. To 
maintain the rotor thrust and forward speed, pitch controls need to be provided to the rotor. A dynamic 
wind-tunnel trim is utilized to provide the time histories of the pitch controls. The pitch controls in steady 
state for several discrete rotor speeds are trimmed for the same thrust and forward speed. During the 
variation of the rotor speed, the pitch controls between two rotor speeds are calculated using the curve-fit 
method. 


2.3.4 Embedded Chordwise Damper Design 

The design of embedded chordwise fluidlastic dampers follows Reference 61. The objective in this 
section is to decrease the transient lagwise loads during the resonance crossing to the level in steady state. 
The weight of the damper is limited to 5 percent blade gross weight, and the stroke is restricted to 
5 percent blade chord length due to the limited space in blade airfoil cavity. For the design of fluidlastic 
dampers, the tuning port area ratio (G) is one of the key parameters. From the view of practical 
application, a suitable G=40 is adopted in the baseline design. The parameters of the fluidlastic damper 
could be optimized and determined to satisfy the design objectives, which include tuning frequency, loss 
factor, prime mass, tuning mass and tuning port area ratio. The design flowchart is shown in Figure 2.36. 
From the general considerations, smaller damper mass, large G and the setup location of the damper at the 
blade tip are preferred. Some parameters of the damper are shown in Table 2.2. 


TABLE 2.2.— EMBEDDED CHORDWISE 
DAMPER PARAMETERS 


Prime mass, kg 1.0 

Tuning mass, kg 1.0 

Tuning port area ratio 40 

Damper radial position Blade tip 
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Figure 2.36. — Design flowchart. 


Figure 2.37 shows the rotating frequencies of the blade and damper, and the blade lag damping with 
respect to the rotor speed. The loss factor is 0.4. The blade lag damping increases initially then decreases 
with the rotor speed. The tuning frequency of the damper has significant influence on the blade damping. 
The peak blade damping decreases with the increase of the tuning frequency. Since the resonance occurs 
at the rotor speed 205 RPM, the damper is tuned to 7.2 Hz. When the damper is tuned to 7.2 Hz, the blade 
lag damping and dynamic stroke with respect to the rotor speed for different loss factors are shown in 
Figure 2.38. When loss factor is 0.3, the peak damping of the blade is more than 8 percent. The bandwidth 
performance is much better than other cases around the resonance rotor speed. Within 180-230 RPM, the 
blade critical damping is above 6 percent. The stroke of the damper decreases significantly with the 
increase of the loss factor. From the above analysis, the damper needs to be tuned to 7.2 Hz with loss 
factor 0.3. 
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Figure 2.37. — Blade damping and rotating frequencies versus rotor speed 
for different tuning frequencies (rj=0.4). 
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Figure 2.38. — Blade damping and the stroke versus rotor speed for 
different loss factors (co n =7.2 Hz). 
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2.3.5 Transient Loads Control Via Embedded Damper 

With the design parameters of the damper, the responses of the blade and damper can be calculated. 
The comparison between the transient response with the damper and the response without the damper can 
be utilized to evaluate the performance of the damper. It is necessary to check whether the displacement 
response of the damper is within the design limit or not. Due to the difference between the two models, 
the analysis of the influence of the tuning frequency and loss factor on the performance of the damper 
using the comprehensive code is needed to verify the design modeling. 

The prime mass, tuning mass and tuning port area ratio of the fluidlastic damper are as specified in 
Table 2.2. Since the blades go through the 2/rev resonance area, the primary excitation frequency in the 
damper modeling is 2.0 times the rotor speed. The moment ratio is defined as the ratio of the peak-peak 
lagwise root bending moment with the fluidlastic damper during 5 to 15 sec to the peak-peak moment 
without the fluidlastic damper. The stroke ratio is defined as the stroke of the damper divided by the blade 
chord length. 

For traditional absorbers, isolators or dampers, tuning frequency is one of the key design parameters. 
For different natural frequencies and loss factors, the moment ratios and stroke ratios are shown in Figure 
2.39 and Figure 2.40, respectively. With the increase of the loss factor, the moment ratio decreases firstly, 
then increases. In the range of the frequency from 5.6 to 8.0 Hz, the moment ratio also decreases firstly 
then increases. The stroke decreases significantly with the increase of loss factor. When the damper is 
tuned to 7.0 Hz with loss factor 0.2, the damper can reduce the peak-peak moment by 65.0 percent. The 
corresponding stroke is 0.72 percent blade chord length. When the damper is tuned to 6.9 Hz with loss 
factor 0.3, the peak-peak moment is reduced by 64.7 percent. The stroke decreases to 0.58 percent chord 
length. From the view of bandwidth, the damper with larger loss factor is obviously better than the one 
with smaller loss factor. With the variation of ±0.5 Hz tuning frequency or ±0.1 loss factor, about 
50 percent peak-peak moment can be reduced. In the following analysis, the damper is tuned to 6.9 Hz 
with loss factor 0.3. These design parameters are very close to the parameters calculated by the two 
degree-of-freedom model, which illustrates the design using the simple modeling is an accurate and 
effective method. 

The comparison of the time histories of the lagwise root bending moment and the stroke of the 
damper are shown in Figure 2.41. The steady loads in the first or last 5 sec increases compared with the 
case without the damper. The peak-peak moment with the damper during the resonance crossing increases 
by 7.61 percent the peak-peak moment at the first 5 sec. The stroke during the transient process is smaller 
than that in steady state. Figure 2.42 shows the comparison and the stroke, when the damper is tuned to 

6.5 Hz with loss factor 0.05. The steady response in the first 5 sec decreases significantly. The damper 
also attenuates the transient loads during the resonance crossing, and the peak-peak moment is reduced by 
27.0 percent. Obviously, the performance of the fluidlastic damper decreases significantly compared with 
the case with the tuning frequency 6.9 Hz and loss factor 0.3. The steady response in the last 5 sec has 
substantially small variation. From the above analysis, the side effect of the fluidlastic damper on the 
steady loads should be paid attention to. It can increase the steady loads, when it is tuned to its optimum 
value. The high transient loads usually exist for tens of seconds, but the blades have to undergo the steady 
loads for tens of hours for long endurance variable speed rotor helicopters. It is necessary to balance the 
benefits and side effects of fluidlastic dampers on transient and steady loads. 
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Figure 2.39. — Moment reduction ratio with tuning frequency and loss factor. 
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Figure 2.41 . — The comparison of moment response and the stroke of the 
damper (6.9 Hz, 0.3 loss factor). 



Figure 2.42. — The comparison of moment response and the stroke of the 
damper (6.5 Hz, 0.05 loss factor). 


2.3.6 Parametric Study 

For the optimization of the parameters of fluidlastic dampers, parametric studies are conducted in this 
section to enhance the performance. The parameters includes tuning port area ratios, loss factors, device 
mass, blade flapping, the duration during resonance crossing, rotor thrust level, and forward speed. 

Figure 2.43 shows the moment ratio, stroke and elastomeric stiffness versus tuning port area ratio. 
The damper is tuned to 6.9 Hz with loss factor 0.3 for all the cases. With the increase of the tuning port 
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area ratio, the variation of the moment ratio is substantially small. The stroke decreases by 90.8 percent, 
when tuning port area ratio changes from 10 to 100. The elastomeric stiffness increases significantly with 
tuning port area ratio. Large tuning port area ratio is preferred in the design of fluidlastic dampers. Figure 
2.44 shows the moment ratio, stroke and elastomeric stiffness versus prime mass. The gross weight of the 
prime mass and tuning mass is kept a constant value 2.0 kg. With the increase of the prime mass, the 
moment ratio and the stroke increase, and the elastomeric stiffness decreases. This means the reduction of 
the performance of the damper. When 0.5 kg prime mass is shifted as tuning mass, additional 23.8 percent 
peak-peak moment can be reduced. For the consideration of the gross weight of the damper, larger tuning 
mass and smaller prime mass are preferred. 

Blade flap motion has vital contribution to blade lagwise loads. For different longitudinal flapping, 
the time histories of the root bending moment with and without the fluidlastic damper are shown in Figure 
2.45. It is obvious that the transient loads increase significantly with the increase of blade flapping. The 
embedded chordwise damper can control the transient lagwise loads excellently during the resonance 
crossing. With 3.0° longitudinal flapping, the damper can reduce 65.7 percent peak-peak lagwise root 
bending moment. The corresponding time histories of the damper displacement are shown in Figure 2.46. 
The stroke of the damper increases significantly with the increase of the flapping. With 3.0° flapping, the 
stroke reaches 2.22 percent blade chord length. 





Figure 2.43. — Moment ratio, stroke and elastomeric stiffness versus G 
(co n =6.9 Hz, /7=0.3). 
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Figure 2.44. — Moment ratio, stroke and elastomeric stiffness versus 
prime mass (m p +mt=2.0 kg). 
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Figure 2.45. — Time responses of lagwise root bending moment for 
different flapping. 
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Figure 2.46. — Time responses of the displacement of the damper for 
different flapping. 
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Figure 2.47 shows the time histories of the lagwise root bending moment for different duration during 
the resonance crossing. With the increase of the duration, the damper can excellently control the loads. 
When the duration is 20 sec, the damper can reduce the peak-peak root bending moment by 66.3 percent. 
The corresponding time histories of the damper are shown in Figure 2.48. The stroke in transient is 
smaller than that in the first 5 sec. The maximum stroke is less than 1 percent blade chord length. 
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Figure 2.47. — Time responses of lagwise root bending moment for 
different duration. 
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Figure 2.48. — Time responses of the displacement of the damper for 
different duration. 
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Figure 2.49 shows the moment ratio and stroke versus forward speed. With the increase of the 
forward speed during the resonance crossing, the variation of the moment ratio reduction is less than 
5 percent. The stroke increases by 31.3 percent, when the forward speed increases from 34 to 136 km/h. 
Figure 2.50 shows the moment ratio and stroke versus rotor thrust. Within 60 percent variation of the 
rotor thrust, the variation of the moment ratio is less than 5 percent. The stroke increases significantly. 
With 30 percent increase of the rotor thrust, the stroke is less than 1 percent blade chord length. For 
different flight states, the variation of the performance of the damper is substantially small, and the stroke 
of the damper is limited to a low value (1 percent chord length). 
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Figure 2.49. — Moment ratio and stroke versus forward speed. 



0.7 0.8 0.9 1 1.1 1.2 1.3 

Thrust Ratio 



0.7 08 0.9 1 1.1 1.2 1.3 

Thrust Ratio 

Figure 2.50. — Moment ratio and stroke versus thrust. 
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2.3.7 Summary 

Embedded chordwise fluidlastic dampers are introduced to control rotor blade lagwise loads during 
resonance crossing for variable speed stiff in-plane rotors. The design of the fluidlastic dampers based on 
the two degree-of-freedom blade-damper system model illustrates that more than 6 percent critical 
damping around the resonance rotor speed can be provided to the blade with less than 5 percent blade 
mass penalty. The analysis based on the elastic blade model illustrates: 

• With the optimum combination of the tuning frequency and loss factor, the transient peak-peak 
lagwise root bending moment can be reduced by 64.7 percent, and the corresponding stroke is 
less than 0.6 percent blade chord length. Within the variation of ±0.5 Hz tuning frequency or ±0.1 
loss factor, about 50 percent peak-peak moment can be reduced. 

• Larger tuning port area ratio and tuning mass need to be utilized to enhance the performance of 
the damper, and control the stroke of the fluidlastic damper. When tuning port area ratio changes 
from 10 to 100, the stroke can be reduced by 90.8 percent. When 0.5 kg prime mass is shifted as 
tuning mass, additional 23.8 percent peak-peak moment can be reduced. 

• For different blade flapping, the damper can excellently control the transient lagwise loads. With 
3.0° longitudinal flapping, the damper can reduce 65.7 percent peak-peak lagwise root bending 
moment, and the stroke is less than 2.5 percent blade chord length. 

• For different period during the resonance crossing, the peak-peak transient loads can be reduced 
to the level in steady states. 

• At different forward speed or thrust, the variation of the performance of the damper is 
substantially small, and the stroke is limited to a low value. 

In a summary, the simple modeling is an accurate and effective method for the design of the 
fluidlastic dampers. The embedded chordwise fluidlastic damper can excellently control the lagwise 
transient loads during the lagwise resonance crossing for the variable speed stiff in-plane rotors. 

2.4 Rotor Modeling for System Integration 

Variable speed rotors change rotor speeds according to flight states to optimize power consumption 
(Ref. 50). For the integration of the rotor modeling with the transmission and engine system modeling, 
blade lag motion is necessarily to be taken into account to capture the transfer of lag dynamics to rotor 
shaft dynamics. Blade flap motion has vital contribution to the lag motion due to the centrifugal force, 
and needs to be included in the comprehensive model. Different from typical helicopter rotors, variable 
speed rotors can vary their rotor speed degree to reduce rotor torque, the terms introduced by the variation 
of rotor speed should be taken into account in the modeling to explore the dynamics associated with the 
variation of rotor speed. 

The modeling method derived by Jaw is utilized in the research (Ref. 63). The modeling can describe 
rigid blade dynamics with coupled flap and lag motions with the variation of rotor speed. The equations 
of motion for the coupled flap and lag motions are 

~ M exti = IxQpi + Ixatyt + fi) 2 c Pi s Pi + m i e y g R2 ^ 2s Pi c ^ ~ rrneygR^iSftSfr + m t gygRc Pl (52) 

+bthfi + ktf + mieygR^iCp.c^. + mieygR^Cp.s^. = -N ext . (53) 

The expressions for the externally applied aerodynamics on the blade are 

-M extl = C R (r-e)F z dr (54) 
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For the general description of variable speed rotors, the term associated with the lag spring is added to the 
equation for the lag motion. 

A nonlinear quasi-steady aerodynamic model is utilized, and the lift, drag, and moment coefficients of 
the aerofoil are calculated by a two-dimensional table-look-up method according to the angle of attack 
and the oncoming air flow (Mach number). The expressions for the lift and drag of the aerofoil are 


where 


F z = ±pV 2 cCi + pnb 2 (h — a h ba ) + pnbVa 

(56) 

F x =\pV 2 cC d 

(57) 

V = JOt + Up 

(58) 

a — 9 — arctan ^ 

Ut 

(59) 

U T = fir + V x sinx/j — V y cosx/j 

(60) 

Up — MIR + (r — e)/? + fl(y x cosip + V y simp ) 

(61) 


Up and U T are the resultant air velocity in the airfoil section plane. The direction of the local air flow 
around the airfoil is shown in Figure 2.51. 

The three state Pitt-Peters dynamic inflow model is utilized to determine the inflow distribution over 
the rotor disk during the steady and transient states (Ref. 46). The inflow with the variation of axial 
position and azimuth is 


A = A 0 + A c ^cosi/j + As^simp (62) 

The calculation of the inflow coefficients is shown in (Ref. 46). The direction of the inflow is shown in 
Figure 2.52. 
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Figure 2.52. — Inflow direction. 


The expression of the rotor torque is 


Qmr ~ Ql^mr 4" “f Q4^mr 4" Qs^mr^ ^6*0 (63) 

where 

Qi = my g R( 1 + 3ey 5 ) (64) 

<J 2 = i “(3y s +i) (65) 

<73 = 1 - 3yj (66) 

q 4 = meRq 3 (67) 

<? 5 = -2 my g R (68) 

<76 = ^3 y 3 +i) < 69 > 


A change in the torque transferred to the rotor system results in an angular variation of rotor speed 

Jh ~ &mr(Qer ~ Qmr ) (70) 

When switching to the transient response calculation of the rotor, initial steady response 
corresponding to specified rotor speed needs to be calculated. The steady response focuses on one blade 
modeling. Then the rotor extracts the initial values for the variables at different azimuth angle. With the 
initial condition, the transient response of the rotor can be conducted with prescribed rotor speed for 
single rotor analysis, as seen in Figure 2.53. For the exploration of the integrated system dynamics, the 
torque equation should be used to describe the coupling dynamics of the system. 

The implicit Newmark integration method (Ref. 47) is utilized to conduct the time integration over 
the duration of the variation of the rotor speed. For the solution of the nonlinear equations, the tangent 
stiffness, damping and mass matrices are utilized in the integration method. 
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Figure 2.53. — Response calculation flowchart. 


2.5 Simplified Rotor Model for Comprehensive Simulation 

The rotor model and analysis presented in the previous sections enable detailed evaluation of rotor 
aeromechanical phenomena during variable speed operation. This section presents a simplified rotor 
model which is utilized in the comprehensive, variable speed, propulsion system simulations presented in 
Section 7.0. This reduced order model includes rotor torque and thrust for hover or a tiltrotor in forward 
flight airplane mode with vehicle speed v c , Here a rigid rotor model based on blade element theory with 
quasi-steady aerodynamics (Ref. 41) is utilized (see Figure 2.54). Here the incremental rotor thrust dF T 
and torque dr at blade radial location r are expressed as 

dF T = N b ( dL cos (j) - dD sin (f) 
dr = N b r(dL sin (/) + dD cos $) 

where dL and dD are blade lift and drag force increments along with the number of blades N b . 
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The tangential velocity of a point on the blade at radial location r is 


U T - ^ MR r 


(72) 


based on this together with rotor axial speed, v c , and induced inflow velocity, v ind , the total inflow angle $ 
is expressed as 


where 


4 + 4ind 


(73a) 


(j) c = tan 1 


^ MR r 


and c/) ind = 


^ind 

A 2 + f W> 2 


(73b) 


here the induced inflow angle, 4) ind , is assumed small. Now, assuming uniform induced inflow velocity v ind 
together with uniform blade pitch, 6 pitc h , the following expression for total rotor thrust is obtained as 

Ft = fVfc +(j)ind -9 pitch )(v c </> ind -m MR )rdr 

2 J° ' (74) 

similarly the corresponding rotor torque is 


F = y c (rQ ;Vffl -vj ind )r 2 dr ... 

~ Pa "' a cL c F h Q. MR r * vj,^ + _ e pitch )(v c + iU MR (f> ind )r 2 dr 

Z 


both with 


'F. 



(75) 


(76) 


with air density p ain rotor radius, R b and blade chord c. Also c D is the blade section drag coefficient and 
a cL is blade section lift coefficient slope. Next, by equating the thrust expression in (74) with an 
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expression based on momentum theory (Ref. 41), we obtain the following equation for the induced inflow 
ratio A md 


with coefficients 


a X Knd + ^>xKnd + C X ~ 0 


(77a) 


ax 


= 48 In 


Yc L 

A, 


A c z + 6a cL <7 In 


l + Tc 

4 



(77b) 


= 24 sin 


-i 




A C J + 6a cL <7\n 


K 


A 2 + (- 24 Yc + 6a cL a tan 1 A c - 6a cL a6 pitch )A c - 3a cl cr (77c) 


( 


c x = a cL (y\ 


n-ie pUch -H 


\_±Yc 

K. 


A 2 - 2a cL y c cjXan '[A C ]A C 2 -a cL y c crA c 


— + 2a cL y c a{y c &p dc h tan A c j 

and with induced and forward speed inflow ratios 


7 _ ^ind 

ind ~ p n ’ 

K b ll MR 


*=- Vc 


RuCl 


and Yc = V 1 + V 


b^^MR 


and with rotor solidity ratio 


cNt 

nR h 


(77d) 


(78) 


(79) 


Based on Equations (77) to (79), the induced velocity, v ind , is determined from the roots of Equation (77a). 
Using this result in Equations (74) to (75), the rotor thrust and torque values can be computed for a given 
rotor speed, forward speed and blade pitch angle as shown in the rotor thrust and torque simulation block in 
Figure 2.55. The rotor parameters used in the comprehensive simulations presented in Section 7.0 are 
summarized in Table 2.3 and Table 2.4. 

Table 2.3 gives the main rotor parameters used in the single rotor helicopter simulation. 



Figure 2.55. — Rotor thrust and torque simulation block. 
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TABLE 2.3.— MAIN ROTOR PARAMETERS USED 
IN VARIABLE SPEED HELICOPTER SIMULATION 


Parameter Value 

Rotor speed (high), Q M R,high 190 RPM 

Rotor speed (low), Cl MR j ow 102.5 RPM 

Number of blades, N b 6 

Rotor radius, R b 11.5m 

Rotor solidity, a 0.115 

Blade mass, m b 386 kg 

Rotor rotational inertia, J MR 4.054x1 0 4 kg m 2 

Section lift coefficient, c D 0.012 

Section drag coefficient, a cL 5.08 


Table 2.4 gives the main rotor parameters used in the tiltrotor simulation. 

TABLE 2.4.— MAIN ROTOR PARAMETERS USED 
IN VARIABLE SPEED TILTROTOR SIMULATION 


Parameter Value 

Rotor speed (high), Cl MR high 190 RPM 

Rotor speed (low), Cl MR j ow 102.5 RPM 

Number of blades, N b 4 

Rotor radius, R b 9.91 m 

Rotor solidity, cr 0.13 

Blade mass, m b 372 kg 

Rotor rotational inertia, J MR 4.874xl0 4 kg m 2 

Section lift coefficient, c D 0.012 

Section drag coefficient, a cL 5.056 


2.6 Additional Information: Transformation Matrices Among Different Coordinate 
Frames 

cos(/? + /?q) 0 sin(/? + /?o) 

[Tfr] =0 1 0 

-sin(fi + j? 0 ) 0 cos(/3 + /? 0 ) 

'cos( —sin( O' 

[Tif] = sinZ sin$ 0 
0 0 1. 

TO 0 ' 

[ T P i] = 0 cos6 p sin9 P 
0 —sind p cosd p 

'—cosxp simp O' 

[T rs ] = simp cosip 0 
. 0 0 1. 

' cos/3 lc sinp lc sinp ls sin(3 lc cosp ls ' 

[■ ftps ] = 0 cosp ls -sinp ls 

—sin/3 lc cos/3 lc sinp! ls cos/3 lc cos/3 ls . 
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1 NJ 


v' 


w' 


[T] = 

—(v'cos9 1 + w’sin9 1 ) 

cos9 1 ^ 

v ,2 \ 

1 — — 1 — l/w'sin#! 

sin9 1 

i 1 - 

w' 2 \ 

~2~) 


v'sin9 1 — w'cos9 1 

— sin9 1 1 

fl — — v'w'cos9 1 
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3.0 Driveshaft Subsystem 

3.1 Introduction 

Rotorcraft power transmission shafting is mounted within flexible structures such as the wings, tail, 
or fuselage, depending on the rotorcraft application. These structures naturally will deform under normal 
operation of the aircraft. To accommodate such motion, the shafting system is comprised of several shaft 
segments connected by flexible couplings, which are generally nonconstant velocity (NCV). A shaft 
system with misaligned NCV couplings gives rise to a phenomenon called Cardan error in which the 
motion of the input shaft is not perfectly matched by the output shaft; for example, if the two shafts are 
misaligned and the input shaft has a constant angular velocity, the velocity of the output shaft will 
fluctuate around the input speed. This Cardan error leads to some interesting dynamic phenomena. 
DeSmidt, et al., 2002, (Ref. 27) for example, shows that a static misalignment angle in the NCV coupling 
creates a coupling between the lateral and torsional vibration of the shaft system, leading to self-excited 
oscillations and even instability in some conditions. Similar parametric instability phenomena is also 
demonstrated by Mazzei and Scott, 2006, (Ref. 64). 

3.2 Modeling 

After many preliminary studies, a model was formulated which captured the essential features of a 
rotorcraft driveline mounted on a flexible support. Figure 3.1 depicts this model, which is based on a 
tiltrotor wing with a shaft running through. 

Figure 3.2 shows a progression of modeling to help explain the final model. The left most figure 
shows a sketch of a typical tiltrotor aircraft, and we wish to model the wing structure and the cross- 
shafting, depicted in the middle pane. For clarity of analysis, the model is kept to a single shaft segment 
of length L coupled to the input shaft by a single NCV coupling, represented by a small circle between the 
two shafts. The torsionally flexible, laterally rigid shaft is connected to a lumped inertia, J 0 , at one end, 
representing a rotor load or the inertia of a shaft vibration mode. This shaft-inertia configuration, viewed 
as an inertia attached to a torsion bar with light torsional damping, has a torsional stiffness and damping 
given by k t and c h respectively. The wing structure is modeled as a flexible cantilevered beam 
(represented by the blue rectangle in the middle figure) and is connected to the shaft by a rigid bearing, so 
the shaft and structure essentially move together. If the structure primarily moves in its first vibration 
mode, then it can be assumed as an equivalent mass-spring-damper oscillator coupled to the shaft as 
shown in the right part of Figure 3.2. Here, the lateral inertia of the support structure, shafting, and rotor 
are all lumped together as the parameter 7 0 . The stiffness and damping of the structure are modeled as 

lumped elements k b and c b , respectively. T L (t) and F L (t) represent arbitrary torsional and lateral forcing, 
for example, from aerodynamic loading. 



Figure 3.1 . — Schematic of the model of a shaft system 
mounted on a flexible support. 


NAS A/CR— 20 13-21 6502 


53 



Figure 3.2. — Progression in modeling from a tiltrotor aircraft to the current model. Black boxes represent ground and 
the blue rectangle of the middle pane represents a flexible support structure. 


The lateral motion of the shaft pivoting about the NCV coupling is given by a{t ), which is a degree of 
freedom. The motion is assumed to be planar since a wing structure is generally less stiff in the vertical 
direction than in the fore-aft direction. The misalignment of the shaft is permitted to vary dynamically, 
which more realistically represents the system dynamics than a static misalignment angle. The input shaft 
has a prescribed motion, and as discussed before, the output motion of the NCV coupling does not 
match the input if misalignment is present. The angular motion directly after the NCV coupling is (j) x {a{t ), 
$)(7)). The motion of the inertia is then given by (j) x +(/) e (t), where (f) e (t) is the elastic motion of the inertia 
and is the second degree of freedom of the system. 

The motion (j) x can be determined using geometric arguments as in Martin, 1969, (Ref. 65). For planar 
misalignment, 


tan^ 0 

tan^ 


= cos a 


(i) 


Note that this is a kinematic constraint relating (j)\ and a, so it can be solved for (j) x and later it will replace 
all instances of (j) x in the equations of motion, leaving just the prescribed parameter and the degree of 
freedom a. 

Now the equations of motion can be derived by formulating the energies and applying Lagrange’s 
Equations. The kinetic and potential energies are given, respectively, by, 


T = \Jo(k+k) 2 +\ha 2 (2) 

V = \k£ + \k h (Laf (3) 


A Rayleigh dissipation function can also be written to account for damping, 

D = \c t $ +\c b (La ) 2 


( 4 ) 


Applying Lagrange’s Equations, 


d_(dJ ^ 
dtydq y 


dT dV dD „ 

+ + = Qa 

dq dq dq 


( 5 ) 


where q represents either (j) c or a and Q q represents the corresponding external loading, either T L and 
LF l , respectively, the equations of motion become 


4> e + 2(Q t % t <j) e + co}(j> e + '(j\=T L 


(6a) 


a + 2 co b % b a + o^a + ^ , <j> e ) = F L 


(6b) 
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Here, the following new parameters have been defined: 


co? = 


K_ 

j 0 


2 


,cot = 


k h B 




2Jo<*>t 


e = C b 0_ T = T J7 = LF L 

0 T id T T ' > ^ L T 

zi 0 a>b r 0 j 0 Y 0 


Also, the function g represents a complicated expression. 


(7) 


3.3 Analysis 

At this point, one can apply the constraint Equation (1) to the equations of motion (6). For the 
purposes of analysis, the terms in these nonlinear equations can be expanded in a Taylor series over the 
degree of freedom a around a = 0, the equilibrium position of the lateral motion of the shaft. The 
trigonometric functions are also converted to exponentials. 

Now, the method of multiple scales (MMS), Nayfeh and Mook, 1979 (Ref. 66), a nonlinear solution 
technique, can be applied to the equations of motion. MMS provides a way to find the approximate 
solution of a set of nonlinear equations by using multiple time scales as independent variables; 
mathematically: 


T 0 — t , 7| — s t, T? — . . . 


( 8 ) 


where s is a small parameter. Then time derivatives become ordered in the small parameter: 


d_ 

dt 


dT 0 d 
dt dT 0 


dT x d 

H — Dc\ + sD\ + . . . 

dt 8T X 


(9a) 


— = Di+2sD 0 D l+ ... (9b) 

dt 1 

Additionally, express the solutions </> e (t) and a(t) as straightforward expansions, neglecting terms of 0(e 3 ) 
and higher: 


</>e (t;e) = m {Tq,T x ) + £ 2 (p 2 (T 0 , T x ) 

(10a) 

cc(t',s ) = SOC\(To,T\) + £ 2 CC2(T(,,T x ) 

(10b) 

Also, the analysis will assume light damping of 0(f, '): 

ca t Zt =stu t ,CQ h % b = m, 

(11) 


Applying Equations (8) to (1 1) to Equations (6), considering the external forcing functions to be 0(s x ), 
and collecting terms of like orders of s (as is typical for MMS) gives the ordered equations of motion. 
Equations of 0(s l ) and 0(s 2 ) are presented here: 

0(e l ): 


Dq<P\ + co}cp x = T l 


(12a) 


£>o a x + cola x + 





= f l 


(12b) 
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(13a) 


Dq<p 2 + (o}(p 2 = -2D 0 D m - 2fj. t D^(p x - (er**> - e 2i *> )(«, (D 0 a, )^ 0 + 

- L ( e -2 to - e 2 uh ) (- a 2 <f>l + \ ( A>«i ) 2 + { «i (7> 0 2 «i )) 


Dfia 2 + colcc 2 + ^ - e 2i( h )^o ^2 = -2D 0 A^i ~~ 2 jU b D 0 a x - {e~ 2i< ^ - e 2i \x\(p\ ( 1 3b) 

Here, i represents the imaginary unit. At this point, the motion of the input shaft $) is still arbitrary. One 
reasonable choice for motion is a constant angular velocity described as: 


</>0 =£V, ^0=^0, </>0=° 

(14) 

Letting the external forcing go to zero, the solution to the 0(s l ) Equations (12) 

is: 

<Pi( T o, T \) = p \ (7) )e im,T ° + P\ (71 )e~ i( °‘ T o 

(15a) 

(r 0 ,T 1 ) = A 1 (T, )e ico >’ T o + A, (71 

(15b) 

Now Equations (14) and (15) can be used to fully define the right hand sides of the 0(s 2 ) equations, (13). 
After this substitution, it becomes clear that terms on the right hand side are secular for certain values of 
the input shaft speed Q 0 (rad/sec) in other words, these terms cause resonance in the 0(s 2 ) equations. This 
would mean that the higher order corrections to the solution grow quite large as time progresses, which 
makes the expansion not uniformly valid. The general procedure for dealing with these secular terms is to 
set them to zero, which provides conditions for finding the functions Pi and A x . 

When the shaft takes one of the following values, a secondary resonance occurs: 

D 

o 

II 

£ 

(16a) 

!0 

o 

II 

£ 

1 

£ 

(16b) 

Q 0 =\o) t + 6) h 

(16c) 


( 1 6d) 


Comparing with Nayfeh and Mook, 1979, (Ref. 66), these secondary resonance conditions appear to be 
cases for internal resonance, in which the two degrees of freedom become coupled nonlinearly, though 
linearly, they are uncoupled. 

3.4 Significance 

The conditions found in the previous section are the result of some interesting analysis, but it must 
also be shown that they have practical implications as well. The goal of this section is to find the natural 
frequencies of the driveline and foundation system in order to show that the shaft speed may coincide 
with one of the conditions in Equations (16). In this section, ball park values for the parameters of a more 
realistic shaft system will be estimated based on preliminary design data for the LCTR2 large civil 
tiltrotor concept (for example, Acree and Johnson, 2008, (Ref. 67)). 
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Based on internal discussions, it was decided to design a driveline for the scenario of a large tiltrotor 
with a total of four engines, one of them being disabled. This shafting should transmit 3750 HP, which is 
half of the power of an engine at full power. The shaft was designed to run at 7000 RPM, nearly the same 
as the cross-shafting of a CH-47, which is also a fairly large and powerful aircraft. The shaft uses 
aluminum material properties. 

The cases of torsional buckling and yielding were investigated as the two failure scenarios. These 
cases were evaluated using maximum shear stress theory and the buckling equations of Shokrieh, et al., 
2004, (Ref. 68). Total driveline length was estimated from Acree and Johnson, 2008, at 65 ft and the 
cases of 6 and 8 driveshaft segments were investigated. Using these constraints and reasonable estimates 
of dimensions and factors of safety (approximately 4), a potential shaft was sized with a wall thickness of 
0.1 18 in. and outer radius of 3.14 in. The stiffness can be calculated by using the equation of a rod in 
torsion, as in Rao, 2004, (Ref. 69), for instance: 


K = 


7lG 


32 L 


(d 4 


(17) 


DL 


where G is shear modulus, L dl is the driveline length, and D and d are the outer and inner shaft diameters, 
respectively. Then the stiffness is 12300 Nm. The shaft properties and dimensions are summarized in 
Table 3.1. 


TABLE 3.1.— SHAFT PROPERTIES AND DIMENSIONS 


Shaft property Value 

Power transmitted, P tr 3750 HP 

Design shaft speed, £l n 7000 RPM 

Material Aluminum 

Yield stress, S y 270 MPa 

Young’s modulus, E 70 GPa 

Shear modulus, G 27 GPa 

Mass density, p 2710 kg/m 3 

Driveline length, L dl 65 ft 

Wall thickness, t 0.118 in. 

Outer radius, r Q 3.14 in. 

Torsional stiffness, k t 12300 Nm 


Now the inertia properties of the rotor must be estimated. From Acree and Johnson, 2008, the mass of 
both rotor disks for the LCTR2 is 8756 lb. The individual rotor blade mass, accounting for the mass of the 
hub, is estimated as one tenth of this value, or 876 lb. Knowing the diameter of the rotor disk ( d rotor = 65 ft) 
from Acree and Johnson, 2008 and estimating the blade as a slender rod rotating about its end, one can 
compute the estimated rotor inertia as 52000 kg-m 2 . 

Due to the gearbox between the shafting and the rotor, the shafting will “see” a different inertia than 
this depending on the gear ratio. Acree and Johnson, 2008 give blade tip velocities for the cases of cruise 
and hover, v tc = 650 ft/s and v th = 400 ft/s, respectively. Using the blade length, one can determine the 
angular velocity of the rotor disk in cruise and hover, and knowing the shaft speed, the gear ratios for 
each case can be determined. The gearing makes the rotor inertia easier to turn from the point of view of 
the shaft, so the effective rotor inertias for hover and cruise, respectively, are J ef p h = 1420 kg-m 2 and 
Jeffx = 877 kg-m 2 . 

Using the estimated effective inertias and the shaft stiffness, now the torsional natural frequency of 
the driveline -rotor system can also be estimated assuming it behaves as a shaft pinned on one end and 
attached to an inertia (the rotor) at the other. For cruise and hover, these turn out to be co t>c = 3.75 rad/s 
and co t>h = 2.95 rad/s. 

From Acree and Johnson, 2008, the first lateral bending natural frequencies of the LCTR2 wing 
structure is about f b = 2.3 Hz or co b = 14.5 rad/s. Now, the shaft operates nominally at 7000 RPM = 735 
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rad/s. In view of Equations (16), the shaft speed will not coincide with a secondary resonance condition 
under a normal speed shift between hover and cruise. 

Considering higher modes of vibration could lead to such a coincidence though. The rotor inertia is 
quite large compared with the torsional inertia of the shaft, so treat the shaft as a pin-pin distributed 
parameter torsion rod. By Rao, 2004, (Ref. 69), natural frequencies of the rod are given by 


®t,n 


U7T 


LD 


-,n = 1,2,3,... 


(18) 


The first natural frequency is given by CD t ,i = 1001 rad/s. In view of Equations (16), during a 50 percent 
speed decrease, taking the shaft speed to 3500 RPM, or 368 rad/s, the shaft speed will cross secondary 
resonances at = 486 rad/s, Q 0 o = 500 rad/s , Q 0 c = 515 rad/s. So in this example, the first shaft 
torsional mode of vibration and the first lateral bending mode of the foundation have natural frequencies 
which allow Equations (16) to be satisfied, and these secondary resonances could be cause for concern 
during aircraft flight. These conditions will be investigated further in the next section. 


3.5 Constant Speed Operation Results 

Operating at these secondary resonances leads to phenomena which are not accounted for in linear 
analysis, so these scenarios could be important in the design of variable speed rotorcraft. If large 
vibrations are involved when the input shaft speed passes through one of these secondary resonances, then 
they should be considered and treated with care by designers, just as linear resonances play a role in 
aircraft design. It makes sense to first explore these conditions at constant speed operation to get an idea 
of how they behave before performing variable speed simulations. 

These secondary resonance conditions were found via some simplifications of the original nonlinear 
equations of motion (6), for example, using the Taylor expansion. It is prudent to run some simulations of 
the system using the full nonlinear NCV joint constraint, (1). For simulation, (1) is put into a slightly 
different form, solving for fa directly: 


fa = tan -1 


r tan fa ^ 
v cos a ; 


+ fa - tan -1 [tan ^ 0 ] 


(19) 


The first term on the right hand side provides for the fluctuation of the motion when misalignment is 
present, and the second and third terms on the right hand side guarantee that the angle fa steadily 
increases as fa steadily increases. 

Now Equations (6) and (19) can be numerically integrated to determine solutions of the system 
dynamics. To clearly showcase the phenomenology, convenient parameters will be chosen for initial 
simulations, taking care to keep the phenomena of these secondary resonance speeds from overlapping. In 
the following simulations, the external forcing will be set to zero, so these will be free responses of a 
parametrically excited system. This is done to avoid any confusion between effects produced by the 
oscillation of system and effects due to the external forcing. 

For comparison with subsequent examples, Figure 3.3 shows an example of the shaft system in free 
vibration with the input shaft speed set away from any of the conditions in Equations (16). The parameter 
set used for this simulation is co t = 100 rad/s, co b = 30 rad/s, Q 0 = 85 rad/s, rj = 0.5, Ct = Cb = 0, along with 
initial conditions y 0 = 0.1 rad, and a 0 = 0.01 rad in the torsional and lateral directions, respectively. The 
result is quite similar to two uncoupled, linear oscillators in free vibration. The left pane shows a very 
dense vibration signal, which is typical of subsequent simulations; for clarity, a zoom-in of the first 2 sec 
of the response is given in the right pane. 

The remainder of this section will present and discuss three examples of simulations under the 
conditions of (16). 
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Figure 3.3. — Nonresonant vibration of the shaft system with c o t = 100 
rad/s, u)b = 30 rad/s, Qo = 85 rad/s, ri = 0.5, £ = = 0, yo = 0.1 rad, 

and cro = 0.01 rad. The right pane zooms in on the first 2 sec of the 
response for clarity. 
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Figure 3.4. — System response (internal resonance) to the condition 
given in Equation (16d) with co t = 100 rad/s, c o b = 90 rad/s, Qo = 40 
rad/s, rj = 1 .2, Yo = 0.1 rad, and ao = 0.01 rad. 


3.5.1 Case 1 

Figure 3.4 shows an example of operating the system under the condition given by (16d). In this case, 
co t = 100 rad/s, co b = 90 rad/s, Q 0 = 40 rad/s, and rj = 1.2, and the system is given in initial conditions of 
y 0 = 0.1 rad and a 0 = 0.01 rad in the torsional and lateral directions, respectively. This would correspond 
to an initial twist in the shaft of about 6° and NCV joint misalignment of about 0.6°. Also, damping is set 
to zero. The response is an internal resonance characterized by strong coupling between coordinates that 
are linearly uncoupled. One can see that the larger amplitude torsional vibration drives the lateral 
vibration from a small value to a large one. 

Looking at the vibration energy in the system over time sheds some more light on this nonlinear 
coupling behavior. Mass specific energy in the torsional (E tors ) and lateral (Ei at ) directions, as well as the 
total specific energy (E tot ) in the system, can be computed by: 


E tors KEfors EE tors ^CC + ^ ^ 

(20a) 

Elat = KE lat +PE lat +\ CO b ( ! ) e 

(20b) 

Elat ~ E tors Ei a t 

(20c) 


These energies are plotted over time in Figure 3.5 for the same scenario as in Figure 3.4. As expected, 
the nonlinear coupling leads to an energy exchange between the torsional and lateral motion. An upshot 
of this is that if the frequency condition is satisfied, then the system could experience motion beyond what 
it was designed to handle, leading to increased fatigue or possibly to outright fracture. 

In addition, the total energy in the system is fluctuating over time. There is no external forcing on the 
system, only the variation in the parameters caused by the NCV coupling. More energy coming into the 
system could, again, lead to high vibration amplitudes. 

Changing the value of the inertia ratio rj also has effects on the system response. Increasing to rj = 3.3 
leads to an increase in the nonlinearity of the system. Figure 3.6 shows that the frequency of the energy 
exchange between modes also increases, as does the amplitude of the lateral vibration. Increasing further 
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to rj = 5.0, shown in Figure 3.7, causes the response to become chaotic. Also, the maximum amplitude of 
the lateral vibration has increased further. This chaotic vibration is qualitatively different than the internal 
resonance, but essentially, the potential harmful effects on the system (large vibrations) are the same. 

Adding a small amount of damping to the system can cause this motion to decay, at least in the case 
of free response under current investigation. Equal lateral and torsional damping of Ct = Cb = 0.01 is set in 
the system. If the system is linearized, this would correspond to 1 percent critical damping. 

The results of adding damping are shown in Figure 3.8. In the top pane, the damping effectively 
causes the motion to decay, preventing any large vibrations for the choice of system parameters and initial 
conditions defined for the previous figure. Increasing the initial conditions to y 0 = 0.2 rad, and a 0 = 0.1 
rad, as in the bottom pane of Figure 3.8, still allows for a substantial growth in the lateral motion of the 
shaft system before the damping takes over and causes the response to die away. Damping could be an 
effective way to deal with this type of response, though it should be chosen judiciously in order to ensure 
a limited vibration response. 

As a note, the condition Equation (16b) has behavior similar to (16d), so the investigation of (16b) 
need not be presented here. 



Figure 3.5. — Vibration energy in the shaft system with cot = 100 rad/s, 
cob = 90 rad/s, Qo = 40 rad/s, q = 1 .2, y 0 = 0.1 rad, and cro = 0.01 rad. 
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Figure 3.6. — System response (internal resonance transitioning 
to chaotic vibration) to the condition given in Equation (16d) 
with (jj t = 100 rad/s, c Ob = 90 rad/s, Qo = 40 rad/s, rj = 3.3, 

Yo = 0.1 rad, and ao = 0.01 rad. 
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Figure 3.7. — System response (chaotic vibration) to the condition 
given in Equation (16d) with co t = 100 rad/s, cob = 90 rad/s, 

Qo = 40 rad/s, rj = 5.0, yo = 0.1 rad, and ao = 0.01 rad. 
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Figure 3.8. — System response to the condition given in Equation (16d) 
with cot = 100 rad/s, c Ob = 90 rad/s, Qo = 40 rad/s, a] = 5.0, and £ = £> = 
0.01 . Initial conditions for the top pane are Yo = 0.1 rad and cro = 0.01 
rad, and for the bottom pane are yo = 0.2 rad and cro = 0.1 rad. 
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Figure 3.9. — System response to the condition given in Equation (16c) 
with u) t = 100 rad/s, c o b = 30 rad/s, Qo = 80 rad/s, ri = 0.5, yo = 0.1 rad, 
and cro = 0.01 rad. Top panel shows zero damping, while the bottom 
panel includes damping (£ = ( b = 0.02). 


3.5.2 Case 2 

The condition (16c) is investigated next. Values for the system parameters are set at co t = 100 rad/s, 
co b = 30 rad/s, Q 0 = 80 rad/s, and rj = 0.5. The initial conditions are the same as the previous example, 
letting y 0 = 0.1 rad and a 0 = 0.01 rad in the torsional and lateral directions, respectively. The result for 
zero damping is given in the top panel of Figure 3.9, which shows a rather severe instability. Clearly, this 
is a worrisome situation during shaft operation. Adding some damping (C t = Cb = 0.02, or 2 percent critical 
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of the linearized system) to the system counteracts this instability quite effectively for this parameter set, 
as shown in the bottom panel of Figure 3.9. 

Increasing the inertia ratio requires an increase in the damping level to counteract the instability. In 
Figure 3.10, the inertia ratio has been increased to rj = 0.5. In the top panel, the lateral damping is C b = 
0.03, and this level of damping keeps the response from growing unbounded, but the response does reach 
a fairly high steady state amplitude of about 0.4 rad or 23°. The bottom panel shows that a much higher 
amount of lateral damping, C b = 0.12, is necessary to cause the vibration to decay. These results show that 
large amplitude can be realized at this shaft speed, but that an appropriate choice of damping may be able 
to limit the vibrations. 
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Figure 3.10. — System response to the condition given in Equation (16c) 
with u) t = 100 rad/s, u) b = 30 rad/s, Qo = 80 rad/s, q = 0.5, £ = 0.02, 

(b = 0, Yo = 0.1 rad, and ao = 0.01 rad. In the top panel, ( b = 0.03, 
while in the bottom panel, = 0.12. 
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3.5.3 Case 3 


Now condition (16a) is investigated with the parameter choices co t = 100 rad/s, co b = 30 rad/s, Q 0 = 50 
rad/s, and r\ = 0.5. In Figure 3.11, applying zero damping and initial conditions y 0 = 0-1 rad, and a 0 = 0.01 
rad, it appears that nothing interesting is happening here, and even that the motion in the torsional 
direction may be decaying. 

Increasing the lateral initial condition to a 0 = 0.03 rad reveals that the torsional response at this shaft 
speed first decreases and then grows unstable, as displayed in the top pane of Figure 3.12. The bottom 
pane shows that decreasing the torsional initial condition to y 0 = 0.01 rad actually results in a larger 
response sooner in time. It is also noted that these adjustments have negligible effects on the lateral 
response. Also, many simulations (not shown here for brevity) indicate that the inertia ratio r\ has 
essentially no effect on the response of the system. As can be observed in the equations of motion (6), rj 
only effects the lateral dynamics. Since the instability occurs in the torsional dynamics, this result makes 
sense. 

Adding a small amount of damping (Q = £ b = 0.01) to the system stabilizes the system quite readily, 
even for somewhat large initial conditions (y 0 = 0.1 rad, and a 0 = 0-1 rad), as presented in Figure 3.13. 
Other choices for initial conditions yield the same qualitative behavior. Simulations indicate that for 
realistic values of damping and vibrations, the condition (16a) may not present problems for shaft 
operation, though behavior for zero damping is noted. 
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Figure 3.11. — System response to the condition given in Equation (16a) with 
cot = 100 rad/s, cj b = 30 rad/s, Qo = 50 rad/s, r] = 0.5, ( t = (b = 0, yo = 0.1 rad, 
and ao = 0.01 rad. 
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Figure 3.12. — System response to the condition given in (Eq. (16a)) 
with u) t = 100 rad/s, c Ub = 30 rad/s, Qo = 50 rad/s, rj = 0.5, £ = £> = 0, 
and cro = 0.03 rad. In the top pane, Yo = 0.1 rad and in the bottom 
pane, yo = 0.01 rad. 
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Figure 3.13. — System response to the condition given in Equation (16a) 
with u) t = 100 rad/s, c o b = 30 rad/s, Qo = 50 rad/s, rj = 0.5, $ = = 0.01 , 

Yo = 0.1 rad, and ao = 0.1 rad. 


3.6 Variable Speed Operation Results 

Now that the secondary resonance conditions of Equations (16) have been investigated at constant 
shaft speed and their chief behavior has been cataloged, the effects of variable speed shaft operation, 
accelerating through these troublesome shafts speeds, can be simulated and studied. First, as a base for 
comparison, the system will be accelerated through a shaft speed region where there is no resonance. This 
result is given in Figure 3.14, taking co t = 100 rad/s, co b = 90 rad/s, and r\ = 1.2, and the system is given in 
initial conditions of y 0 = 0.1 rad and a 0 = 0.01 rad in the torsional and lateral directions, respectively. The 
shaft speed profile is shown in the top pane, starting at an initial speed of Q, = 70 rad/s for 3 sec, 
accelerating to a final speed of Q/= 80 rad/s over 10 sec, and then holding the final speed for an 
additional 3 sec. The initial region of constant speed operation allows the system to settle before the speed 
change. The final region of constant speed allows the system to settle from any effects of the speed 
change. This result here is essentially a free response of two uncoupled oscillators, which is expected. 

Next, the parameters corresponding to Example 1 are implemented, in particular, co t = 100 rad/s, co b = 
90 rad/s, and 77 = 1.2, and the system is given in initial conditions of y 0 = 0.1 rad and a 0 = 0-01 rad. The 
system will be accelerated through the condition Equation (16d), or Q 0 = 40 rad/s for this parameter set, 
using a shaft speed schedule as defined by the top pane of Figure 3.15. The result is that the lateral 
amplitude of the shaft becomes significantly larger as the shaft speed passes through the troublesome 
speed at a simulation time of 8 sec. This transient growth occurs because the shaft system transitions from 
a speed regime where the coupling between the lateral and torsional dynamics is rather weak to a regime 
where the coupling is strong and the lateral dynamics are excited. Such a significant amplitude increase 
could be cause for concern as an aircraft changes speeds in flight. 

This scenario will be simulated again for a slower acceleration to examine the effects of dwelling near 
the secondary resonance for a longer time. In the previous case, the shift time was 10 sec; Figure 3.16 
gives the result of shifting over 20 sec. Note the longer time of simulation on the horizontal axis. 
Curiously, in spite of the longer dwell time, the amplitude decreases in the lateral response. 
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Figure 3.14. — System response to a variable speed input with co t = 100 rad/s, 
c Ob = 90 rad/s, rj = 1 .2, yo = 0.1 rad, ao = 0.01 rad, Q, = 70 rad/s, and Of = 80 rad/s. 
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Figure 3.15. — Accelerate the shaft through the speed condition Equation (16d) 
in 10 sec. Parameters are co t = 100 rad/s, c Ob = 90 rad/s, ri = 1 .2, yo = 0.1 rad, 
ao = 0.01 rad, O, = 35 rad/s, and Of = 45 rad/s. 
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Figure 3.16. — Accelerate the shaft through the speed condition Equation (16d) 
in 20 sec. Parameters are co t = 100 rad/s, c Ub = 90 rad/s, 77 = 1 .2, yo = 0.1 rad, 
ao = 0.01 rad, O, = 35 rad/s, and Of = 45 rad/s. 

This can be explained by another simulation given in Figure 3.17. The only change in this simulation 
is to extend the time of constant speed operation at the beginning of the simulation by 1 sec from 3 sec in 
Figure 3.16 to 4 sec in Figure 3.17. As the system transitions through this secondary resonance regime, a 
big factor in the amplitude of the motion is the “initial condition” of the response as the shaft speed 
approaches a secondary resonance condition. If, in the course of oscillation, the motion is nearly zero as 
the condition is approached, then the amplitude will naturally be small as in Figure 3.16, since in a 
nonlinear system, small amplitudes of vibrations don’t engage the nonlinear effects very much. 
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On the other hand, if the oscillation is relatively large as the shaft speed approaches the secondary 
resonance condition, the motion has the potential to grow large, as in Figure 3.17. In addition, the lateral 
response is larger than in Figure 3.15, confirming that a longer dwell time also can lead to a larger response. 

Accelerating the shaft through the secondary resonance condition given by Equation (16a) does not 
lead to any notable results. Responses typically look like those in Figure 3.14. This instability grows 
much more slowly in amplitude than the growth produced by the other conditions, so it follows that the 
response would not grow much by passing through this speed. 
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Figure 3.17. — Accelerate the shaft through the speed condition Equation (16d) in 
20 sec. Parameters are u) t = 100 rad/s, cob = 90 rad/s, q = 1 .2, y 0 = 0.1 rad, 
cro = 0.01 rad, Q, = 35 rad/s, and Q f= 45 rad/s. The only change from the results 
of Figure 3.16 is to extend the initial constant speed of constant speed operation 
by 1 sec. 
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Accelerating the shaft through the other conditions in Equations (16) give qualitatively similar results; 
an amplitude change which depends on oscillation amplitude as the shaft speed approaches the secondary 
resonance and also a dependence on dwell time near the speed of secondary resonance. 

Effects of damping on variable speed operation were not presented here. Since these simulations are 
free vibration, the damping would only lead to a decrease in amplitude as the shaft speed approaches a 
secondary resonance. At that point, it is equivalent to simply starting the system at constant, but lower, 
amplitude. 

3.7 Conclusions 

Due to the secondary resonances brought on by the nonlinearity of the flexible couplings, large 
vibrations could arise in the shaft system. Since in variable speed operation, the shaft may undergo speed 
changes over a wide range, it makes sense to consider this nonlinearity in the design of variable speed 
aircraft. The nonlinear modeling and analysis allows one to find these shaft speeds that may cause 
problems. Selecting an appropriate amount of damping and choosing a proper speed schedule for 
performing a rotor speed shift are two possible strategies for minimizing shaft vibration levels. 

To further this investigation, more in-depth nonlinear analysis possibly could be performed to support 
the existence of internal resonance, chaotic vibration, and instability. More analysis may lead to closed 
form or semi-closed form ways to determine what parameter values cause transition from internal 
resonance to chaotic vibration and what values of damping are needed to subdue effects of operating near 
a secondary resonance condition. More realistic modeling, such as implementing a distributed parameter 
flexible foundation or using multiple shafts and NCV couplings, would add fidelity to the results. Free 
vibration of the system, as is studied here, is great for an initial study, but is not very realistic in terms of 
an aircraft operating in the real world; realistic external loading should be incorporated. Adding some 
unbalance to the shaft or rotor would also generate some interesting results, as one could compare the 
effects of operating near an ordinary linear resonance to operating near a nonlinear secondary resonance. 
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4.0 Gas Turbine Engine Modeling 

4.1 Introduction 

In order to simulate the variable speed operation of the rotorcraft driveline, a transient gas turbine 
engine model has been developed. In general the gas turbine will undergo transient operations due to 
startups, changing loads, and sudden shutdowns. Here the main focus is to capture the effects due to 
sudden load changes which will be encountered during gear changes. 

In the area of transient gas turbine engine modeling several researchers have explored techniques 
based on linearization models about the nominal operating conditions. Here the model Jacobian is 
estimated and the engine response about the nominal is computed (Refs. 70 and 71). However, such 
approaches are not valid for significant deviations about nominal which could be expected during severe 
or sudden operational load changes. More recently engine models based on the full nonlinear aero- 
thermal relations have been developed (Refs. 72 to 78). In these models, the nonlinear aero-thermal 
equations have been solved using various Newton-Raphson based methods. Such iterative solving 
procedures have been further explored by (Ref. 79) where it was found that the iterative nature was 
problematic for real-time simulations due to slow converges times. To address this, others have proposed 
non-iterative simulations techniques based on forward substitution methods (Ref. 80). Still others have 
utilized simplified compressor and turbine models which do not rely on engine maps allowing explicit 
solutions without the need for iteration at the expense of reduced simulation fidelity (Ref. 81). 

The gas turbine engine model developed in this investigation is a nonlinear mechanical-aero-thermal 
model with lumped representation of the gas turbine components. The solution approach is based on non- 
iterative forward substitution similar to (Ref. 80) together with a new analytical method to explicitly 
compute the compressor and turbine performance maps. One key feature for a transient model is that the 
mass flow through each component is not necessarily the same. For example, at any instant the 
compressor and turbine could induce different mass flows depending on the instantaneous upstream and 
downstream pressures across each component. To account for this plenum chamber elements are included 
between each component to accommodate these transient differential mass flows and density changes. 

The mechanical-aero-thermal equations are implemented in the Matlab Simulink environment for use as a 
module in the comprehensive variable speed drive system simulation. 

4.2 Transient Gas Turbine Engine Model 

To simulate the effect of the engine and fuel control system dynamics during variable speed rotor 
operation, a generic twin spool gas turbine model is developed. The components of this model are shown 
in Figure 4.1. The components included in this analytical model are axial compressor, combustor, gas- 
generator turbine, power turbine and two mass accumulators (plenums). The torque balances between the 
compressor and gas generator turbine as well as between the power turbine and the load have also been 
included and the physical and thermodynamic laws (Ref. 75) have been utilized to describe the system 
dynamics. This analytical model is then implemented in the MATLAB/Simulink environment. The 
thermodynamic model is based on the standard idealized Joule/Brayton cycle (Ref. 70). The main 
assumptions made in this analytical model can be summarized as follows: 
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Figure 4.1. — Two spool gas turbine engine system components and analysis stations. 


Engine Model Assumptions: 

1 . The air mass flows through the system are all one-dimensional at the mean line. 

2. Air with ideal gas law behavior is assumed. 

3. All expansion and compression is adiabatic. 

4. The inlet stagnation temperature and pressures are constant. 

5. All power transfer occurs in the rotor stages, the stators do no work. 

6. The instantaneous compressor mass flow, pressure ratio and efficiency are determined via pre- 
calculated beta maps. 

7. Gas generator and power turbines are in a state of choking flow 

8. No air bleed occurs from the compressor. 

9. The combustion is modeled as isobaric heat addition. 

The power required to drive the compressor is 


W c =m c Cp(T in -T (tl ) 


( 1 ) 


Where m c is the compressor air mass flow rate, c p is the specific heat of air and r 0 1 and T 02 are the 

stagnation temperatures at the compressor entrance and exit respectively. Furthermore the stagnation 
temperature ratio is related to the compressor pressure ratio by 


?o i tic 


f p \ 
r 02 


y - 1 


v^oiy 


( 2 ) 


Where r/ c is the compressor isentropic efficiency, y is the specific heat ratio, and P 0 i and P 02 are the 
stagnation temperatures at the compressor entrance and exit respectively. Based on (1) and (2), the torque 
r c required to drive the compressor is 
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Where Q c is the compressor rotation speed and J c is the compressor rotational inertia. Also, in the above 
equations the stagnation temperature and pressures at the compressor entrance are assumed to be the 
ambient atmospheric pressure, P 0 i=P a , and temperature T 0 i = T a values. Furthermore compressor pressure 
ratio, n c - P<)2 ! P<3\ , and efficiency, r/ c , are related to m c and Q c the via the compressor characteristic 
maps, 


ft c 3-iid tj c f cr j (7qi , /qi , tH c , F2 C ) (4) 

In this study, the compressor pressure and efficiency maps are generated numerically based on the 
nominal compressor design values as detailed in Section 4.3. 

Furthermore, the compressor mass flow rate, m c , pressure ratio, 7r c , and efficiency characteristics, r / c , 
are computed by way of so-called beta maps with functional dependencies, 


n c = fpr(T a ,P a ,£l c ,P c ) 

rile — fm (Pa Pc ) 


(5a) 


tlc=MT a ,P a ,n c ,P c ) 


Where J3 C is the compressor beta parameter which maps the mass flow to pressure ratio as 

7i c = p c m c 2 (5b) 

In this study, the compressor pressure and efficiency maps are generated numerically based on the 
nominal compressor design values (see Section 4.5). 

To drive the compressor and load, separate gas generator and power turbines are present in the two- 
spool gas turbine engine (Figure 4.1). The power generated from these turbines is 

w ggt = riiggtC p (T 04 - T 05 ) and W pt = m pt c p (T 06 - T m ) (6) 


Where m ggt and m pt are gas generator and power turbine mass flows, respectively, and T 04 , T 05 , T 06 and 

T 0 7 are turbine entrance and exit stagnation temperatures as defined in Figure 4.1. Furthermore, the 
turbine pressure ratios based on isentropic expansion thought the turbines are 
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and the turbine mass flows are computed based on the assumption that both turbines are operating in the 
condition of choking flow. Therefore 


riiggt=<P 


A 4 P 04 



and m pt - cp 


^6^06 
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(8a) 


with 
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(8b) 


<p = 
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Where R is the ideal gas constant. Next, since the gas generator turbine and compressor are on a common 
spool, we have Q c = Cl ggt , Also, the torque difference between the compressor and gas generator turbine is 

Uc + Jggt)^c = -^(W ggt - w c ) (9a) 

and the torque difference between power turbine and load is 

Up, + J L )& pt = -^-(Wp - W L ) (9b) 


Where W L is the power consumed by the load and J c , J ggh J pt and J L are rotational inertias of the 
compressor, gas generator turbine, power turbine, and the load, respectively. 

During steady state operation the mass flow rates through the compressor and turbines are equal, 
m c = m ggt = m pt . However, during transient operation, differential mass flows can develop between the 

components due to volume packing effects. To account for these transient effects, two plenum chambers 
(Plenum 1, between compressor and combustor, and Plenum 2 between the turbine units) are included in 
the model as depicted in Figure 4.1. 

The plenum chambers are modeled as ideal mass accumulators with no losses or heat transfer. 

Figure 4.2 shows a schematic of the plenum chamber element with input and output mass flow rates m in 
and m out generated by the upstream and downstream compressor and turbine subsystems. Due to the 
differential mass flows, the rate of change of air density inside a plenum, p p , is given by 


P P 


Win ~ Win 


( 10 ) 


Where V p is the plenum chamber volume. Furthermore, through differentiation of the ideal gas law and 
assuming adiabatic compression, the following differential equation for plenum chamber pressure P p is 
obtained 


PpU) = r(m '" m ' M,) P p (t) with M (t) - [ (m, n - m out )dt + M p init (1 1) 

M p (t ) J 

Here M p is the mass of the air accumulating inside the plenum and M p init is the initial air mass. 



Figure 4.2. — Plenum chamber element. 
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It is assumed that the plenum pressure is uniform such that input and output pressures P in = P p and 
P oul = P p equal the internal plenum pressure. Furthermore, based on energy conservation inside the 
plenum we have 


d^pMpTpp) 

dt 


Cpdii n Toin Cpdi ou tTo out 


Where T 0in , and T 0out are entrance and exit stagnation temperatures and T 0p is the internal plenum 
stagnation temperature. Using (10)-(12) the exit stagnation temperate can be approximated as 


( 12 ) 


L 0 out 


m n 


rn in T 0in P p 


Yl 

R 


(13) 


Finally, the combustion chamber (entrance station 3, exit station 4 in Figure 4.1) is modeled as a pure 
isobaric heat addition process. With heat addition, Q cc , based on an idealized combustion model 


Qcc = m f LHV 


(14) 


Where rhf is the fuel mass flow rate and LHV is the fuel lower hearting value. 

As a result of an energy balance and the input stagnation temperature, T 0 3 , the combustion chamber 
exit stagnation temperature T 04 is determined as a solution to 

Tcc ^dt L= niccCp{Tm ~ t ° 4) +rh f LHV ( ,5 ) 


Where r cc is the combustion chamber time-constant and m cc is the air mass flow through the combustion 
chamber which equals the gas generator mass flow m cc = m ggt . Furthermore, the combustion chamber exit 
stagnation pressure, P 0 4 , is computed from the input stagnation pressure P 03 using, 

^04 =r lccPo3 (16) 

Where ij cc < 1 represents the combustion chamber pressure drop. Also, at steady-state the combustion 
temperature rise AT 0cc = ^04 “ 7q 3 is 


AT r LHV m /• 

Ar 0cc = / with / = — 


m r 


Where /is the fuel-air ratio. 


Qcc = MfLHV w cc =0 



1 04 


04 
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Figure 4.3. — Combustion chamber. 


(17) 
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4.3 Engine Nominal Design 

Next the overall engine design is determined by the choice of several nominal operating parameters. 
These parameters, selected based on the NASA LCTR-2 heavy lift concept (Refs. 82 and 83) are 
summarized in Table 4.1. Here the subscript refers to nominal conditions. 


TABLE 4. 1 .—ENGINE NOMINAL DESIGN PARAMETERS 


Parameter Value 

Ambient temperature, T a 288° K 

Ambient pressure, P a 1.023 bar 

Nominal output (load) power, w L 7500 Hp 

Nominal compressor mass flow, m cn 13.5 kg/s 

Nominal compressor pressure ratio, pr cn 28 

Nominal compressor axial velocity, C an 150 m/s 

Number of compressor stages, N s 8 

Compressor mean radius, R c 0.119m 

Compressor inlet hub-to-tip ratio, ht c 0.35 

Nominal compressor speed, Cl cn 26,500 RPM 

Power turbine mean radius, R pt 0.1811m 

Power turbine inlet hub-to-tip ratio, ht pt 0.65 

Nominal power turbine speed, Cl ptn 12,500 RPM 

Compressor rotational inertia, J c 0.641 kg m 2 

Gas generator turbine rotational inertia, J ggt 0.144 kg m 2 

Power turbine rotational inertia, J pt 1 .04 kg m 2 

Plenum chamber volume, V p 0.1 m 3 

Combustion chamber loss coefficient, ij cc 0.85 

Combustion chamber time constant, r cc 5.0><10' 4 sec 

Fuel lower heating value, LHV 43.1 MJ/kg 

Nominal fuel rate, m fn 0.225 kg/sec 


Assuming 100 percent isentropic compression efficiency at nominal conditions ( rj cn = 1) we have 


Y - 1 

^ 02 n — TaP^cn Y 


and A T 0cn = T r 


y- 1 

pr cn Y -1 


(18) 


Where A T 0cn = T 02n - T a is the nominal stagnation temperature rise in the compressor. Under steady-state 
conditions we have 


W ggt =W c and W pt =W L (19) 

Thus using (1), (5), (16), (18) and (19) the required combustion chamber temperature rise at nominal 
operating conditions is 
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And the nominal gas generator and power turbine temperature drops are 


AT Oggta — T a 


y- i A 

pr cn r “I 


and AT 


0 ptn 


w Ln 




Finally, the nominal fuel flow rate and fuel/air ratio are 

m fn = 


^cn^pAT Q ccn j' — 


LHV 


m n 


( 21 ) 


( 22 ) 


4.4 Compressor and Turbine Design 

The compressor design is based on a multistage constant mean radius design. Figure 4.4 shows the 
eight stage (N c = 8) axial compressor modeled in this study which also includes an inlet guide vane (IGV) 
stage and an exit guide vane (EGV) stage. 

Figure 4.5 shows a typical rotor/stator stage and the corresponding air flow velocity triangles. Here 
U c is the mean line rotor speed of the compressor 

U c =R m Q c (23) 

Where R m is the compressor mean radius. As shown in Figure 4.5, C u C 2 and C 3 are the absolute air 
velocities and a l9 a 2 and a 3 are the corresponding flow angles. Also, V u V 2 and V 3 and P 2 and p 3 are 
the flow velocities and flow angles relative to the rotor. Here, the stage axial flow velocities C a i, C a2 and 
C a3 are 


C a \ — G] cos CC\ 

C a2 = C 2 cosa 2 (24) 

G «3 - C 3 cos cc 3 


Here the entrance flow angle is equal to the upstream stage exit angle. Based on the flow velocity 
triangles, the flow angles also satisfy 


— — = tan a x + tan p x 

C a \ 


Cal 


= tma 2 +tanp 2 


and 


— — = tan^ 3 + tan p 3 

0*3 


(25) 
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Furthermore, the rotor entrance, rotor exit and stator exit stagnation temperatures T 01s . T$ 2s and T^ s are 

c ^ c ^ c ^ 

Ton = Ty s + —2— , T 02s = T 2s + — — and T 03s = T 3s + —2— (26) 

p AC p AC p 

with corresponding flow temperatures T Xs , T 2s and T 3s . Also, the compressor mass flow rate, m c , which is 
equivalent at every stage, satisfies 

^ c — P\A\sC a \ — Pl^-isC al ~ P3^-3s^a3 (^ 7 ) 

Where p x , p 2 and p^ are the flow densities and^, A 2s and^ are stage cross-section areas. In each stage, 
the rotor increases the energy of the flow by increasing the swirl velocity. The power added in each 
compressor stage is 

W cs = m c c p AT 0s = m c U c ( C a2 tan a 2 - C al tan a x ) (28) 


Where A T 0s = T 02s - T 0Xs is the stage stagnation temperature rise. After each rotor the swirl velocity is 
reduced by the stators. However, the stators do not add energy to the flow thus ^03 5 - T()2s . Likewise, 
neither does the IGV or EGV stages. 

Based on the total compressor stagnation temperature rise A T 0cn in (18), the nominal stage stagnation 
temperature rise, A T 0sn , for each stage is designed to be identical, therefore 

AT 

Ar 0OT =— fL (29a) 

and the stagnation temperature at the k th stage is 

T(ft sn ,k = T u + (k - \)AT 0sn and T 02sn k = T 03sn k =T a + kAT 0sn (29b) 
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The compressor stage cross-sectional areas A ls , A 2s and A 3s are determined such that the axial flow 
velocities though each stage remain constant at the nominal operating condition. Thus, at nominal 
conditions we have 


Can C a \ C a 2 C a 3 . 


( 31 ) 


In this case, the nominal design axial flow speed is selected as C an =150 m/s. Using (28)-(30), the 
nominal flow angles for the N c identical stages are 


a \n ~ a IGV j 


a 2n = tan - 1 


A T 0sn 
U C 

w cn^an 


+ tan a Ul 


a 3 n ~ a \n ~ a IGV 


(32) 


Where a IGV is the IGV exit angle which is used as the free design variable (within limits). Note, for 
identical stages, since the stator exit angle is the same as rotor entrance angle of the next stage, we have 
a 3 n = oc\ n . Also, using (25) together with (3 1), fan and jB 3n can be determined. Here the range of 
allowable a IGV \ s determined by checking the so called de Haller number (Ref. 70). Next, using, (26), 
(27), (29) and (30), together with the ideal gas law expression for density and the adiabatic compression 
relations starting from the compressor inlet, the compressor cross-section areas at the k th stage are 


A\s k 


r 

Mcn RT ar-\ 

l ? ^2 s,k 

C an R> aJ'\sn,k 


^cn RrR aY~^ 


^cn RR a Y~\ 


A lv t and A 3s9k = " [k= 1 - N c ] (33) 


Can R a R 2sn,kY 1 


CanPa R 3sn,kY -1 


Where R is the ideal gas constant. Finally, the compressor mean radius is 


| As,i i + 

V 471 l ~J U ht 


(34) 


Where ju ht = R tip /R hub is the hub-to-tip ratio at the inlet, which is a geometric design variable. Here the hub- 
to-tip ratio is selected as ju ht = 0.35. The resulting compressor stage hub and tip radii for this design are 
show in Figure 4.6, and the corresponding compressor stage angles are summarized in Table 4.2. 
Furthermore, the gas generator and power turbine angles as well as the rotational inertias computed based 
on the component mean radius values are also given in Table 4.2. 
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TABLE 4.2.— COMPRESSOR AND TURBINE DESIGN PARAMETERS 



Temperature rise or 
drop 

Absolute 

angles 

Relative 

angles 

Mean 

radius 

Rotational 

inertia 

Compressor 
(8 Stage) 

Ar 0c „ = 458.2 1°K 

cc\ n — 40.0° 
a 2n = 63.43° 
a 3n = 40.0° 

P i„= 53.79° 

Pm =11-60° 

A»= 53.79° 

R mc = 0.12 m 

J c =1.54 kg m 2 

Gas Generator Turbine 
(2 Stage) 

AT 0ggm = 458.21 °K 

a ]n = 10.0° 
a 2n = 58.38° 
a 3n = 10.0° 

Pu= 54.96° 
Pm= 20.49° 
P in = 54.96° 

Rmg g t = 0.13 m 

J ggt = 0.26 kg m 2 

Power Turbine 
(2 Stage) 

AT 0ptn = 41 2.22 °K 

a ]n = 10.0° 
a 2n = 53.38° 
a 3n = 10.0° 

Pln= 51.17° 
p2n= 24.23° 
p3n = 51.17° 

R mpt = 0.24 m 

J pt = 0.98 kg m 2 



Figure 4.7. — Assumed flow angle deviations for off-nominal operation. 


4.5 Analytical Calculation of Off-Nominal Characteristic Maps 

Under transient conditions such as during gear ratio changes or during sudden main rotor power 
demands, the engine will undergo deviations about the nominal design operating point which is regulated 
by the engine fuel control. Therefore, it is necessary to model the off-nominal operation of compressor 
and turbine units. The off-nominal characteristics are typically described by compressor and turbine maps 
which show pressure ratio and isentropic efficiency as a function of mass flow for a given rotational 
speed as shown in (4). The compressor map can either be obtained experimentally or estimated 
analytically. In this investigation, compressor map data was not available, therefore a novel procedure is 
developed to analytically generate these maps for a given compressor design. 

Since the stage areas and blade angles are fixed, the constant axial velocity condition (31) is only 
valid at the nominal operating point. For off-nominal operation, the axial velocity and the flow angles will 
deviate from nominal and vary from stage to stage. However, as discussed in (Ref. 70), since the flow 
exiting a rotor or stator is essentially attached to the blade, the exit flow deviation angles are negligible 
compared to the entrance flow deviation angles (See Figure 4.7). 
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Based on this assumption, the k th stage flow angles in the cascade are 

a \,k = a 3,k-l Pl,k = Pin + A Pl,k 

®2,k = n + A^ 2 ,k and p2,k = P2n [k = \ ~ N c ] (35a) 

^3 ,k ~ ^3n P3,k ~ P3n ^P3,k 

Where A/? u and A p 3 k are the flow deviation angles during off-nominal operation. Also, since 

the initial entrance angle into the cascade is equal to the fixed exit angle of the IGV we have 

a l,k = a 3,k-l = = a IGV (35b) 

Next, in order to compute the velocities and corresponding stagnation temperatures and pressures in a 
given stage, the stage entry conditions (which are equal to the exit conditions of the previous stage) are 
assumed to be known. In this way the results from stage k can be utilized in the computation of stage k + 1 
in the cascade. Starting with stage work and mass flow expressions in (27) and (28) and utilizing (24)- 
(26) together with the ideal gas law relations, we arrive at the following polynomial equation for the k th 
stage turbine exit axial velocity C a2 ^ 


With coefficients 


C a 2 ,/ +1 + a \Ca2,k 7 + a 2 ~ 0 


ai=~(U c -2C ahk U c tana ln +2c p T 0ls , k )cos(fi 2 n) 2 


a 2 =2c cos(/? 2 „) : 


m r R 


ls,k^2s,k J 


7 - 1 


L 0ls,k 


(36a) 


(36b) 


for [k= 1 - N c ]. with 


Cal,k ^'a3,k-\'> ^0 ls,k ^03^,A:-1? P)ls,k P)3s,k-l (^7) 

Likewise utilizing the fact that U c = 0 in the stator section (or equivalently that T 03s k = T 02s k ) we obtain 
the following polynomial equation for stator exit axial velocity C a3 ^. 

C a 3 ,k +l +b\C a 3 ,k 7 +b 2 =0 (38a) 


With coefficients 


b\ - -2* c p T 02s k cos (a 3n ) 2 


b 2 = 2c „ cos(a 3 „) : 


f ■ n V" 1 

m c R 

y^02s,k^3s,k j 


rri Y 

1 02 s,k 


(38b) 


for [k = 1 - N c ] . The cascade computation is initialized at the first stage, k= 1, for a given ambient inlet 
condition T a , and P a , mass flow rate, m c , and rotation speed, Q c . The subsequent stages are computed 
using the output of the previous stage as input for the next stage. Here the k th stage rotor exit stagnation 
pressure P 0 )2sk is computed as 
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Y 

y - 1 


(39) 


A 


02 s,k 


P()ls,k 

11s, k 

f T ^ 

7 02 s,k | 

T 7 

+i 



^01 s,k ) 



with tj sk 


= \-£ 


NYk_ 

An 


Where r/ sk is the stage isentropic efficiency with scaling parameter e. Here, ij sk varies in proportion to 

the rotor entrance flow deviation angle which is directly computed from the stage input conditions based 
on (25). 


A/? u = tan 


-l 


Ur. 


C, 


- tan a Xn 


a\,k 


'An 


(40) 


Thus, at nominal conditions we have A J3 l k = 0 and i] sk =] . Finally, it is assumed that the stator stage is 
lossless, therefore, since Po3s,k ~ Po2s,k , the stator exit stagnation pressure Am,* is just 


r 03s,k 


= A 


02s, /t 


(41) 


After performing the cascade calculations (36)-(41) and proceeding through to the final stage k=N c . The 
overall compressor pressure and temperature ratios n c and r c are 


x c 


P<)2 
Po 1 


Af* and rc =?°L 

P, 7i, 


^03 s,N c 

T 

± a 


(42) 


and the overall isentropic efficiency is 


Vc 


r- 1 

n c Y 

r c ~l 


(43) 


By performing this analysis over a range of mass flows and compressor rotation speeds we numerically 
obtain the compressor pressure and efficiency characteristic maps. 

Figure 4.8 shows the pressure-ratio versus mass flow characteristic curves and beta lines for the 
7500 Hp gas turbine engine considered in this study and Figure 4.9 shows the compressor adiabatic 
efficiencies ?j c . Note r/ c =1 at the nominal operation design point. These maps are computed for the 
particular engine design and then stored digitally as lookup tables for use in the transient engine 
simulation (see Section 4.7). Finally, Figure 4.10 shows the computed stator exit axial velocities at two 
different off-nominal operating conditions (80 percent under speed and 10 percent over-speed conditions 
respectively). 
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Compressor mass flow rate, m c (kg/s) 

Figure 4.8. — Compressor pressure ratio versus 
air mass flow. 



Compressor mass flow rate, m c (kg/s) 

Figure 4.9. — Compressor adiabatic efficiency 
versus mass flow. 




Compressor Stage 

Figure 4.10. — Stator exit axial velocities at two different off- 
nominal operating conditions; (a) [Q c =0.8Q CW , m c = 0.6m cn ] ; 
(b) [Q C =1.10Q C „, m c = 1.05 m cn }„ 


NAS A/CR— 20 13-21 6502 


86 


4.6 Engine Closed-Loop Fuel Control 

To maintain the engine output shaft set-point speed during operation, a closed-loop fuel controller is 
required to regulate the fuel flow into the combustion chamber. A schematic of this arraignment is 
depicted in Figure 4.1 1. Here a Proportional Integral Derivative (PID) control law 

th f =k P / e pt + k i f J e P i dt + k d , e pt with e pt = Cl ptsei - Cl pt (44) 

is utilized for this purpose with instantaneous fuel flow into the combustion chamber m y . Where the PID 
control, with gains k pf , k If and k df , determines the fuel flow rate in order to regulate the power turbine 

set-point speed error e pt to zero. It is also necessary to account for the maximum and minimum fuel flow 
rate constraints rhf m in and my max as well as compressor maximum speed limit Q cmax which is 
accomplished through the use of two saturation blocks. Finally, to prevent windup of the PID integrator 
block, the back-calculating gains k b and k bci on the fuel saturation and compressor speed saturation 

blocks are utilized. The block diagram of the complete engine fuel control system is shown in 
Figure 4.12. 

The stability and performance characteristics of the closed-loop engine system are strongly affected 
by the gain values k pf , k df and k If . The value of these gains were selected by trial-and-error until a 

stable engine response was achieved. The engine fuel control system parameters are given in Table 4.3. 



Gas Power 

Generator Turbine Load 

Figure 4.1 1 . — Engine with closed-loop PID fuel control 
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Speed 

Limiter 


Figure 4.12. — Engine fuel control block diagram. 


TABLE 4.3.— ENGINE FUEL CONTROL SYSTEM PARAMETERS 


Parameter Value 

Proportional fuel control gain, k Pf 4.0xl0~ 3 kg 

Integral fuel control gain, k If 5.0><10' 3 kg/sec 

Derivative fuel control gain, k d 4. Oxl O' 4 kg-sec 

Compressor speed back-calculating gain, k bQ 0.018 kg/sec 

Fuel rate back-calculating gain, k b 100/sec 

Minimum fuel rate, m fmin 0.03 kg/sec 

Maximum fuel rate, rn /max 0.337 kg/sec 

Compressor speed limit, Q cmax 27,500 RPM 


4.7 Transient Gas Turbine Engine Simulation 

The dynamic simulation model of the gas turbine consists of a mixed set of ordinary differential and 
nonlinear algebraic equations. During each simulation time step, the mass flow through each component 
(compressor, gas generator turbine and power turbine) is interpolated using the numerically generated 
characteristic maps for the instantaneous rotation speed and upstream and downstream pressure states. As 
was shown in (Ref 80) a non-iterative solution to the equations can be obtained through forward 
substitution provided that the equations are solved in the proper sequence. Here equations are arranged in 
order to explicitly generate the output values by forward substitution. At each time step, the differential 
equations are integrated by means of an explicit scheme, thus obtaining the new values of the system 
plenum pressures, fuel control input and shaft rotational speeds. As was shown in (Ref 80) when the 
equations are ordered in accordance with the stream direction of the gas turbine model this results in the 
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proper solution order. Figure 4.13 to Figure 4.20 summarize the non-iterative transient gas turbine engine 
simulation model which is implemented in the Matlab Simulink environment. Finally, Figure 4.21 and 
Figure 4.22 show the gas turbine engine response to a simulated shock load event. In this loading scenario 
(Figure 4.21(a)) the main rotor drag is suddenly dropped by 50 percent at time *=5.0 sec and then 
suddenly increased *=15.0 sec. This scenario could result due to a sudden drop (or increase) in main rotor 
collective pitch. Similar shock load events can also occurs during upshift and downshift events (see 
Sections 5.0 and 6.0). This simulation shows that the closed-loop fuel control gains indeed render the 
system stable. The simulation also verifies correct functioning of the forward substitution method since 
the mass flows through each component converge to identical at steady state values. 



Rotor Drag Load 


Figure 4.13. — Gas turbine engine with fuel control and rotor load — main simulation block. 
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Figure 4.15. — Compressor simulation block. 



dP/dt 


Figure 4.16. — Inter compressor/combustion chamber plenum. 
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Figure 4.17. — Combustion chamber simulation block. 



Figure 4.18. — Gas generator turbine simulation block. 
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Time, sec 

Figure 4.21 . — Simulated gas turbine engine response to shock load; 
(a) imposed rotor drag; (b) compressor speed; (c) power turbine. 



Figure 4.22. — Simulated gas turbine engine response to shock load; 
(a) stagnation temperatures; (b) fuel flow rate. 
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5.0 Gearbox Dynamics Modeling 

5.1 Introduction 

It is well known that gear tooth meshing stiffness fluctuations and tooth backlash effects are a 
significant source of vibration and instability in geared rotor systems. The overall objectives of this 
portion of the project can be summarized as follows; 

• Develop gear system dynamics models to predict torsional vibration levels during variable speed 
operation 

• Conduct experimental evaluation of gear dynamics model under variable speed conditions 

• Develop computationally efficient tools to asses gear vibration levels during variable speed 
operation 

• Explore gearbox/shaft dynamic interactions and stability behavior via variable speed dynamic 
simulations 

The topic of gear induced vibration and its mitigation has been widely studied and remains an 
important topic in propulsion system research and design. For meshing gear pairs, the number of teeth in 
contact is kinematically a function of the gear rotation angle. Since gear mesh stiffness due to elastic tooth 
bending is related to the instantaneous number of teeth in contact, this results in time -varying stiffness 
fluctuations which can excite vibrations and induce parametric instability (Refs. 84 to 95). For a spur gear 
pair, References 84 to 87 experimentally demonstrated that gear mesh stiffness parametric effects excited 
large amplitude period-1 (17) and period-n (nT) subharmonic responses. For gear mesh harmonics in the 
vicinity of natural frequency multiples, these vibrations resulted tooth loss-of-contact which produced 
nonlinear chattering and jump phenomena (Refs. 85 and 86). Several reduced-order lumped-parameter 
spur gear dynamic models including time-varying mesh stiffness and clearance nonlinearities were 
analyzed by References 85, 86 and 88. Furthermore, References 89 and 90 explored the effects of mesh 
stiffness variation on two-stage gear-trains and planetary gear transmissions. Here the undamped 
equations-of-motion with an idealized square wave type mesh stiffness was considered and conditions for 
primary, secondary and combination parametric instabilities were derived via perturbation analysis. 
Additionally, the roles of mesh damping, involute contact ratio, and tooth sliding friction on spur gear 
dynamics have been explored (Refs. 91 to 95). 

Another body of gear vibration research has focused full scale system modeling of coupled 
gear/rotor/bearing vibrations (Refs. 96 to 101). However, in most of these studies gear backlash and mesh 
stiffness variations were not considered since the scope was primarily on modal analysis and coupling 
between transverse and torsional motions (Refs. 97 and 98), and torsion, lateral and axial motions 
(Refs. 99, 100, and 12) in large degree-of-freedom gear/rotor/bearing systems. These studies typically 
considered a constant (nominal) value of gear tooth mesh stiffness together with shaft and bearing 
compliance. The results offered valuable insights into the effects of gear, shaft and bearing placement on 
modal coupling and natural frequencies. However, the influence of the coupled substructure interactions 
on gear-mesh induced parametric instabilities was not determined. Finally, Reference 101 considered a 
periodically stepped alloy shaft for reducing torsional gear vibrations via wave propagation stop-bands. 
Here however, the focus was on blocking transmitted vibrations and not on suppressing the gear tooth 
mesh instabilities. 

Despite this large body of gear dynamics research, the majority of the work has focused on constant 
speed operation where the system is linear periodically time-varying. The exceptions to this are 
investigations into nonlinear backlash effects (Refs. 85 to 88). However, these studies focused on quasi- 
steady response under slow-speed RPM changes and did not characterize transient gear vibrations nor the 
effects of spinup or spindown acceleration rates. The following sections will address these issues as well 
as the four main objectives stated at the beginning of this section. 
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Figure 5.1 . — Spur gear/shaft system dynamic model. 


5.2 Spur Gear/Shaft Structural Dynamics Model 

To analyze the gear system dynamics, a structural dynamic model of the spur gear/shaft system 
depicted in Figure 5.1 is developed. Here, a torsionally flexible input shaft, length^, drives a load inertia, 
J L , and a resistive load torque, 7), through as single spur gear pair. In this model, the shaft elastic twist 

(p s (x, t) is defined over the domain [0 < x < L s ] . 

On the input side (at x = 0). The shaft is rigidly fixed to the input coupling thus has boundary 
condition 


^(o,0 = o (i) 

Furthermore during operation, it is assumed that the input shaft rotation, $>, is driven with angular 
velocity Q 


<f>o = ( 2 ) 

Based on this, the total rotations and ^ g2 (0 of gears 1 and 2 are 

0 g i(t) = Qt + <4(4,0 and </> g2 (t) = -n n Q.t + 0 g2 (t) (3) 

where gear 1 is fixed to the shaft at x = L s and is the elastic twist deviation (due to tooth mesh 

compliance) from the nominal rigid-body rotation dictated by the gear ratio n u 

«12=| L (4) 

where Z x and Z 2 are number of teeth on gear 1 and 2 respectively. 

5.2.1 Lumped Spur Gear Model 

Here the gears are modeled as lumped rigid disks with base radius, R hh mass, m gi , and rotational 
inertia, J gh [for i= 1,2]. The tooth mesh compliance is accounted for by a spring stiffness element, k m (t ), 
acting along the line-of-action (at pressure angle a PA ). See Figure 5.2. 
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Figure 5.2. — Lumped gear-pair model, inertias and tooth mesh stiffness. 

For gears with non-integer contact ratio, cr, the number of teeth in mesh fluctuates between the 
integers ceil(cr) and floor(cr) as a function of gear rotation. 

cr - J (Rp\ + a ) 2 - Rb\ + a / (.Rpl + a ) 2 - Rb2 (Rpl + Rp2)^ na PA 

cp 

Where cp is gear circular pitch, R pi [7=1, 2] is pitch radius 

Z- cp 

R p . = -4 — and R bi = R pi cos a PA [ i = 1,2] 

2 K 


(5) 


( 6 ) 


and tooth addendum 


n 

Since k m (t) is proportional to the number of teeth in contact at any instant, at constant rotation speed, k m (t ) 
switches between maximum and minimum stiffness values £ max and k min with square wave-like variation 
which has Fourier series representation 


N h 

k m (t ) = k 0 + k sn sm(nZ x £lt) + k cn cos (nZ^t)] (8-a) 

n - 1 


with 


*o = (2 - cr)k mm + (cr - 1 )k max 


(8-b) 


and 


7 n 1 - cos \2n{cr — 2)n\ 7 _ n 7 sin \2n{cr - 2)n\ 

ksn (^max ^min) 3Iid k cn (^ max ^ m in ) (^"^) 

no- 
where Nh is the number of expansion harmonics and Z x £l is the gear mesh frequency (and associated 
fundamental period T p = 2tt/ZiQ). Similar to Reference 102, the gear mesh stiffness is estimated by 
considering each tooth as a tip loaded cantilever beam with average tooth dimensions in series with gear 
body stiffness values k hx and k b2 . For [1 < cr < 2] 


NAS A/CR— 20 13-21 6502 


97 



1 


1 


]r . — 

^min j 


2 

H h 

kbl K k b 2 


— and k mm = — 


1 1 

1 1 

kb\ K k b 2 


and tooth stiffness, k h and gear body stiffness k bi 

E cp 3 d 


K = 


32cos a PA (a + b) 3 


and k bi = 


4 G7rR hi d 

Rhi — Rhi 


[/=1,2] 


( 9 ) 


( 10 ) 


with gear elastic and shear moduli, E and G, tooth dedendum, 6=1.25*a, face width d , and hub radius R hi 
[i= 1, 2]. 

5.2.2 Work and Energy Expressions 

Based on the above lumped gear model approximation, the gear/shaft system kinetic energy, 7, can be 
expressed as 


T = \psJs\ L q ' [kx,t) + nfdx + \vJ 2 gX (t) + (. J g 2 + 2(0] ( 11 a) 

where “* indicates differentiation with respect to time, t , Also, p s is the shaft material density and J s is 
the shaft polar moment of area defined in terms of shaft inner and outer radii, R si and R so , 

J s =^(Rj-R si 4 ) (lib) 

Furthermore, the total system strain energy, J7, due to both shaft torsional elasticity and gear tooth mesh 
compliance is 


1 5 1 

u = ~\ G s J s fc(x,t) 2 dx + —k m (t)[R b i</> g i(t) + R b2 <f> g 2 (t)] 2 (12) 

2 o 

where “ ? indicates differentiation with respect to the shaft axial coordinate x, and G s is the shaft material 
shear modulus. Note, in the second term of the above strain energy expression, the quantity 


u m (t) = R b \</> g \(t) + Rbi&giit) 


(13) 


is the gear tooth mesh deflection u m (t ) along the line-of-action. In addition, dissipation effects are 
accounted for via the Rayleigh dissipation function 

D = U S \ Ls G s jj' s {x,tf dx + U m k m {t)[Rj gl (t) + R b2 i g2 (t)? (14) 

2 z 

where ^ and are shaft and gear mesh proportional viscous damping parameters. Finally, the virtual 
work due to the driven load torque T L and prescribed angular acceleration rate Q is 


SW -~(T l +n n (J L + J 


(15) 
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Figure 5.3. — Two-DOF shaft torsional finite element. 

5.2.3 Gear/Shaft Finite Element Model 

In order analyze the coupled gear/shaft system structural dynamics, a finite element model (FEM) of 
the gear/shaft system in Figure 5.1 is developed. The shaft is discretized into n el , uniform torsional 
elements each with length L d =L s ln e i. Each element, as shown in Figure 5.3, has two nodes with 1 torsional 
degree of freedom (DOF) per node. The j th element DOF vector is 


qelj(t) 


' fj{t) 
^■+i(0 


(16) 


where (/>/(! ) is the shaft elastic twist DOF at the j th nodal location. 

Using the standard linear interpolation element, the shaft twist, <j> s j(x,t ) , within j th element domain 
[0<x<4,] is 


4,7 CA 0 N e i (0 

with corresponding elemental shape function vector 

N d (x) = \l-x/L el x/L el ] (17b) 

After fusing the adjacent shaft nodes and applying the constraints (1) and (3), the global DOF vector of 
the discretized gear/shaft system becomes 


q(t) = 


Qsif) 

j > g 2 (0 


Which is composed of shaft twist DOF vector with N n nodes 

= k,i (0 ••• A,n„ ( t)Y 

and the output gear elastic rotation (/) g2 (t ) . 


(18a) 


(18b) 
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Next, utilizing the system work, energy and dissipation expressions in (1 1), (12), (14) and (15) 
together with the FEM discretization (17), the overall gear/shaft system equations-of-motion are 
synthesized via Lagrange’s Equations 


d_ ST 
dt _ dq 


dT dU 8D 

1 1 

dq dq dq 


F 


The resulting linear time-varying equations-of-motion are 

M<?(0 + C(t)q(t) + K(t)q(t) = F (t) 


With system inertia matrix 


M = 


M s + J gi N^N gi 

0 


Jg2 + JL 


damping matrix 


c(0 = 


Zs K s + %mk m (t) R b\ N ll N g\ Zmk m (t)R bl R bl N T gl 
‘vmkm (t )R b \ R b2^ g] ‘ymkm (f)R b 2 


stiffness matrix 


K(0 = 


K s + k m mZX gl N gl k m (t)R h[ R h2 N T gl 
k m (t)R bl R b2 N gl k m (t)R b2 


and generalized force vector F obtained from 5W = F T Sq 

-J gl QN/ 

Jl +n u(^L + Jg2)& 


( 19 ) 


(20a) 


(20b) 


(20c) 


(20d) 


(20e) 


Here, the DOF selection matrix, N gl , defined in (21), extracts the gear 1 twist DOF (i.e., shaft twist at 
x = L s ) from the shaft twist DOF vector q s 

i(L s ,t) = [ 0 o ••• 0 l]q s (t) = N gi q s (t) (21) 


Also, the mass, stiffness and damping submatrices, M s , K s and C s in (20) are formed according to the 
FEM approach from the elemental mass stiffness and damping matrices given as 


and 


ei 

M s>el = PsJ s J Nli ( x )N el (x)dx and K s el = G S J S ^ ' N'J (x) N' el (x) dx 
o 


(22a) 


C s ,ei = ZsG s J s \l S N'J(x)N' el (x) dx 


(22b) 
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Figure 5.4. — Gear backlash; (a) gear pair with tooth clearance (b) gap function. 


5.2.4 Gear-Tooth Clearance and Backlash Effects 

Next the model is modified to include nonlinear backlash effects due to tooth gap clearance h gap . 
Specifically similar to References 85 to 88, the so called gap function is introduced (see Figure 5.4). 
Based on this, the gear/shaft system equations-of-motion including nonlinear backlash become 




<7s( 0 

g[4 2(oi 


= F 


with gap function 


gWgiit)]- 


gap 


2R 


+ ■ 


b 2 


* g a P 


2 R 


-b 2 
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gap 


2 R 


I $g2 I — 


b 2 

h 

n gap 

2 R 


< — 


b 2 

h 

'gap 


2 R> 


b 2 


(23a) 


(23b) 


Where the system mass, stiffness and damping matrices and generalized force vector are identical to those 
defined in (20). 


5.3 Single Gearbox/Shaft System Variable RPM Response 

To explore the variable speed response of gearbox/shaft system presented in the previous section, the 
parameters for a representative rotorcraft driveshaft subsystem are utilized for analysis. These parameters 
are summarized in Table 5.1. 

Figure 5.5 shows the [Q = 0 to 3000 to 0] RPM spinup-spindown vibration response of the gear mesh 
deflection (13) for three load torque levels. 

Here the variable tooth stiffness results in vibratory behavior under constant load torque conditions 
due to parametric excitations. The dotted lines at ±10//m in Figure 5.5 indicate the backlash limits 
[-hgap/ 2 <u< h gap /2\ also shown Figure 5.4(b). Tooth deflection greater than h gap l 2 indicates front 

side tooth contact and tooth deflection less than - h gap / 2 means back-side contact. Deflections within 
[-hgap /2<u< h gap / 2] indicate loss-of-contact. Figure 5.5(a) shows that for lower torque values the tooth 

remains in single-sided contact (deflection remains near upper limit). As the load torque is increased, 
Figure 5.5(c) shows that the tooth begins to chatter back and forth between the upper and lower limits 
(double sided contact). This double sided contact occurs for certain rotational speeds only. The magnitude 
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of the chattering is proportional to the load torque level. It is also noted that the spinup and spindown 
vibration amplitudes are not the same. This difference is again a result of the nonlinear dynamics nature 
of the system. 

Utilizing the above nonlinear dynamic model for the gear/shaft system, the variable speed dynamic 
behavior is further investigated in Section 5.4. To assess the validity and performance of this model, the 
analytical transient vibration response predictions will be compared with measurements from a variables 
speed gearbox/shaft test rig experiment. 


TABLE 5. L— GEARBOX-SHAFT SYSTEM PARAMETERS 


Parameter 

Value 

Number of teeth, Z { 

Z 1 =42, Z 2 =30 

Pressure angle, a PA 

25° 

Circular pitch, cp 

0.0157 m" 1 

Contact ratio, cr 

1.4954 

Face width, d 

18 mm 

Gear Inertia, J t 

..J! =0.243, J 2 =0.063 kg-m 2 

Load Inertia, J L 

J L = 2.0 kg-m 2 

Load Torque, T L 

Various* 

Backlash gap clearance, h 

20 jam 

Shaft Length, L 

6 m 

Density (shaft) 

2800 kg/m 3 

Elastic modulus (gear), E 

200 GPa 

Shear modulus (shaft), G 

27 GPa 

Cross-section polar area moment, J s . 

4.469xl0" 6 m 4 


* Various values used in the investigation 


Tooth 

Deflection, 

pm 



Spinup 


Time, sec 
>< 


Spindown 


(+100 RPM/Sec) 
[Q = 0 to 3000 RPM] 


(-100 RPM/Sec) 
[Q = 3000 to 0 RPM] 


Figure 5.5. — Gear mesh vibration response for 0-3000-0 RPM ramp input; 
(a) T l = 100 Nm; (b) 71 = 200 Nm; (c) 71 = 300 Nm. 
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5.4 Variable Speed Gear/Shaft Model Experimental Validation 

To account for effects of gear meshing on the system behavior, a nonlinear time varying dynamic 
model including the effects of gear tooth flexibility, non-unity contact ratios, backlash clearance, gear 
inertias and load torque has been developed (see Section 5.2). To validate this model, an experimental 
gear dynamics test rig at the University of Tennessee has been utilized, see Figure 5.6. 

The test rig consists of a DC electric motor with PID speed controller which drives a single-stage spur 
gearbox (also capable of multistage). The input shaft is fitted with an inductively coupled torque and 
RPM sensor and an optical encoder. The gearbox output drives a magnetic brake which has a manual 
brake torque setting. Table 5.2 and Figure 5.7 summarize the experimental test setup. 



Brake 


Stage 
Gearbox 
(Spur Gears) 


Torque 

Sensor 


DC Motor 


Figure 5.6. — Spur gear dynamics test rig (Univ. Tenn.). 




Z.,=30 

(input) 


Z 2 =42 

(output) 


Figure 5.7. — Spur gear dynamics test rig experimental setup 


TABLE 5.2.— GEARBOX DYNAMICS TESTRIG PARAMETERS 


Gears 

Input shaft 

Gear type Spur 

Material Steel 

Length 18 in. 

Diameter (solid) 5/8 in. 

Number of teeth (input) Zi=30 

Number of teeth (output) Z 2 =42 

Tooth module 5 mm 

Diametral pitch 200 

Pressure angle 25° 

Addendum 5 mm 

Dedendum 6.25 mm 

Face width 25 mm 

Gap clearance (nominal) 50 pm 

Motor speed controller 

P-gain ,k p 1^10^ 

I-gain, hi 1><10”6 

D-gain, k d 20 

Load 

Rotational inertia, J L 0.0632 kg m 2 

Brake settings, T L 10, 20 and 30 in-lb 
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Finally, the gear pair contact ratio of this spur gear pair is calculated to be c r = 1.495. 

To perform the experimental validation, several variable RPM tests under three different brake 
settings (load torque, T L = 10, 20 and 30 in-lb) were performed to validate the high frequency gear mesh 
dynamics model. In each test the motor set-point RPM was ramped linearly from [0 - 3000 - 0] RPM over 
60 sec while the input shaft speed and torque data were collected. Figure 5.8 shows the ramped RPM 
profile and Figure 5.9 to Figure 5.11 show the resulting input shaft load torque. Here Figure 5.9 to Figure 
5.1 1(a) give the experimentally measured load torque while Figure 5.9 to Figure 5.1 1(b) show the 
corresponding predicted values from the model. 





time, sec 


Figure 5.9. — Input shaft torque versus time during 0-3000-0 RPM 
ramp test 1 , J L = 10 in-lb, (a) UT gearbox test rig (Exp), (b) gearbox 
dynamics model (SIM). 
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Figure 5.10. — Input shaft torque versus time during 0-3000-0 RPM 
ramp test 2; 71= 20 in-lb, (a) UT gearbox test rig (Exp), (b) gearbox 
dynamics model (SIM). 




time, sec 


Figure 5.1 1 . — Input shaft torque versus time during 0-3000-0 RPM 
ramp test 3; 71= 30 in-lb, (a) UT gearbox test rig (Exp), 

(b) gearbox dynamics model (SIM). 


In each test, the torques predicted by the model agree with the measured dynamic torque behavior for 
the variable RPM condition. During both spinup and spindown, the model predicts jumps in dynamic 
torque amplitude. The main discrepancies occur just after the transition from spin-up to spin-down (near 
t = 30 sec). Here the measured transient torque spikes are larger than the model predictions. This could be 
attributed to some nonmodeled system looseness or backlash in the gear-shaft spline connections. 
However, as seen in Figure 5.9 to Figure 5.1 1, the dynamic torque amplitude predicted by the model are 
comfortably within 25 percent of the measured values during the constant spinup (0 < t < 30 sec) and 
spindown portions (30 < t < 60 sec). Furthermore, the jumps occur nearly at the same RPM values. This 
jump phenomena is due to both linear and nonlinear effects. As the mesh harmonics nZ x £l x pass through 
the system natural frequencies resonances are excited however nonlinear behavior due to gear backlash 
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and tooth clearance are also a critical feature of the model and necessary to accurately predict the 
measured values. 

5.5 Dual Gearbox/Shaft System Interactions Under Variable RPM 

In this section the nonlinear dynamic interactions between two gearboxes coupled by a long 
torsionally flexible shaft. Such configurations are common in rotorcraft propulsion systems including 
tailrotor and cross-shafting drivelines. The shaft flexibility induces dynamic phasing between the two 
upstream and downstream gearboxes. This dynamic phasing together with the kinematic dependence of 
the gear tooth mesh stiffness on the gear rotation induces a nonlinear dynamic interaction not explored in 
previous single-stage or rigid-shaft gear system investigations. 

For constant speed operation, Reference 90 showed that gear mesh stiffness variations in a single 
stage gear system can be effectively modeled as a square wave with frequency N t Q where N t is the 
number of teeth in mesh and Q is the rotation speed. In the case of multistage gear systems, the critical 
importance of relative gear mesh phasing was demonstrated in Reference 13. Here, the concept of static 
mesh phasing, where individual gear mesh pairs are phased to cancel out or suppress the vibration from 
previous stages, was explored in the design of a multistage planetary gear system. For constant speed 
operation, the “built in” or static design mesh phasing was effective at reducing mesh stiffness induced 
gear system torsion vibrations. 

It is clear from the existing research that the relative mesh phasing between different gear sets in a 
multistage gear train is a key design parameter which has a significant impact on the system vibration and 
stability characteristics for constant speed operation. However, the case of variable speed operation has 
not been explored. Furthermore, previous investigations have only considered the case of constant mesh 
phase between two gear sets, however in certain situations (such as when long cross-shafting and large 
dynamic windup is involved) this assumption may not be adequate. This paper explores nonlinear time- 
varying torsional dynamic behavior of a gear driven system involving two gearboxes coupled by a 
torsionally elastic shaft. See Figure 5.12. 

Such configurations are common in many rotorcraft propulsion systems including tailrotor and cross- 
shafting drivelines. In this situation, if the flexible shaft connecting the two gear-sets has a nonnegligible 
windup, the relative phase becomes a dynamic variable related to the shaft twist resulting in a nonlinear 
time-varying gear mesh stiffness. Nonlinear dynamic phase interactions of this nature between gear sets 
has not been investigated in the literature, hence the aim of this study is to explore this effect. 
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5.5.1 Dual Gearbox/Shaft System Model 

To facilitate analysis, the equations-of-motion of a dual stage, spur gear/shaft system including gear 
mesh stiffness variations and shaft torsional dynamics is derived. See Figure 5.12. Here, gears in mesh 
pair a-b are modeled as rigid disks, with rotations c/) a and rotational inertia J a and J b , and number of 
teeth Z a and Z b . Gear tooth meshing is accounted for by spring stiffness, k mab , acting along the line-of- 
action in series with backlash gap clearance h. 

Due to the periodic nature of the mesh stiffness with respect to gear rotation, the gear mesh stiffness 
functions, k mn (</\) , h m34 (^) for the upstream and downstream gearboxes / and II (see Figure 5.12) can 
be expressed in a Fourier series as 


k mn (<h) =k p cr + 2k p Y l Mn7Z — cos(wZi^) 
n = 1 nn 

k m , 4 (<fa) = k p cr + 2 k p ^ sm(n7r cr) cos (nZ 3 <fa + y i n ) 


(24) 


n = 1 


T171 


Where y I>n is the relative mesh phase angle between gearboxes / and II (which is a free design parameter) 
and N h is the number of expansion harmonics. By considering the intermediate shaft elastic twist, <j) s (x,t ) , 
the elastic gear mesh deflections, U\ 2 {\) and u 34 (t), and gear ratios n u and n 34 


R b \ Z\ R b 3 Z 3 

n n = = ~ and « 34 = 

Z 2 


Rfi 2 Z 2 Rb 4 Z 4 

the total gear rotations, fait) [i= 1-4], of the dual gearbox-shaft system in Figure 5.12 are 

1 


(25) 


fa(t) = Clt, (/> 2 (t) = ——u n (t)-n 12 nt 

R b2 


<h (0 = -T- “12 (0 - n n&t + (4 (L, t ) , <j) A (0 = l -b22l2 + „ 34 
^02 


“12(0 

Rb 2 


+ n n Q.t-<j) s (L,t ) 


(26a) 

(26b) 


with shaft elastic constraint 


(4(0,0 = 0 (26c) 

Here Q is the driven input speed, L is the intermediate shaft length, and </> s (L,t) is the total shaft elastic 

windup between the upstream and downstream gearboxes / and II. Note, the mesh deflections can be 
expressed in terms of gear rotation and base radii as 

u n(t) = Rbl<A(0 + Rblfait) an d M 34(0 = ^63^(0 + ^64^4 (0 (27) 

where zero mesh deflection corresponds to perfectly rigid gear rotations dictated by ratios n n and n 34 . The 
total system kinetic energy T is expressed as 

2 o 2 
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where p is the shaft material density, J s is the shaft cross-sectional polar moment-of-area (units m 4 ), and 
J L (units kg-m 2 ) is the load inertia. Furthermore, the system strain energy, V 9 due to shaft flexibly and 
mesh deformations is 


i L 


V = ^\ GjJ' s (x,t) 2 dx + ^k mn (<A)[gn(Un)] 2 + ^k m34 (<k)[g 34 (u 34 )] 2 


(29) 


Where “ ' “ indicates differentiation with respect to the shaft axial coordinate x, G is the shear modulus 
and gn(ui 2 ) are £34(1/34) are functions with account for the backlash gap clearance, h, as 


^12(^12) - 


u u (t)-h/2 u n (t) > hi 2 
0 \u u (t)\<h/2 and g 34 (u 34 ) = 

U\i (t) + h/ 2 u n (0 <-hl2 


u 34 {t) -hi 2 u 34 (t ) > hi 2 

0 | u 34 (t) |< hi 2 

u 34 (t ) + hi 2 u 34 {t) < —h!2 


Also, the virtual work due to a driven load torque T L is 

Su 34 (t) 


SW - -T L 8(j) 4 = —T l 


R 


+ ?L n 34 


b 4 


Su n (t) 


R, 


+ 8(fi s {Lp) 


b 2 


(30) 


(31) 


To obtain equations-of-motion, the shaft twist angle deformations are first expanded using an 1 *N m vector 
of assumed modes 0(x) as 


0 s (x,t) = Q>(x)7] s (t) with O(x) = [sin j5 \ y sin ... sin^y] (32a) 

i-/ -Z-/ X-/ 

Where r/Jf) is the N m x 1 shaft modal coordinate vector and the modal parameters [i= 1 - N m ] are 

determined via numerical solution of 


cot p n = 


./ 3 + (,/ 4 + J i ) 

pLJ s 


Pn 


(32b) 


Which is the exact characteristic equation corresponding to the case k mn > co and k,„ > go . 

Using the work and energy expressions in (28), (29) and (31) together with the kinematic relations 
(26), and the modal expansion (32), the equations-of-motion are synthesized via Lagrange’s Equations as 

d_\dT_ 

dt dq 


8T 8V „ 

+ — = 2 

dq dq 


(33) 


with the (N m + 2) x 1 generalized coordinate vector q(t) = \>] s (t) T U\ 2 (<) u 3A {t)J which is comprised of 

the shaft modal coordinates and the gear mesh deflections. 

The resulting nonlinear time varying equations-of-motion of the dual gearbox-shaft system written in 
terms of q and the and £12(1*12) are £34(1/34) are 


ds 


Mq + Cq + K(t,q)\ 


gu(u l2 ) 


£34(1/34) 


F(t,q) 


(34a) 
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With mass matrix 


M = 


pJs [ >dx + ( J 3 + K34 2 i-h + Jl ))® l T ®l [ ® T dx + <F 

0 R b2 j 

M 

Rb2 


^2 ^4 


- J O dx 


+ 0 


♦^ 3+^34 id 4 + Jl ) 


pJ s L + J 2 +J 3 +n 34 (J 4 +J L ) 


n 34 


b2 

K b4 
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n 34 


Rh?Ri 


-Va+Jl) 


^b2^bA 

1 


Va+Jl) 


b2*b4 


Ri 


(J 4 + J L ) 


b 4 


(34b) 


with nonlinear time varying stiffness matrix K 


K(t,q) 



0 


k mi2 (0 


dgn 

dtl\2 


0 

0 


0 


0 


k m}4 (<h(t,u n ,ij s )) 


dgM 

du 3 4 


and with nonlinear time varying force vector F 


F{t,q) 


Jl_ 

R b2 


n 34 R b2 4> T {L) 
n 34 

-1 


1 ^m 34 (^3(^12,^)) 

2 dtjh, 


[^ 34 ( m 34)] 2 


4> r (I) 

1 

R b 2 

0 


Furthermore, the gap function derivatives become the so-called separation functions 

dgn _ [ 1 \un(t)\>h!2 ^ Jg 34 [l \u 34 (t)\> h/2 

du l2 [0 | w 12 (0 1< /2 du 24 [O \u 24 (t)\<h/2 


and the mesh stiffness derivative in terms of the Fourier expansion in (4) is 


dkm 34 (&) 

d(j > 3 


. 2t z 3 y!2a££2 

71=1 ^ 


sin(nZ 3 ^,(t,M l2 ,/ 7 s ) + y /// ) 


(34c) 


(34d) 


(35) 


(36) 


In this system, nonlinearly arises from the tooth backlash effects and from shaft elastic windup. The latter 
effect is due to the dependence of the downstream gearbox input rotation = </h,(t,u u ,ri s ) on the shaft 

modal coordinates, r/ s (t) and the stage / mesh deflection Un(t). Hence the shaft windup essentially acts as 
a dynamic mesh phase term. The nonlinear coupling arises due to the dependence of the mesh stiffness 
k m 34 ($3) on ^3 . Previous nonlinear gear dynamic investigations have only considered the backlash as a 
source of nonlinearly which is valid for a single-stage gearbox or with a rigid shaft. 

5.5.2 Nonlinear Harmonic Balance and Continuation Analysis 

To account for the effect of dynamic mesh phasing due to shaft flexibility and backlash the full 
nonlinear time -varying system is analyzed using a harmonic balance and numerical continuation 
approach, (Refs. 92 and 93). 
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To proceed with the analysis, it is recognized that the stiffness variations are periodic with respect to 
the input rotation angle. Thus harmonic and (subharmonic) steady-state solutions of the form 


q(t) = q 0 + 


N 


n= 1 L 


nZ x . nZ x ^ 

q cn cos Q t + q sn sin Q t 


(37) 


Where the vectors q 0 , q cn are the q sn unknown harmonic solution amplitudes and the positive integer s is 
the subharmonic index which is assumed (Refs. 87 and 88). The assumed value of s will generate sT p 
periodic solutions where T p = 2x1 Z X Q is the fundaments gear mesh period. Furthermore, based on (34), 
the backlash functions can be expanded as 


£12(^12) - §12,0 + ^ 

n - 1 

N 

g3 4 («3 4 ) = g 34 0 + ^ 


n - 1 L 


nZ, . «Z, 

gu,cn cos Qf + g u , sn Sin fit 

s s 

nZ x . nZi _ 

§ 34 , cn cos Qt + g 34;i „ sin Of 

s s 


(38) 


where the backlash function Fourier coefficients are expanded via a Discrete Fourier Transform (DFT) by 
sampling (37) at NDFT points and subsiding the result into (38) we have 


0 k =^P^,k = [O,l ... ,N dft -\] 
JS! dft 


1 Ndft 1 Ndft 

§12,0 = — T.gl2(um) §34,0 = — y,g34(»34./c ) 

M DFT k=0 M DFT k =0 

2 Ndft 2 Ndft 

§\2,cn ~~ ^g'l2,i( M 12,*)C0S^ : gs4,cn =~ T.g 34,k( u 34,k) cos &k 


(39a) 


DFT k =0 


N, 


(39b) 


DFT k=0 


2 y v DFT 2 iy DFT 

§12, sn = — ^ gl2,k ( u l2,k) s i n @k §34, sn = ~TZ ^ g34,k ( u 34,k ) s i n 

^ DFT k=0 M DFT k=0 

Upon substituting expansions (37) to (39) into the equations of motion, (34), and equating like harmonic 
terms, the a nonlinear algebraic system is obtained, 

S(X, Q) = 0 (40a) 

where Xis the 9U e (2N h +\) NDOF x 1 solution vector of steady-state Fourier expansion vectors. 


j 1 j 1 j 1 T 1 t 

q sl q cX - q sn q cn ] 


(40b) 


Next, the solution of the nonlinear algebraic system in (40) for various speeds, Q, is obtained via the 
numerical arc-length continuation approach (Ref. 103), outlined below. 

Essentially, the arc-length continuation approach is a prediction-correction strategy for tracing the n- 
dimensional solution curves in the coordinate space with respect to a continuation parameter. In this case, 
the solution curve y e 9U+1 is composed of the system degrees of freedom, X , and the rotation speed Q. 

y = [X r ,Q] (41) 
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By starting at a known solution point, y k , a prediction of the neighboring solution, y k+x , is made along the 
local unit tangent vector, y h of the solution curve. 


yk+ 1 -yk + h v k (42) 

Where h k is a chosen step-size parameter and the unit tangent vector, y k , is be computed from the system 
Jacobian J(X, Q) = and the previous tangent vector v k - 1 as 



-1 

o“ 

T 

L v k - 1 


l 


(43) 


Finally, the prediction, y k+1 , is refined via Newton-Raphson iteration indicated by the subscript i. 


-U 


yk+l,i+\ yk+\,i J aug ^ aug 

Where the augmented system and Jacobian, S aug and J aug , 


u aug 


S(y k ) 

and J m .„ = 


.( yM-yifv k _ 

aug 

. V k 


(44) 


(45) 


ensure that the Newton corrections are constrained along a path normal to v k . This nonlinear solution 
procedure is shown schematically in Figure 5.13. 

Once the nonlinear system steady-state n/rev solution amplitudes, X, are calculated, the stability of the 
solutions can be determined via Floquet Theory the by assuming perturbations about X. 


prediction 
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5.5.3 Dual Gearbox/Cross Shaft Response Results and Observations 

Utilizing the harmonic balance and continuation analysis methodology developed in the previous 
section, the sT p period steady-state solutions are computed for a constant load torque level 7^=100 N-m. 
Table 5.3 summarizes the dual gearbox/shaft system parameters utilized in this numerical investigation. 

Figure 5.14 shows the resulting \ T P ,2T P and 2T P period solutions versus operating speed for to values 
of gap clearance. 

Here the effect of the tooth gap clearance is clear illustrated. For the case with no gap (i.e., the h gap = 0 
jum case), the response harmonics experience ordinary resonance-like behavior. However, once the 
backlash is included, (i.e., the h gap = 20 jum case), the response is significantly different. The presence of 
the gap clearance produces a softening “backbone” behavior of the resonance curve. This is a 
characteristic of a system which undergoes jump phenomena during spinup or spindown resulting in 
different limit cycle amplitudes depending on direction of the RPM change. Despite this behavior, 
backlash the nonlinearity has the beneficial effect of significantly reducing the vibration amplitudes 
compared with the h gap = 0 jum case. 

TABLE 5.3.— DUAL GEARBOX-SHAFT SYSTEM PARAMETERS 


Parameter Value 

Number of teeth, Z 1 =42, Z 2 =30, Z 3 =42, Z 4 =30 

Pressure angle, a PA 25° 

Circular pitch, cp 0.0157 m' 1 

Contact ratio, cr 1 .4954 

Face width, d 18 mm 

Gear Inertia, J t J^O.243, J 2 =0.063, J 3 =0.234, J 4 =0.063 kg-m 2 

Load Inertia, J L J z =2. 0 kg-m 2 

Load Torque, T L 500 Nm 

Backlash gap clearance, h Various * 

Shaft Length, L 10 m 

Shaft torsional stiffness, k t Various * 

Density (shaft) 2800 kg/m 3 

Elastic modulus (gear), E 200 GPa 

Shear modulus (shaft), G 27 GPa 


* Various values used in the investigation 



Q, RPM 

Figure 5.14. — Dual gearbox/shaft vibration response harmonics versus Q with kt = 6000 Nm, 
for two values of gap clearance; (a) h gap = 0 //m; (b) h gap = 20 jum. 
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u 34 \, jum 
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Q, RPM 

Figure 5.15. — Dual gearbox/shaft vibration response harmonics versus Q with h gap = 20 jum 
for; (a) k t = 6000 Nm; (b) k t = 2000 Nm; (c) k t = 1000 Nm. 


Next, Figure 5.15 shows the effect of the dynamic mesh phase by comparing the dual gearbox/shaft 
system vibration response for three different shaft torsional stiffness values [k t = 6000 Nm, k t = 2000 Nm, 
and k t = 1000 Nm]. 

Here the effect of dynamic mesh phase is seen by the increased softening of the response backbone 
curves at each harmonic. To explore the onset of chattering behavior, the steady-state tooth impact 
velocities are plotted versus the shaft rotation speed in Figure 5.16. Tooth impact velocity v c is defined as 
the tooth velocity at the instant when either front-side or back-side tooth contact occurs. This is evaluated 
from time-domain simulation as 


V c = u with \\u(t c )\ - h \ < s (46) 

14 l c I 1 1 I 

where t c is the time of impact determined from the £ criteria in (46) where sis small positive number. In 
this simulation the both gears pairs are operating in full contact for [1000 < Q <1247] RPM. Above this 
threshold operating speed, chattering begins. In the range [1247 < Q < 1328] RPM, the chattering 
response executes alT p orbit. However over [1328 < Q < 1823], the behavior appears chaotic. For [1823 
< Q < 2305] RPM, 2T P behavior emerges with and IT P orbit reappearing for Q > 2305 RPM. 

Finally, Figure 5.17 shows the gear mesh deflection (13) vibration response of the dual gearbox/shaft 
systems during a [Q = 3000 to 0 RPM] spindown for two shaft torsional stiffness values \k t = 15000 Nm 
and k t = 6000 Nm]. 

Here the increased shaft torsional flexibility (and increased dynamic mesh phase) has both stabilizing 
and destabilizing effects depending on the RPM range. Specifically, by comparing Figure 5.17(a) and (b), 
over the shaft speed range [Q « 1250 - 1700 RPM], the reduced shaft stiffness case (Figure 5.17(b)) tends 
to destabilize the higher amplitude 27^, limit cycle causing the response to jump down to the lower 
response amplitude \T p , limit cycle compared to the high torsion stiffness case (Figure 5.17(b)). However, 
the opposite behavior occurs over the shaft speed range [Q « 2400 - 2650 RPM]. In this range the 
response jumps from the 2T p , limit-cycle to the higher amplitude 3 T p limit-cycle. Thus, the nonlinear 
dynamic mesh phasing effect between the upstream and downstream gearboxes, which is strongly related 
by the shaft torsional stiffness, is seen to have both beneficial and negative effects on the variable RPM 
torsional vibration response amplitudes. 
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Figure 5.16. — Tooth impact velocity at steady-state versus shaft 
speed; h gap = 20 //m; k t = 6000 Nm. 



Figure 5.17. — Gear mesh vibration during shaft spindown from 
[Q = 3000 to 0 RPM]; with h gap = 20 //m for; (a) k t = 15000 Nm; 
(b) kt = 6000 Nm. 


5.6 Summary and Conclusions 

This section studies the nonlinear interactions between two gearboxes coupled by a torsionally 
flexible shaft. Such configurations are common in rotorcraft propulsion systems including tailrotor and 
cross-shafting drivelines. Shaft flexibility induces dynamic phasing between two coupled gearboxes. This 
dynamic phasing together with the kinematic dependence of the gear tooth mesh stiffness on the gear 
rotation induces a nonlinear dynamic interaction not explored in previous single-stage or rigid-shaft gear 
system investigations. The nonlinear time -varying equations-of-motion of a dual-spur gearbox/flexible 
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shaft system are derived via an assumed mode technique. Steady-state harmonic and subharmonic 
solutions are obtained via a harmonic balance and arc-length continuation strategy. Under zero backlash 
conditions, the dynamic mesh phase produces a softening-type nonlinearity in the resonance curves. 
Furthermore, under chattering conditions, increasing the shaft compliance tends to destabilize the single 
and double sided contact response amplitudes. The equations and analysis presented for the dual-gearbox 
shaft system demonstrate the importance of including nonlinear mesh phase interactions in multigearbox 
drive systems. 

The equations and analysis presented for the dual-gearbox shaft system demonstrates the importance 
of including nonlinear mesh phase interactions in multigearbox drive systems. It is found that, similar to 
backlash effects, the mesh phase nonlinearity produces a softening-type resonance behavior. Furthermore, 
for higher load torque levels where gear chattering vibrations occur, reducing the intermediate shaft 
flexibility has the effect of destabilizing 2 T p solutions which naturally reduce the chattering response 
amplitudes over certain RPM ranges. This study also demonstrates the power and effectiveness of the 
harmonic balance and continuation approach for finding steady-state solutions to multiple degree-of- 
freedom nonlinear time- varying systems. One key factor is the location of the system structural modes in 
relation to the gear mesh harmonics. For variable speed operation, it is necessary to minimize the 
frequency crossings during spinup and spindown. When integer multiples of gear mesh frequency become 
close to structural modes involving gear mesh deflections chattering can occur. This chattering is 
exacerbated by the static load torque level. The analysis and results give insights into the relative 
importance of dynamic mesh phasing phenomena in multistage gear/shaft systems and should provide 
some design guidelines for flexible shaft/gearboxes systems. 
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6.0 Two-Speed Dual Clutch Transmission 

6.1 Introduction 

Aircraft sizing design studies have identified that a two-speed rotor configuration is required (Refs. 1 
and 104). Rotor speed at cruise is required to be reduced, possibly down to 50 percent of hover speed, to 
keep rotor blade advancing tip speeds reasonable. One early study of a multispeed rotorcraft drive system 
consisted of a high speed traction drive variator and a planetary differential gear unit (Ref. 4). A unit was 
tested and performed as designed. However, the mechanical efficiency was lower than that of a 
conventional drive system. Also, the use of traction drives for a main drive path is not an accepted practice 
for U.S. rotorcraft manufacturers. Previous studies have been sponsored by NASA to investigate multispeed 
drive systems (Refs. 3 and 105). These studies looked at earlier tiltrotor applications where 15 percent speed 
reduction for cruise was required. Dual path configurations and compound planetary configurations looked 
promising. However, a major concern of the shifting process as well as drive system weight penalties 
resulted from the studies. An electro-mechanical infinitely variable transmission, comprising a pair of 
planetary trains interconnected with two electric machines and clutches, has been proposed (Ref. 106). 
Again, the mechanical efficiency suffered as well as the added complexity of two planetary systems. A 
unique concept called a pericyclic continuously variable speed transmission is being investigated (Ref. 107). 
It uses pericyclic kinematics to achieve single speed reduction ratios between 1.05:1 to 50:1 with variable 
speed capability in one hardware stage. This concept, however, had yet to be tested. In summary, the results 
from the previous studies indicate that the incorporation of multispeed concepts in rotorcraft application is 
not a trivial process and adds complexity and weight. Thus, further research and development in this area is 
required, such as described in (Ref. 108). In addition, dynamic modeling of multispeed drive systems can 
help define the characteristics and limitations of the drive systems and the shifting process. 

Since very few multispeed rotorcraft studies have been published, very little work on dynamic 
modeling for multispeed rotorcraft drive systems has been done. There has been, however, numerous 
publications on dynamic modeling of automotive drive systems. Much of the work addresses the 
automatic transmission shifting control with emphasis on maximizing fuel economy and maintaining 
performance (Refs. 109 to 115). For these, the typical system being modeled is the gasoline combustion 
engine, torque convertor, multispeed gearbox, tires and vehicle dynamics, as well as throttle and shifting 
controls. Dynamic speed, gear ratio, and torque are amongst the parameters predicted. This strategy has 
been applied to hybrid powertrains (Refs. 116 and 1 17), continuously variable transmissions (Refs. 1 18 to 
121), and even heavy trucks and tanks for the military (Refs. 122 to 126). The modeling in the automotive 
field has had much success in improving drive train development. The approach and some modeling tools 
can be extended to the rotorcraft arena (Refs. 127 and 128). However, specifics of integration of the gas 
turbine engine, multispeed drive systems, clutches, and the dynamics of the rotor and rotorcraft vehicle 
have yet to been addressed in the rotorcraft application. 

6.2 Two-Speed Dual Clutch Transmission 

In order to implement the variable speed rotorcraft driveline confutations, a two-speed Dual Clutch 
Transmission (DCT) similar to that considered by Kish, 2002 is utilized in this study. Figure 6. 1 shows the 
two-speed transmission considered in (Ref. 3) and Figure 6.2 shows the two-speed DCT considered in this 
study. 

One key advantage of the DCT gearbox configuration is that it provides speed change via gear shifting 
though a split-path differential planetary arrangements. Here, in Figure 6.2, the high-speed ratio is 
accomplished by disengaging a ring clutch and engaging a set of planet clutches. For the low-speed ratio, 
the planet clutches are disengaged and the ring clutch is applied to lock the ring gear. In the low speed case, 
the power is split between two parallel paths which has the advantages of allowing an overall lower 
transmission weight. The corresponding kinematic diagram of the two-speed DCT in Figure 6.2 is shown in 
Figure 6.3. 
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Figure 6.2. — Two-Speed DCT, current study. 
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Figure 6.3. — Two-Speed DCT, kinematic diagram. 


The rotation speed kinematics of the DCT are as follows 
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where Q, is input shaft speed, Q cg i, Q g 2 , and Cl cg3 are control gear speeds, Q s , Q r and Q p are sun, ring and 
planet gear speeds, Q cr is the carrier speed and Cl 0 is the output shaft speed, and where Z h Z cgU Z cg2 , Z cg3 , 
Z r , Z s and Z p are the number of gear teeth respectively. Here the carrier is the DCT output which is related 
to the input by 


~ n DCT^i (2) 

Where the n DC r is the total DCT reduction ratio. The high ratio, n DC r,hi g h is achieved by engaging the 
multiple clutch disks that couple the planet shafts (clutch 1). 


n DCT , high 


ZiZcg 2 Z r 
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Furthermore, low ratio, n DC T,iow , is achieved by fully engaging, the cylindrical ring clutch (clutch 2) which 
clamps the ring gear. 


n DCT,low 


Z, 

2 (Z p +Z s ) 


(3-b) 


Clutch engagement is controlled by applied pressures, p dx and p d2 , normal to the clutch fiction 
surfaces. Due to the DCT kinematics, the only condition under which both sets of clutches can be 
engaged and locked simultaneously corresponds to the zero speed or fully stalled condition. When both 
clutches are applied simultaneously there is a negative recirculating power ratio and clutches act as a 
braking system. Thus, during steady-state constant speed operation, only one clutch is fully engaged and 
locked (with maximum applied clutch pressure p di =p dm ax ) while the other clutch remains open (with no 
applied clutch pressure p dj = 0). The parameters of the two-speed DCT utilized in this study are 
summarized in Table 6.1. 


TABLE 6. 1 . — TWO-SPEED DUAL CLUTCH TRANSMISSION PARAMETERS 
Parameter Value 

Number of gear teeth Z,=40, Z cgX = 42, Z cg2 = 29, Z cg3 = 52, Z r =74, Z= 14, Z= 46 

Gear radii, m R,=0.1, R cgl =0.106, i? cg2 =0.074, R cg 3 = 0.132, R=0A9, R p = 0.036, 7^=0. 12, R c = 0.152 

Rotational inertias, kg m 2 .7=0.064, «/ cgl =0.078, J cg2 = 0.039, J cg3 =0.402, J= 1.766, J p = 0.002, .7=0.264, .7^=0.848 

Number of planets N p = 8 

Planet mass, kg m p = 3.576 

Planet clutch radius (clutch 1), m R clx = 0.238 

Number of planet clutches (clutch 1) N cg = 2 

Ring clutch dimensions (clutch 2), m R c n = 0.197, L cl2 = 0.1 

Clutch friction coefficient ju d = 0.45 

Clutch viscous torque parameter, N m sec r d = 1 x 1 0 2 


Based on these parameters, the low ratio is n DCT j ow = 0.383 and the high ratio is n DCT high = 0.71 1 
which means about 47 percent speed reduction from the high to the low ratio. 

6.3 Dual Clutch Transmission Dynamics Model 

To account for the DCTs in the driveline simulation, a rigid body dynamics model is developed which 
accounts for the DCT kinematics, the gear rotational inertias, and the nonlinear clutch friction torques. 

The equations of motion have the form 


M = F(£2, p c . n , p c / 2 ) (4-a) 

with corresponding angular velocity degrees-of-freedom 

n(0 = [n,-(0. n r (f)] r ( 4_b ) 

consisting of the DCT input and ring gear rotation speeds. Also, the system inertia matrix is 
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Where P and P, are input and output shaft the torques on the DCT and where Pn and P c a are the 
transmitted clutch torques which are computed based on a combined Coulomb/viscous model (Refs. 127, 
129 and 130 


fell 
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\pei 2 Pd IxRcn L c n sgn(Q r ) when |Q r | > s 


(sticking) 
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Here ju d and r d are clutch friction and viscous torque coefficients (Ref 129). Note, clutch 1 is a disk-type 
clutch therefore r dl is a function of disk radius R cn , and clutch 2 is a cylindrical ring clutch thus r d2 
depends on clutch radius, R c 12 , and axial length L cll . This nonlinear torque model accounts for the 
possibility of clutch sticking (lockup) and slipping. 

In this clutch model, (5), the individual clutch torques, r di , transition between viscous (sticking) 
mode and Coulomb friction (slipping) mode based on the clutch slip-speed threshold which is evaluated 
at each simulation time-step. The slip-speed threshold parameter, £, is a small positive number in the 
simulation. If the clutch slip-speed exceeds £then the torque transitions from sticking to slipping mode. 
Likewise, if clutch slip-speed drops below £then the clutch torque transitions back to sticking mode. By 
proper selection of s and r ch the viscous mode numerically approximates clutch lockup without no need 
for explicit kinematic constraints or model DOF reduction (Ref. 127). The DCT system parameters used 
in this investigation are summarized in Table 6.1. Also, it is important to note that the clutch torque, T di , 

depends on the applied pressure clutch p di only during the slipping mode. 

Figure 6.4 shows the Matlab Simulink block diagram representation of the dual clutch DCT system 
utilized in the comprehensive simulation models presented in Section 7.0. Here the saturation blocks are 
utilized to allow the transition from sticking to slipping torque modes. 
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Figure 6.4. — Two-Speed DCT dynamics model. 
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7.0 Comprehensive Variable Speed Rotor craft Propulsion System Model 
and Simulation 

7.1 Introduction 

One of the main goals of this research project is to develop a comprehensive first principles-based 
drive-system modeling tool to predict the complex transient dynamic response of variable speed rotorcraft 
propulsion systems including engine, transmission, clutch, cross-shafting, and rotor dynamic interactions. 
By accounting for the various propulsion subsystem dynamic interactions, this new Comprehensive 
Variable Speed Rotorcraft Propulsion System Modeling (CVSRPM) tool enables performance prediction 
to give future design insights and enable system level understanding. 

Furthermore, another goal is to exercise this newly developed computational tool to explore the 
behavior of various variable ratio rotorcraft drive-system configurations and ratio changing strategies 
under simulated operational conditions. To address these goals, a comprehensive rotorcraft propulsion 
system simulation tool based on the subsystem models developed in Sections 2.0 to 6.0 is constructed. 
Simulations are then performed for several different drive system configurations to evaluate their variable 
speed performance characteristics. 

7.2 Variable Speed Rotorcraft Propulsion System Model 

The comprehensive dynamic simulation tool assembled in this project is based combining several of 
the key rotorcraft driveline subsystem models developed for in Sections 2.0 to 6.0. Figure 7.1 shows as 
schematic of the overall CVSRPM architecture along with the various subsystem model interaction paths. 
For reference, Figure 7.1 also gives the corresponding section number of the particular subsystem model 
when appropriate. Specifically, the main elements and assumptions included in this comprehensive, 
variable speed, drive system simulation tool are as follows 

• Two-spool gas turbine engine transient dynamics model and closed-loop PID fuel controller to 
regulate power turbine speed (Section 6.0). 

• Dual-Clutch Transmission (DCT) gearbox model with stick-slip friction and clutch lockup effects 
(from Section 6.0). 

• Simplified main rotor model based on momentum theory (from Section 2.0) 

• Torsionally flexible shafting and cross-shafts (misalignment and bending not considered) 

• All rigid gearing (tooth mesh compliance and mesh stiffness variations neglected) 

• A two DOF flight dynamics model for fuselage drag, vehicle forward flight and elevation 
dynamics with vehicle states [ x, x, y, y ]. 

• A PID based main rotor blade pitch controller to regulate rotor thrust to track a desired forward 
speed setpoint (controls main rotor pitch 9 pitch ) 

• A PID based wing elevon pitch controller to maintain level flight or track a desire climb rate 
setpoint (controls wing flap angle 6 wing ) 

The comprehensive dynamic simulation code is implemented under the Matlab Simulink software 
environment. Here, each subsystem block of the CVSRPM code as shown in Figure 7.1 is detailed in the 
Appendix. 
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Figure 7.1. — Comprehensive Variable Speed Rotorcraft Propulsion System Model. 

Due to the inherently high frequency nature of the gear mesh stiffness variation and backlash effects 
(as is considered in Section 5.0), difficulties related to excessive simulation times and storage space were 
encountered when these effects were included in the comprehensive simulation model. For example, the 
first fundament drive system torsion mode is typically on the order of 1 to 5 Hz while the gear mesh 
stiffness variations may be on the order of 5000 Hz, thus requiring excessively small numerical 
integration time-steps compared to the overall drive-system maneuver timescales. Here, the decision was 
made to neglect the gear mesh effects in the comprehensive simulation model and consider rigid gearing 
to reduce the simulation runtimes and data storage requirements. Apart from these issues, the simulations 
can be run within a “reasonable” about of time on a modern PC or notebook computer. For example, for a 
typical dual-path drive system (tiltrotor) performing an upshift followed by a downshift maneuver over a 
simulated time span of 280 sec, the simulation can be can be completed in approximately 15 min of CPU 
execution time. In the subsequent sections, several two-speed gear change maneuvers for a single-path 
and dual-path rotorcraft drive system configuration is explored using the CVSRPM code. 

7.3 Two-Speed Shift Rotorcraft Drive-System Case Studies 

In this section, the newly developed CVSRPM code is used to investigate the variable speed 
operational behavior of two different rotorcraft drive system configurations each equipped with two speed 
dual clutch transmissions (see Section 4.0). Here both the conventional helicopter and the tiltrotor drive- 
system configurations are considered. In each case, the system executes a fully trimmed simulated in- 
flight maneuver as regulated by the three different closed-loop control systems (i.e., the engine fuel 
control, blade pitch forward speed control, and the auxiliary wing lifting surface controller). During flight, 
open-loop clutch shift commands are applied to the dual-clutch gearbox transmission(s) and the resulting 
variable speed system response is computed. 
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7.3.1 Two-Speed Helicopter Drive System 

Figure 7.2 gives a schematic of the two-speed helicopter drive system simulation model implemented 
with the CVSRPM code. 

In this helicopter drive system there are three reduction stages from the power turbine speed Cl pt to the 
main rotor speed Cl MR . 


^ MR ~ n f n DCT n i^pt (1) 

Here n f is the initial reduction ratio, n DC r is the DCT reduction ratio (see Section 4.0) and n f is the final 
reduction. 

In this analysis, the variable rotor speed propulsion system response is computed for the helicopter in 
a hover condition. The free-body-diagram of the simplified hover model is shown in Figure 7.3 and the 
hover flight dynamics model is given by Equation (2). 

m v y = F T (Fl MR 

•> @ pitch ? n)-^vg with V v =y (2) 

Where m v is vehicle mass, g is acceleration of gravity, y is vertical position coordinate of the vehicle mass 
center, v v is vehicle vertical speed, and F T is main rotor thrust. Here F T is a function of main rotor speed, 
Cl MR , main rotor blade pitch, 6 pitch , and v v . The main rotor thrust computation is shown in Section 2.5. In 
this study, the parameters for the hypothetical two-speed rotor helicopter system are based on a CH-53 
helicopter platform (Ref. 131). All of the helicopter system parameters used in this study are given in 
Table 2.3, Table 4.1, Table 4.3, Table 6.1, Table 7.1, and Table 7.2. 




T ™ v g 

Figure 7.3. — Helicopter hover flight dynamics 
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TABLE 7. 1 .—HELICOPTER PARAMETERS 


Parameter Value 

Vehicle mass, m v 20,865 kg 

Initial ratio, rij 0.5347 

Final ratio, n d 0.04 



Figure 7.4. — Blade pitch control block diagram. 


In order to maintain hover with a desired vertical speed setpoint, v v , blade pitch is regulated with a 
PID-based blade pitch controller of the form 

^ pitch ~ kp p e v + k Ip J e v dt + k dp e v with e v = v Vj- - v v (3) 

With blade pitch control gains k Pp , k Ip and k dp , to regulate the setpoint speed error e v . It is also 

necessary to account for the physical blade pitch limits O p i tc h m in and @ P itch max which is accomplished 
through use of a saturation block and back-calculating gain, k bp , to prevent integrator windup. Figure 7.4 

shows a block diagram of the blade pitch control. 

The gains k Pp , k Ip and k dp were selected by trial-and-error until a stable closed-loop system was 

achieved. The blade pitch control system parameters are shown in Table 7.2. 


TABLE 7.2.— HELICOPTER BLADE PITCH 
CONTROL SYSTEM PARAMETERS 


Parameter 

Value 

Proportional gain, k p 

0.004 rad-sec/m 

Integral gain, kj 

0.005 rad/m 

Derivative gain, k d 

4. 0x1 0' 5 rad-sec 2 /m 

Back-calculating gain, k b 

100 sec' 1 

Minimum blade pitch, 6 pitch 

0° 

Maximum blade pitch, 6 pitch 

17° 
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Figure 7.5. — Main rotor downshift/upshift during hover. 


The first variable rotor speed mission profile studied is a main rotor downshift/upshift procedure 
during hover conditions. This variable rotor RPM maneuver consists of a downshift followed by an 
upshift, all during hover with a zero vertical speed setpoint, v Vset = 0 kn, as shown in Figure 7.5. 

The resulting system response is detailed in Figure 7.6 to Figure 7.9. Here the helicopter is initially 
operating at 100 percent rotor speed, then, at t = 40 sec, the downshift is initialed. The high-speed clutch 
axial engagement pressure is commanded to ramp-down while the low-speed clutch pressure (originally 
disengaged) is simultaneously commanded to ramp-up. The pressures are commanded to keep increasing 
(or decreasing) until clutch becomes fully locked (or fully disengaged). The clutch command engagement 
pressures are plotted in Figure 7.6(c). 

Figure 7.6(b) shows corresponding response of the main rotor speed which shows an essentially linear 
decrease in rotor RPM until the low-speed clutch becomes locked. Then a very slight RPM dip is observed 
immediately after lockup which is due to system torsional compliances. Furthermore, Figure 7.6(a) shows 
the engine power turbine speed which is strongly influenced the choice of PID fuel control gains (see 
Section 4.0). Next at t = 80 sec, the upshift command is initiated. Here the process is essentially the revere 
of the downshift and the corresponding rotor speed and engine speed responses are also well behaved. 

Still more insight is gained by examining Figure 7.7, which shows the system torque responses. Here 
it is observed that the main rotor torque (Figure 7.7(a)) actually increases as a result of the lower rotor 
RPM. At first pass this seems counter intuitive, however, the rotor torque increase is a result of the pitch 
controller trying to maintain a constant rotor thrust in order to maintain the hover. As the rotor RPM 
drops, the pitch controller must add more collective blade pitch to generate the thrust, thereby increasing 
the rotor drag torque. The blade pitch control input is shown in Figure 7.9(c). 

Next the, engine output torque, shown in (Figure 7.7(b)), is also seen to significantly rise during both 
the downshift and the upshift events. Due to the magnitude of the engine torque increase, this cannot be 
explained by the effect of the main rotor blade pitch control. Rather, it is the result of the opposing split path 
power flows inside the DCT gearbox. For the applied shift command inputs (Figure 7.6(c)), both the high 
and low speed clutches are simultaneously engaged and slipping during up and downshifts. The resulting re- 
circulating power must be counteracted by the engines. However, it is this same effect which enables the 
smooth, linear rotor speed transitions observed in Figure 7.6(b). Without some level of clutch overlap and 
power re-circulation, the rotor RPM may go into a neutral condition and lose speed rapidly. This is further 
illustrated by observing the transmitted clutch torques in Figure 7.7(c). Here, also high frequency stick-slip 
torque oscillations are observed in the low-speed clutch during its disengagement. Also Figure 7.8 shows the 
corresponding clutch power dissipated during the shift. Here again both clutches dissipate heat during the 
transitions but then, of course, the power dissipation is zero after it either becomes fully locked or fully 
disengaged. It is also interesting to observe that the power dissipated during the downshift is significantly 
more than during the rotor upshift. One explanation for this is that a significant portion of the main rotor’s 
kinetic energy is being dissipated by the clutches during the rotor speed downshift. Likewise, during upshift, 
power is flowing into the rotor. Here again, it is all down to choice of shift command inputs or shift strategy. 
If the rotor was allowed to go into a pure neural condition it would naturally execute a decrease in RPM due 
the rotor drag torque, with zero power dissipated in the clutches. However, the designer must find the 
correct balance between clutch overlap (with smooth ratio changes) and a neutral drop (no clutch heating 
but less rotor RPM control during shifts). 
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Finally Figure 7.9 summarizes the corresponding helicopter flight parameters during the hover 
downshift/upshift simulation. Here, it is seen that the rotor downshift and upshift causes minimal 
disturbances to the vertical speed and accelerations (Figure 7.9(b)). 
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Figure 7.6. — Two-speed helicopter drive-system response for 
downshift/upshift during hover; (a) engine power turbine 
response; (b) main rotor speed; (c) clutch engagement 
pressures 
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Figure 7.7. — Two-speed helicopter drive-system response for 
downshift/upshift during hover; (a) main rotor torque; (b) engine 
output torque; (c) clutch torques. 
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Figure 7.8. — Two-speed helicopter drive-system response for 
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power dissipation. 
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Figure 7.9. — Two-speed helicopter drive-system response for downshift/upshift 
during hover; (a) vertical speed; (b) vertical acceleration; (c) blade pitch input, 
(d) rotor induced flow. 


7.3.2 Two-Speed Tiltrotor Driveline System 

As mentioned in References 1,3, and 104, one of the potential beneficiaries of variable speed rotor 
drive-system technologies is in the area of Heavy-Lift helicopter and tiltrotor configurations. Therefore, 
the next case explored using the CVSRPM simulation code is a tiltrotor configuration based on the NASA 
LCTR-2 Civil Tiltrotor shown in Figure 7.10 (Ref. 104). 

Specifically, Figure 7.1 1 is the schematic of the two-speed tiltrotor drive-system simulation model 
implemented with the CVSRPM code. 
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Figure 7.11. — Two-speed tiltrotor dual path driveline simulation model. 


In each branch of the dual path tiltrotor driveline there are three reduction stages from the power 
turbine to the main rotor. Here n t is the initial reduction, n DC r is the DCT reduction ratio (see Section 6.0) 
and rif is the final reduction stage. Furthermore, as shown in Figure 7.1 1, the left and right branches of the 
drive-system are linked via a cross-shaft. In this analysis, the cross-shaft is considered to be torsionally 
flexible with torsional stiffness k cs . This dual path driveline is implemented in the CVSRPM code by 
utilizing two identical single path drivelines similar to the conventional helicopter and then linking them 
via a torsional spring element to simulate the cross-shaft. 
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In this study, variable rotor speed operation is analyzed for the tiltrotor operating in forward flight 
airplane mode as shown in Figure 7.12. Furthermore, the forward flight dynamics model of the tiltrotor is 
given by Equation (4). 


m v x - 2F T (Cl MR 

? @ pitch ? ^ c ) F D (Vc ? @wing ) 

Ky = F L ( Owing > V c ) - m vg (4) 

V c =X 

Where x and y are horizontal and vertical vehicle position coordinates. F T is main rotor thrust, F D is 
vehicle aerodynamic drag and F L is the total wing lift force. In the tiltrotor case, F T is computed according 
to Section 2.5 and is a function of main rotor speed, Fl MR , main rotor blade pitch, 0 pitch , and forward speed 
v c . The lift a drag forces are computed as 

F l = P^A w c l v c 2 and F D =^-A w c D v c 2 (5) 

Where p air is air density, A w is total wing lift area and c L and c D are wing lift and drag coefficients 
respectively. Finally, the wing lift coefficient is computed in terms of the wing flap pitch 6 wing as 

~ C Lapwing (^) 

Where c La is wing lift coefficient slope. The tiltrotor system parameters assumed in this investigation are 
based on the NASA LCTR-2 (Ref. 104). All of the tiltrotor system parameters used in this study are given 
in Table 2.4, Table 4.1, Table 4.3, Table 6.1, Table 7.3, Table 7.4, and Table 7.5. 



Figure 7.12. — Tiltrotor forward flight airplane mode flight dynamics. 
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TABLE 7.3.— TILTROTOR PARAMETERS 


Parameter Value 

Vehicle mass, m v 48762 kg 

Initial ratio, 0.5347 

Final ratio, n f 0.04 

Cross-shaft torsional stiffness, k cs 21219 N m 

Wing lift coefficient slope, c La 12.578 rad' 1 

Vehicle drag coefficient, c D 0.054 

Wing lift area, A w 93.04 m 2 


In order to maintain forward flight two PID control loops are utilized. The first control loop adjusts 
main rotor pitch 0 pitch to maintain a forward speed setpoint v Cset See Figure 7.4 for a block diagram. 

Q pitch = k Pp e v + k Ip J e v dt + k dp e v with e v = v Csa - v c (7) 

The other control loop adjusts wing flap angle 6 wing to maintain a vertical speed setpoint v Vset . For level 

fli g ht = 0. 


Gyring = k iK e v + k i„ ]e v dt + k dw e v with e v = v Vsa - v v (8) 

Similar to the rotor pitch control, the wing flap control also utilizes a saturation block and a back- 
calculating gain, k bw , to account for the wing pitch limits. Figure 7.13 shows the wing pitch control block 
diagram. 

The wing and pitch control gains were selected by trial-and-error until a stable closed-loop system 
was achieved. The blade and wing pitch control system parameters for the tiltrotor are given in Table 7.4 
and Table 7.5. 
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TABLE 7.4.— TILTROTOR BLADE PITCH 


CONTROL SYSTEM PARAMETERS 


Parameter 

Value 

Proportional gain, k p 

0.004 rad-sec/m 

Integral gain, k d 

0.005 rad/m 

Derivative gain, k d 

4.0* 10' 5 rad-sec 2 /m 

Back-calculating gain, k b 

100 sec' 1 

Minimum blade pitch, 0 pitch 

0° 

Maximum blade pitch, 6 pitch 

80° 


TABLE 7.5.— TILTROTOR WING PITCH 
CONTROL SYSTEM PARAMETERS 


Parameter Value 

Proportional gain, k Pw 0.1 rad-sec/m 

Integral gain, k Iw 0.05 rad/m 

Derivative gain, k dw 0.02 rad-sec 2 /m 

Back-calculating gain, k bw 1.0 sec' 1 

Minimum wing pitch, 0 wingmin 0° 

Maximum wing pitch, 0 wingmm 10° 
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Figure 7.14. — Variable rotor speed tiltrotor mission profiles; (a) steady forward cruise downshift/upshift; (b) variable 
forward cruise downshift/upshift. 


Another factor which becomes important in coupled multipath drive-systems is component 
dissimilarly. In the tiltrotor case (Figure 7.14), if each component is identical with symmetric rotor loads, 
then the cross-shaft torque is zero and there will be no interactions between the two drive-system paths. In 
the case of the two-speed tiltrotor configuration, one unique and potentially significant source of 
component dissimilarly is with the clutches. Differences in clutch frictional characteristics between the 
left and right-side tiltrotor driveline branches could result in clutch torque interactions across the cross- 
shaft during shift events. To assess this possibility, a 1 percent dissimilarly between the left and right side 
DCT clutch friction coefficients /u d is assumed in the analysis. 

The simulated mission profiles considered in this investigation for the variable rotor speed tiltrotor 
are shown in Figure 7.14. 
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The first scenario (Figure 7.14(a)) is an upshift followed by a downshift while maintaining level flight 
with constant forward setpoint speed. The second scenario (Figure 7.14(b)) is an upshift and a downshift 
in level flight with variable forward speed. This is the mission profile discussed in (Ref. 104) (referred to 
as the slowed rotor cruise scenario) for the LCTR-2 when it is operating in airplane mode. In this 
situation, in order to reduce the blade tip Mach number and associated losses, the rotor is slowed to permit 
higher a forward flight cruising speed. In this case, at the start of the speed change maneuver, the rotor 
speed is reduced from Cl MR 190 to 102 RPM (a 47 percent reduction) while the forward speed setpoint v cset 
is simultaneously commanded to increase from 150 to 300 kn. 

These mission profiles are simulated under two different shift-control strategies; (a) one using so- 
called Paralleled-Shift Control (PSC) and the other (b) using so-called the Sequential-Shift Control 
(SSC). Under the PSC strategy both the left and right side Dual Clutch Transmissions (DCTs) in Figure 
7.1 1 are shifted simultaneously while under the SSC strategy the left and right side DCTs are shifted one 
after another to avoid shifting under full load. Both the PSC and the SSC simulation results are presented 
in the next two subsections. 

7.3.2. 1 Parallel-Shift Control (PSC) 

During PSC shifting both left and right side DCTs are commanded to shift simultaneously. To 
execute the upshift in each DCT, the low-speed clutch pressure is steadily decreased to zero while the 
high-speed clutch pressure is steady increased from zero to some finite value sufficient to cause full 
lockup of the high speed clutch. The downshift is executed in the same manner, just in reverse. As 
explained in the previous section, this process results in a cross-over from low to high ratio (or from high 
to low) without passing through a neutral condition which may be detrimental to rotor speed. However, 
during this transient handoff process when both high and low speed clutches are slipping under load, a re- 
circulating power loop is setup resulting in some amount of power being dissipated into heat. The DCTs 
must be designed to thermally cope with and absorb this transient fictional power without overheating. 

7.3.2. 1.1 Steady Forward Cruise PSC Downshift/Upshift 

Figure 7.15 to Figure 7.21 give the simulation results for the two-speed DCT tiltrotor executing a 
PSC downshift followed by a PSC upshift in steady forward cruise. By examining Figure 7.15(a) and (b), 
it is apparent that main rotor speed is successfully lowered with little deviation in forward velocity. Figure 
7.15(c) and (d) show the open-loop DCT clutch pressure command inputs. Also, the engine power turbine 
speed is shown in Figure 7.15(e) and the corresponding combustion chamber fuel rates commanded by 
the closed-loop engine fuel control are shown in Figure 7.15(f). 

Furthermore, the flight dynamics are summarized in Figure 7.16. Figure 7.16(a) shows that both the 
horizontal and vertical vehicle accelerations are minimal. Thus, no forward speed or altitude is lost during 
the PSC downshift or upshift. In order to maintain the constant forward setpoint speed under changing 
rotor speed, the closed-loop blade pitch controller automatically adjusts the main rotor collective pitch to 
keep a constant thrust. Specifically, as seen in Figure 7.16(b), collective blade pitch increases as a result 
of the downshift and decreases again in reaction to the upshift. The corresponding rotor induced velocity 
is seen in Figure 7.16(c). Finally, as seen in Figure 7.16(d), under this constant thrust condition, there is 
no appreciable change in wing flap angle. 

The drive system torques are summarized in Figure 7.17. As observed in the engine fuel rate response 
in Figure 7.15(f) as well in the engine torque response in Figure 7.17(a) the engine torque increases 
during both the PSC downshift and the upshift. At first this seems counter intuitive. It may be obvious 
why engine torque increases during the upshift phase since extra torque is requited to accelerate the rotor 
however, it is not initially obvious why engine torque also increasers during downshift. The reason for 
this is due to the re-circulated power loop setup during the PSC shifting as described previously. During 
the handoff from the disengaging to the engaging clutch there is no neutral condition, hence when both 
clutches are slipping they are essentially acting as a brake which the engines must resist. 
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Figure 7.15. — Steady cruise PSC downshift/upshift; (a) forward speed, (b) rotor speed, (c) left 
DCT shift commands, (d) right DCT shift commands, (e) power turbine speeds, (f) fuel 
rates. 
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Figure 7.16. — Steady cruise PSC downshift/upshift — flight dynamics; (a) vehicle 
acceleration, (b) collective pitch, (c) induced velocity, (d) wing flap. 
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Figure 7.17. — Steady cruise PSC downshift/upshift — torque response; 
(a) engine torques, (b) cross-shaft torque, (c) rotor torque. 
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Figure 7.18. — Steady cruise PSC downshift/upshift — clutch response; (a) clutch 
torques, (b) power dissipation (clutch heating). 
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Figure 7.19. — Steady cruise PSC downshift/upshift — engine response; 

(a) compressor speed, (b) mass flow, (c) power turbine inlet temperature. 
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This is further seen in Figure 7.18 which shows the high-speed and low-speed clutch torques of the 
DCT and the corresponding frictional power dissipation of each clutch. This dissipated power essentially 
goes directly into clutch heating. Based on these values it is seems that the PSC shifting strategy is not a 
viable scenario for this size rotorcraft due to the high power levels which must be thermally absorbed. 
Also, Figure 7.17(c) shows the increasing rotor torque which is a result from the increasing collective 
pitch (see Figure 7.16(b)). 

Furthermore, the internal details of the two-spool gas turbine engine dynamic response are shown in 
Figure 7.19 to Figure 7.21. Figure 7.19(a) to (c) shows compressor rotational speed, air mass flow rate 
and turbine inlet temperature respectively. Also, Figure 7.20 shows the compress running line in terms 
pressure ratio and mass flow. Here it is noted that the entire operating line remains below the surge 
compressor surge margin during both PSC downshift and upshift at constant forward speed. 

In addition, since Figure 7.20 does not show temporal information, the engine surge margin (SM) 
versus time is examined in Figure 7.21. Here, the working definition of SM used is 

SM = x 1 00% (9) 

P r c , surge 

Therefore SM > 0 indicates no compressor surging. 


compressor 
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Figure 7.20. — Steady cruise PSC downshift/upshift — two-spool gas 
turbine engine transient running lines. 


Downshift Upshift 



Figure 7.21. — Steady cruise PSC downshift/upshift — engine compressor 
surge margin. 
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Finally, the effect of DCT shift-rate is examined under the PSC upshift/downshift scenario. Table 7.6 
summarizes the four shift-rate cases examined and Figure 7.22 and Figure 7.23 show the corresponding 
engine compressor pressure ratio dynamics and clutch frictional power dissipation for each case. Figure 
7.22 and Figure 7.23 illustrate the tradeoff between frictional power dissipation (which corresponds 
directly to clutch heating) and engine compressor surge. In the PSC shifting strategy slower shift rates are 
less prone to cause engine compressor surge but also result in increased clutch heating. 


TABLE 7.6.— CLUTCH PRESSURE ENGAGEMENT AND DISENGAGEMENT RATES FOR PSC 


Case 

Upshift 

Downshift 

Slow 

Clutch 1 engagement rate, p dx = 2 psi/sec 
Clutch 2 disengagement rate, p d2 = -2 psi/sec 

Clutch 1 disengagement rate, p dx = -2 psi/sec 
Clutch 2 engagement rate, p d2 = 2 psi/sec 

Medium- Slow 

Clutch 1 engagement rate, p dx = 4 psi/sec 
Clutch 2 disengagement rate, p d2 = -4 psi/sec 

Clutch 1 disengagement rate, p dx = -4 psi/sec 
Clutch 2 engagement rate, p d2 = 4 psi/sec 

Medium-Fast 

Clutch 1 engagement rate, p dx = 8 psi/sec 
Clutch 2 disengagement rate, p d2 = -8 psi/sec 

Clutch 1 disengagement rate, p dx = -8 psi/sec 
Clutch 2 engagement rate, p d2 = 8 psi/sec 

Fast 

Clutch 1 engagement rate, p dx =16 psi/sec 
Clutch 2 disengagement rate, p d2 = -16 psi/sec 

Clutch 1 disengagement rate, p dx = -16 psi/sec 
Clutch 2 engagement rate, p d2 =16 psi/sec 
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Figure 7.22. — Steady cruise PSC downshift/upshift — compressor pressure ratios versus mass flow rate for 
four shift-rate cases; (a) slow, (b) medium-slow, (c) medium-fast, (d) fast. 
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Figure 7.23. — Steady cruise PSC downshift/upshift — clutch frictional power 
dissipation for four shift-rate cases; (a) slow, (b) medium-slow, (c) medium- 
fast, (d) fast. 


7.3.2. 1.2 Variable Forward Speed Cruise PSC Downshift/Upshift 

Figure 7.24 to Figure 7.30 give the simulation results for the two-speed DCT tiltrotor executing a 
PSC downshift followed by a PSC upshift during the variable forward speed cruise condition depicted in 
Figure 7.14(b). By examining Figure 7.24(a) to (b), starting at t = 20 sec it seen that the forward cruising 
speed is increased from 150 to 300 kn while the rotor speed is simultaneously downshifted from about 
190 to 102 RPM. Then, starting at t = 105 sec it seen that the forward speed is increased back to 300 kn 
while the rotor is simultaneously upshifted back to 190 RPM. Figure 7.24(c) to (d) show the open-loop 
DCT clutch pressure command inputs. 

Also, the engine power turbine speed is shown in Figure 7.24(e) and the corresponding combustion 
chamber fuel rates commanded by the closed-loop engine fuel control are shown in Figure 7.24(f). During 
the rotor downshift/forward speed acceleration phase Figure 7.24(e) shows that the power turbine speed 
drops below the engine speed setpoint. This is a result of the increased demand for thrust for forward 
acceleration. As the rotor speed drops, the forward speed control demands more blade pitch in order to 
satisfy the forward speed setpoint. Lower rotor speed means higher blade pitch angles and larger torque. 
The previous case was less severs since the forward speed was constant. In this case, extra thrust is 
demanded to produce the required forward speed acceleration. This could be mitigated by placing more 
stringent constraints on the forward speed acceleration setpoint. 
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Figure 7.24. — Variable forward cruise PSC downshift/upshift; (a) forward speed, (b) rotor 
speed, (c) left DCT shift commands, (d) right DCT shift commands, (e) power turbine 
speeds, (f) fuel rates. 

Furthermore, the flight dynamics are summarized in Figure 7.25. Figure 7.25(a) shows that the 
vertical vehicle accelerations are minimal indicating essentially no change in altitude during both the PSC 
upshift or downshift phases. In the horizontal direction, positive and negative forward accelerations are 
seen due to the forward speed setpoint ramp-up and ramp-down command. The forward acceleration is 
not constant however the deceleration is essentially constant. Figure 7.25(b) shows the closed-loop blade 
pitch angle response. During forward speed acceleration the rotor collective pitch is increased in order to 
provide the acceleration thrust. Once the increased forward speed setpoint is reached, the blade collective 
pitch remains at a higher level in order to provide the required thrust at the lowered rotor speed. Figure 
7.25(d) shows closed-loop wing flap angle response. Here the flight control system automatically reduces 
blade flap angle to maintain level flight in response to the increasing forward speed and increased lift. 

The drive system torques are summarized in Figure 7.26. As observed in the engine fuel rate response 
in Figure 7.24(f) as well as in the engine torque response in Figure 7.26(a) the engine torque increases 
during the forward speed acceleration/rotor downshift phase however, during the deceleration/rotor 
upshift phase, the torque fluctuates between negative and positive values. This is due to two competing 
factors; on the one -hand, in the deceleration phase, the flight control system automatically reduces rotor 
collective pitch (Figure 7.25(b)) to reduce forward thrust. This action results in a negative rotor inflow 
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which attempts to drive the rotor (i.e., the so-called “windmilling” effect). On the other hand, at some 
point during the upshift, the re-circulating power due to the PSC shifting strategy momentarily outweighs 
the windmilling effect briefly resulting in positive torque before again becoming negative due to the 
windmilling. This can be observed in both the engine torque plot (Figure 7.26(a)) as well as in the fuel 
rate plot (Figure 7.24(f)) over the interval [100 < t< 130] sec. 
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Figure 7.25. — Variable forward cruise PSC downshift/upshift — flight dynamics; 
(a) vehicle acceleration, (b) collective pitch, (c) induced velocity, (d) wing flap. 
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Figure 7.26. — Variable forward cruise PSC downshift/upshift — torque response; 
(a) engine torques, (b) cross-shaft torque, (c) rotor torque. 
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Figure 7.27. — Variable forward cruise PSC downshift/upshift — clutch response; 
(a) clutch torques, (b) power dissipation (clutch heating). 


Figure 7.26(b) shows the tiltrotor cross-shaft torque. In this analysis a 1 percent clutch friction 
dissimilarly between the right and left side DCTs is assumed. Prior to the start of the downshift phase, the 
cross-shaft carries no torque. However, once the shift begins, the cross-shaft begins to experience torque. 
This is exclusively due to the differential of clutch frictional characteristics between the left and right 
sides of the tiltrotor driveline. 

Figure 7.27 shows the DCT clutch torques transmitted by the high and low speed clutches and the 
corresponding dissipated power. Compared to the constant forward speed cruise/downshift case (Figure 
7.18(b) [15<f<40] sec), the power dissipated during the forward speed acceleration/PSC downshift 
(Figure 7.27(b) [15</ L <40] sec) is lower. This due to the fact that, in the latter case, rotor deceleration is 
assisted by the increased rotor drag resulting from the blade collective pitch is increase required for the 
forward speed acceleration. This has the effect of reducing overall downshift torque carried by the DCT 
clutches. One way of viewing this is that some the energy that would normally have to dissipate as heat is 
actually sent into increasing the vehicle forward kinetic energy. Likewise During the forward speed 
deceleration/rotor upshift phase (Figure 7.27(b) [100<Kl30] sec) some of the vehicle forward speed 
kinetic energy is transferred into the drive system which is then dissipated as heat during the PSC upshift. 
This additional power dissipated can be seen by comparing Figure 7.27(b) from [100<£<130] sec with 
Figure 7.18(b) from [75<7<90] sec and noting that more heat is dissipated in the forward speed 
deceleration case. 

The two-spool gas turbine engine response for this scenario is shown in Figure 7.28 to Figure 7.30. 

Figure 7.28(a) to (c) shows compressor rotational speed, air mass flow rate and turbine inlet 
temperature respectively. Figure 7.29 shows the compress running line in terms pressure ratio and mass 
flow and Figure 7.30 shows the corresponding surge margin, SM. Here, the compressor operating line 
momentarily crosses compressor surge line at the initial instant of the forward speed ramp-up and at the 
final instant of the forward speed ramp-down. This can be avoided by smoothly increasing the forward 
speed acceleration rather than the linear ramp which has “step like” acceleration. 
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Figure 7.28. — Variable forward cruise PSC downshift/upshift — engine response; 
(a) compressor speed, (b) mass flow, (c) power turbine inlet temperature. 


pressure 

ratio 



Figure 7.29. — Variable forward cruise PSC downshift/upshift — two- 
spool gas turbine engine transient running lines. 
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Figure 7.30. — Variable forward cruise PSC downshift/upshift — engine compressor 
surge margin. 
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7.3.2.2 Sequential-Shift Control (SSC) 

To avoid shifting under full power, the sequential shift control strategy, first presented in (Ref 132), is 
now explored to execute the two mission scenarios in given in Figure 7.14. Fundamentally, the sequential 
method utilizes the freewheel clutches at the power turbine output to disengage each engine sequentially 
wherein shifting along the particular driveline path can be performed in an unloaded condition. Since the 
FWU prevents the power turbine from being driven by the load, it will disengage when the load speed 
becomes larger that of the power turbine. In a multiengine driveline equipped with a variable ratio gearbox 
there are essentially two situations which can cause a specific FWU to disengage; 1) by lowering the power 
turbine setpoint speed of that engine or, 2) by downshifting the DCT gearbox in that engine’s driveline path. 
In the first situation, the remaining engines maintain the system speeds as the particular engine speed is 
lowered. In second the situation, the commanded downshift acts as an upshift on the DCT gearbox input 
side thus increasing the load speed at the FWU. Since the engine setpoint speed remains fixed, the FWU in 
this driveline path will disengage. Both of these situations arise in the SSC shifting method. 

7.3.2.2.1 Steady Forward Cruise PSC Downshift/Upshift 

The timeline of the SSC downshift procedure along with a brief description of each event is given in 
Table 7.7. 

Figure 7.31 to Figure 7.36 give the forward cruise SSC downshift response of the two-speed DCT 
tiltrotor system. By examining Figure 7.31(a) to (b), it is apparent that the rotor speed is successfully 
lowered with little deviation in forward velocity. 

Figure 7.31(c) to (d) show the open-loop DCT clutch pressure command inputs used to sequentially 
downshift the left and right side DCTs (steps 2 and 4 in Table 7.7). Furthermore, the engine power turbine 
speeds are shown in Figure 7.31(e) to (f). 


TABLE 7.7.— SSC DOWNSHIFT TIMELINE* 


Time 

Commands 

Comments 

£=12 sec 

1) Disable engine torque sharing control loop 

2) Start left side DCT downshift (disengage clutch 1 and 
engage clutch 2). 

Pd\ - -30psi/sec 
p d 2 = +4psi/sec 

Start of SSC downshift procedure 
Increases left DCT input shaft speed which 
disengages left engine FWU 

£=32 sec 

3) Start right engine power turbine setpoint ramp-down 
(-600 RPM/sec) 

Main rotor speed reduces from 190 to 102.5 RPM 
Left engine re-engages FWU at end of ramp-down 

£=52 sec 

4) Start right side DCT downshift (disengage clutch 1 
and engage clutch 2). 
p cl i = -30psi/sec 
Pd 2 = +4psi/sec 

Increases right DCT input shaft speed which 
disengages right engine FWU 

£=72 sec 

5) Start right engine power turbine setpoint ramp-up 
(+600 RPM/sec) 

Right engine re-engages FWU at end of ramp-up 

£=82 sec 

6) Re-enable engine torque sharing control loop 

End of SSC downshift procedure 


*Under constant forward speed setpoint v c>set = 228 kn 


In Figure 7.3 1(e) it is observed that, even though the left side engine setpoint speed remains constant, 
the power turbine transiently increases due to the left side DCT downshift. This is a result of the rapid 
unloading due to FWU disengagement (see step 2 in Table 7.7). According to steps 3 and 5 in Table 7.7, 
Figure 7.3 1(f) shows the ramp-down and ramp-up of the right engine power turbine. Similarly, when the 
right side DCT is downshifted, the right engine power turbine experiences at transient speed increase above 
its setpoint as its FWU disengages (see step 4 in Table 7.7). Furthermore, Figure 7.31(g) shows the 
corresponding fuel rates commanded by the closed-loop engine fuel control for both the left and right side 
engines. Note, the left and right side engine fuel rates are equivalent prior to disabling (£=12 sec) and after 
re-enabling (£=82 sec) the multiengine torque sharing control loop (see steps 1 and 6 in Table 7.7). 


NAS A/CR— 20 13-21 6502 


143 





t, sec 

Figure 7.31 . — Steady cruise SSC downshift; (a) forward speed, (b) rotor speed, (c) left 
DCT shift commands, (d) right DCT shift commands, (e) left power turbine speed, 

(f) right power turbine speed (g) fuel rates. 


The flight dynamics are summarized in Figure 7.32. Figure 7.32(a) shows that both the horizontal and 
vertical vehicle accelerations are minimal. Thus, no forward speed or altitude is lost during the downshift. In 
order to maintain the constant forward speed setpoint under the decreasing rotor speed, the closed-loop 
blade pitch controller increases the blade collective pitch to keep constant thrust as shown Figure 7.32(b). 
This results in a corresponding increase in rotor induced velocity is seen in Figure 7.32(c). Finally, as seen 
in Figure 7.32(d), under this constant thrust condition, there is no appreciable change in wing flap angle. 

The system torque response is summarized in Figure 7.33. As observed in the engine fuel rate 
response in Figure 7.31(g), the left and right side power turbine torques shown in Figure 7.33(a) are 
identical prior to *=12 sec and after t= 82 sec. During the downshift, the engine torques become 
dissimilar and are due to both rotor aerodynamic torque as well as acceleration torque demands. As a 
result of the dissimilar left and right engine torques, the cross-shaft also becomes loaded as shown in 
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Figure 7.33(b). Also, Figure 7.33(c) shows the increasing rotor torque which is a result from the 
increasing collective pitch and constant thrust requirement (see Figure 7.32(b)). 
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Figure 7.32. — Steady cruise SSC downshift — flight dynamics; (a) vehicle acceleration, 
(b) collective pitch, (c) induced velocity, (d) wing flap angle. 
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Figure 7.33. — Steady cruise SSC downshift — torque response; (a) engine torques, 
(b) cross-shaft torque, (c) rotor torque. 
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Finally, more details of the two-spool gas turbine engine response are given in Figure 7.34 and Figure 

7.35. In particular, by plotting compressor pressure ratio versus air mass flow, the left and right engine 
transient operating lines during the SSC downshift are shown in Figure 7.35. Figure 7.35 shows that both 
engines experience a wide range of operation during the SSC process and that both engines transiently 
pass over the surge line. In this analysis, three different DCT shift-rate cases (slow, medium and fast) as 
given in Table 7.8 are considered. Note the medium case corresponds to the values used thus far in Figure 
7.31 to Figure 7.35. 

The left and right engine SM values for the slow, medium and fast shift-rates are shown in Figure 

7.36. In each case, the surge margin is crossed briefly. Here, Figure 7.36(a) shows that the left engine 
crosses the surge line immediately after t= 52 sec as the right DCT begins to downshift (step 4, Table 7.7). 
Thus, left engine surge is a result of the transient shock due to unloading the right engine FWU. 

Furthermore, Figure 7.36(b) shows that the right engine experiences a brief surge near the start and 
end of the downshift process for similar reasons as the left engine. Specifically, the first surge of the right 
engine is due to the initial disengagement of the left engine FWU (step 1, Table 7.7) and the second surge 
occurs when the left engine re-engages the FWU at the end of its ramp-up (step 5, Table 7.7). Figure 7.36 
shows that slower DCT shifting rates are more beneficial for the engine surge margins in the SSC 
downshift procedure. 
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Figure 7.34. — Steady cruise SSC downshift — engine response; (a) compressor 
speed, (b) mass flow, (c) power turbine inlet temperature. 
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Figure 7.35. — Steady cruise SSC downshift — gas turbine engine 
transient running lines. 


TABLE 7.8.— DUAL CLUTCH TRANSMISSION DOWNSHIFT RATE CASES FOR SSC 


Case 

Clutch pressure engagement and disengagement rates 

Slow 

Clutch 1 disengagement rate, p dx = -l5 psi/sec 

Clutch 2 engagement rate, p d2 = 4 psi/sec 

Medium 

Clutch 1 disengagement rate, p dx = -30 psi/sec 

Clutch 2 engagement rate, p d2 = 4 psi/sec 

Fast 

Clutch 1 disengagement rate, p dx = -300 psi/sec 

Clutch 2 engagement rate, p dl = 4 psi/sec 



Figure 7.36. — Steady cruise SSC downshift — engine compressor surge margins for 
three DCT shift-rates. 
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Next, the simulation results of the SSC upshift procedure, which is essentially the reverse process of 
the SSC downshift, are described. The SSC upshift timeline is given in Table 7.9. Figure 7.37 to Figure 
7.42 show the constant forward cruise SSC upshift response of the two-speed DCT tiltrotor. By 
examining the forward speed and the rotor speed in Figure 7.37(a) and (b), it is apparent that rotor speed 
is successfully increased with little deviation in forward velocity. Figure 7.37(c) and (d) show the open- 
loop DCT clutch pressure command inputs used to sequentially upshift the right and left side DCTs (steps 
3 and 5 in Table 7.9). Furthermore, the engine power turbine speeds are shown in Figure 7.37(e) to (f). In 
Figure 7.37(e) it is observed that the left power turbine speed transiently increases due to the right engine 
ramp-up. This is due to the unloading of the left engine FWU. (See step 4 in Table 7.9). According to 
steps 2 and 4 in Table 7.9, Figure 7.37(f) shows the ramp-down and ramp-up of the right engine power 
turbine. Figure 7.37(g) shows the corresponding fuel rates commanded by the closed-loop engine fuel 
control for both the left and right engines. Note, the left and right side engine fuel rates are equivalent 
prior to disabling (t=l sec) and after re-enabling (t = 92 sec) the multiengine torque sharing control loop 
(see steps 1 and 6 in Table 7.9). Also note that the fuel rate is essentially flat during the rotor speed 
acceleration phase (step 4, Table 7.9). This is due to the fact that as rotor speed increases, the collective 
pitch is automatically decreased by the forward speed pitch controller to maintain constant thrust as 
further shown in Figure 7.38. 

Similar to the downshift case, Figure 7.38(a) shows that both the horizontal and vertical vehicle 
accelerations are minimal thus with no forward speed or altitude lost during the upshift. To maintain 
constant forward speed under the increasing rotor speed, the closed-loop blade pitch controller decreases 
the blade collective pitch to keep constant thrust as shown Figure 7.38(b). This causes the drop in induced 
velocity as seen in Figure 7.38(c). Finally, as seen in Figure 7.38(d), under the constant thrust condition, 
there is no appreciable change in wing flap angle. 

The driveline torques are summarized in Figure 7.39. As observed in the engine fuel rate response in 
Figure 7.37(g), the left and right side power turbine torques are identical prior to t=l sec and after t = 92 
sec. During the downshift, the engine torques become dissimilar and are due to both rotor aerodynamic 
torque as well as acceleration torque demand 

As a result of the dissimilar left and right engine torques, the cross-shaft also becomes loaded as 
shown in Figure 7.39(b). Here, the cross-shaft is subjected to a fully reversing torque load. This is also 
true in the downshift case. Finally, Figure 7.39(c) shows the deceasing rotor torque which is a result from 
the decreasing collective pitch. 


TABLE 7.9.— SSC UPSHIFT TIMELINE* 


Time 

Commands 

Comments 

t=l sec 

1) Disable engine torque sharing control loop 

2) Start right engine power turbine setpoint ramp-down 
(-600 RPM/sec) 

Start of SSC upshift procedure 
Disengages right engine FWU 

t=21 sec 

3) Start right side DCT upshift (disengage clutch 2 and 
p r n = +4 psi/sec 

engage clutch 1). 

P c i 2 = 30 psi/sec 

Decreases right DCT input shaft speed. 
Right engine FWU becomes re-engaged 

t=Al sec 

4) Start right engine power turbine setpoint ramp-up 
(+230 RPM/sec) 

Main rotor speed increases from 102.5 to 190 RPM 
Speed increase disengages left engine FWU 

t= 82 sec 

5) Start left side DCT upshift (disengage clutch 2 and 
b rU = +4 psi/sec 

engage clutch 1). 

p d2 = -30 psi/sec 

Decreases left DCT input shaft speed which re- 
engages left engine FWU 

t = 92 sec 

6) Re-enable engine torque sharing control loop 

End of SSC upshift procedure 


*Under constant forward speed setpoint v c>set = 228 kn 
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t, sec 

Figure 7.37. — Steady cruise SSC upshift; (a) forward speed, (b) rotor speed, (c) 
left DCT shift commands, (d) right DCT shift commands, (e) left power turbine 
speed, (f) right power turbine speed (g) fuel rates. 


NAS A/CR— 20 13-21 6502 


149 


0.05 

vehicle 

accelerations 0 
(g) 


- 0.05 



Horizontal - 
— Vertical 

(a) 

1 1 1 1 1 1 1 1 1 


rotor 

collective 

pitch 


rotor 

induced 

velocity 

(m/s) 

wing 

flap 

pitch 


0 


10 20 30 40 


50 


60 


70 


80 90 


100 



t , sec 

Figure 7.38. — Steady cruise SSC upshift — flight dynamics; (a) vehicle 
acceleration, (b) collective pitch, (c) induced velocity, (d) wing flap angle. 
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Figure 7.39. — Steady cruise SSC upshift — torque response; (a) engine torques, 
(b) cross-shaft torque, (c) rotor torque. 
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Other important aspects of the two-spool gas turbine engine response are shown in Figure 7.40 and 
Figure 7.41. As seen in Figure 7.41, similar to the downshift case, the transient running lines of each 
engine undergo small excursions past the compressor surge line during the SSC upshift process. 
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Figure 7.40. — Steady cruise SSC upshift — engine response; (a) compressor speed, 
(b) mass flow, (c) power turbine inlet temperature. 
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Figure 7.41. — Steady cruise SSC upshift — gas turbine engine transient running lines. 
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Next, as in the downshift case, three different DCT shift-rates (slow, medium and fast), given by 
Table 7.10, are examined. 

The resulting SM values for each shift-rate are shown in Figure 7.42. In each case, the SM briefly 
crosses the surge line. Figure 7.42(a) shows that the left engine crosses the surge line when the FWU re- 
engages due to the left side DCT upshift (step 5, Table 7.9). Furthermore the right engine experiences a 
surge near t = 30 sec and t = 48 sec. The first surge of the right engine is due to re-engagement of the right 
engine FWU as a result of the right side DCT upshift (step 3, Table 7.9). The second surge, near t = 48 
sec, is due to the right engine power turbine setpoint ramp-up (step 4, Table 7.9) which causes the left 
engine to disengage from its FWU. Finally, Figure 7.42 shows that faster DCT shifting rates are more 
beneficial for the engine surge margins in the SSC upshift procedure. 



t , sec 

Figure 7.42. — Steady cruise SSC upshift — engine compressor surge margins for three 
DCT shift-rates; (a) left side, (b) right side. 


TABLE 7.10.— DUAL CLUTCH TRANSMISSION UPSHIFT RATE CASES FOR SSC 


Case 

Clutch pressure engagement and disengagement rates 

Slow 

Clutch 1 engagement rate, p dl = 4 psi/sec 

Clutch 2 disengagement rate, p d2 = -15 psi/sec 

Medium 

Clutch 1 engagement rate, p dl = 4 psi/sec 

Clutch 2 disengagement rate, p d2 = -30 psi/sec 

Fast 

Clutch 1 engagement rate, p dl = 4 psi/sec 

Clutch 2 disengagement rate, p d2 = -300 psi/sec 
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73.2.2.2 


Variable Forward Speed Cruise SSC Downshift/Upshift 

Figure 7.43 to Figure 7.49 gives the simulation results for the two-speed DCT tiltrotor executing an 
SSC downshift followed by an SSC upshift during the variable forward speed cruise condition depicted in 
Figure 7.14(b). 

By examining Figure 7.43(a) to (b), starting at t = 30 sec it seen that the forward cruising speed is 
increased from 150 to 300 kn while the rotor speed is simultaneously downshifted from about 190 to 102 
RPM. Then, starting at t = 195 sec it seen that the forward speed is increased back to 300 kn while the 
rotor is simultaneously upshifted back to 190 RPM. Figure 7.43(c) to (d) show the open-loop left and 
right side DCT clutch pressure command inputs. Here note that the left and right sides are shifted 
sequentially rather than simultaneously. 
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Figure 7.43. — Variable forward cruise SSC downshift/upshift; (a) forward speed, (b) rotor 
speed, (c) left DCT shift commands, (d) right DCT shift commands, (e) left power turbine 
speed, (f) right power turbine speed (g) fuel rates. 
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In addition, the left and right side engine power turbine speeds are shown in Figure 7.43(e) and (f) 
and the corresponding combustion chamber fuel rates commanded by the closed-loop engine fuel control 
are shown in Figure 7.43(g). During the SSC downshift/forward speed acceleration phase, Figure 7.43(e) 
shows that the left side power turbine speed transiently increases above the setpoint due to the left side 
DCT downshift and resulting FWU disengagement. Furthermore, Figure 7.43(f) shows the ramp-down 
and ramp-up of the right engine power turbine as required by the SSC approach. Similarly, when the right 
side DCT is downshifted, the right engine power turbine experiences at transient speed increase above its 
setpoint as its FWU disengages. Also, Figure 7.43(g) shows the corresponding fuel rates commanded by 
the closed-loop engine fuel control for both the left and right side engines. Note, the left and right side 
engine fuel rates are equivalent prior to disabling and after re-enabling the multiengine torque sharing 
control loop. Figure 7.43(g) also shows that the engine fuel rate become saturated at their upper-limit 
during the forward speed acceleration phase. 

Furthermore, the flight dynamics are summarized in Figure 7.44. Figure 7.44(a) shows that the 
vertical vehicle accelerations are minimal indicating essentially no change in altitude during both the SSC 
upshift or downshift phases. In the horizontal direction, positive and negative forward accelerations are 
seen due to the forward speed setpoint ramp-up and ramp-down command. Unlike the PSC case, SSC 
shifting achieves nearly constant forward vehicle accelerations and decelerations. 

The drive system torques are summarized in Figure 7.45 and the DCT clutch torques and power 
dissipation levels are shown in Figure 7.46. By comparing clutch power dissipation levels under the PSC 
and SCC strategies, (see Figure 7.27(b) and Figure 7.46(c) to (d)), it is apparent that the SSC approach 
generates significantly less frictional heating. Finally, Figure 7.47 to Figure 7.49 give more details of the 
two-spool gas turbine engine response. Once again, the engine experience brief transient excursions 
across the surge line during transitions. 
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Figure 7.44. — Variable forward cruise SSC downshift/upshift — flight dynamics; 
(a) vehicle acceleration, (b) collective pitch, (c) induced velocity, (d) wing flap. 
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Figure 7.45. — Variable forward cruise PSC downshift/upshift — torque response; (a) 
engine torques, (b) cross-shaft torque, (c) rotor torque. 
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Figure 7.46. — Variable forward cruise SSC downshift/upshift — clutch response; 

(a) left side clutch torques, (b) right side clutch torques, (c) left side clutch power 
dissipation, (d) right side clutch power dissipation. 
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Figure 7.47. — Variable forward cruise SSC downshift/upshift — engine response; 
(a) compressor speed, (b) mass flow, (c) power turbine inlet temperature. 



Figure 7.48. — Variable forward cruise SSC downshift/upshift — 
two-spool gas turbine engine transient running lines. 
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Figure 7.49. — Variable forward cruise SSC downshift/upshift — engine compressor 
surge margin; (a) left side, (b) right side. 


7.3.2.3 Effect of Tiltrotor Driveline Topology 

In this section the effect of cross-shaft location on the two speed shift response is briefly explored. In 
general it is possible to have the cross-shaft upstream or downstream from the speed changing gearbox or 
DCTs as shown in Figure 7.50. In all the tiltrotor results of the previous subsections in this section, the 
cross-shaft is located downstream from the dual-clutch gearboxes, (see Figure 7.50(a)). In the upstream 
configuration the cross-shaft will experience torques produced by differentials in left and right side DCT 
clutch friction properties See Figure 7.17(b) and Figure 7.26(b). Note, it is only possible to implement the 
SSC strategy on the upstream configuration since the DCTs must be able to be decoupled via the engine 
FWUs. However, it is possible to implement the PSC strategy with either an upstream or downstream 
cross-shaft configuration. 

Figure 7.51 compares the cross-shaft torques under the steady cruise PSC upshift/downshift scenario 
for an upstream and downstream cross-shaft configuration for two values of clutch friction dissimilarly. 
Prior to the start of the downshift phase, the cross-shaft carries no torque. However, once the downshift 
begins the cross-shaft begins to experience torque. This is exclusively due to the differential clutch 
frictional torques between the left and right sides of the tiltrotor driveline. Most notably, in the 
downstream configuration, after the shift is completed, a static torque remains in the cross shaft. This is a 
result of the differential windup between left and right side clutch systems prior to clutch lockup. Once 
the clutches on both sides become locked, the static torque is now built into the system until the next 
shifting event or until the rotors becomes unloaded. 
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Figure 7.50. — Two-speed tiltrotor design cases; (a) cross-shaft downstream from DCT gearboxes; (b) cross-shaft 
upstream from DCT gearboxes. 
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Figure 7.51 . — Steady cruise PSC downshift/upshift cross-shaft torques for two 
levels of clutch friction dissimilarity; (a) Downstream configuration, 

(b) Upstream configuration. 


7.4 Summary and Conclusions 

The Comprehensive Variable-Speed Rotorcraft Propulsion Modeling (CVSRPM) tool developed in 
this research accounts for coupled rotor/engine/fuel control/gearbox/shaft/clutch/flight control subsystem 
dynamic interactions under hover and forward speed conditions. In this investigation, a prototypical two- 
speed Dual-Clutch Transmission (DCT) is proposed and designed to achieve 50 percent rotor speed 
variation. Utilizing the modeling and simulation tools developed in this investigation, variable-speed 
operation and transient shift response of both single path (single rotor helicopter) and dual path (twin 
rotor tiltrotor) configurations are analyzed. Here, the two-speed DCT is utilized as the variable speed 
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elements in each drive-system. In the case of the tiltrotor configuration, the simulation parameters are 
based on the NASA LCTR-2 prototype operating in forward-flight airplane mode. 

Here both Parallel Shifting Control (PSC) and Sequential Shifting Control (SSC) strategies are 
explored. In the PSC strategy, all branches of the multipath drive-system are shifted simultaneously, 
while under SSC, each branch is shifted in sequence to avoided shifting under full power conditions. The 
chief advantage of PSC is that gear ratio changes can made relatively rapidly without need to involve 
engine control during the shift. However, for large rotorcraft, the clutch power dissipation and thermal 
energy dissipation requirements may be prohibitive. In SSC, power to a specific branch is removed by 
reducing the setpoint speed of the associated engine thereby disengaging the overrunning clutch and 
enabling the transmission to shifted under minimal load. Thus, SSC significantly reduces clutch power 
dissipation and hence clutch heating. However the SSC strategy requires a multiengine/multi DCT drive 
system topology whereas the PSC strategy can be implemented on singe DCT and multi DCT 
architectures. 

One main contribution of this investigation to the existing body of research is the evaluation of both 
the PSC and SSC shift strategies in a comprehensive rotorcraft drive system simulation model. In 
particular, this investigation builds upon previous variable-speed rotorcraft propulsion studies by 1) 
including a fully transient gas-turbine engine model, 2) including clutch stick-slip friction effects, 3) 
including driveshaft flexibility, and 4), incorporating a basic, fully trimmed, flight dynamics model to 
account for interactions with the flight and engine control systems. 

In PSC, the engaging clutch pressure is linearly ramped up while the disengaging clutch pressure is 
simultaneously commanded to ramp down all while maintaining the engine setpoint speed of all engines. 
One of the key parameters in PSC is the clutch engagement/disengagement ramp rates. Under the PSC 
strategy, selecting clutch shift-rate results in a design tradeoff between transient engine surge margins and 
clutch frictional power dissipation. 

In SSC, the sequential shifting provided smooth rotor speed changes while maintaining a constant 
forward speed setpoint and constant altitude. However, it was found that transient shocks produced by 
engine engagements and disengagements from their freewheeling clutches tended to produce brief 
excursions of the engine compressor surge limit. It was found that proper choice of dual-clutch 
transmission shift rate could minimize these excursions. In the sequential downshift case, it was found 
that slower shift rates tended to reduce the engine surge effects while the opposite trend was found for the 
pshift case. 

Given the current trend toward tandem, tilt-rotor, multirotor, and co-axial rotor/pusher-prop 
configurations with numerous cross-shafting and multipath power flow arrangements (e.g., Boeing CH-47 
and V-22, Sikorsky X-2 High Speed Lifter and Heavy Lift quad tiltrotor concepts), It is hoped that the 
analysis tools and results presented from the CVSRPM code will be utilized as a guide for investigating 
and evaluating new configurations and designs for fixed and variable speed rotorcraft under a variety of 
operating conditions. 
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Appendix. — Comprehensive Variable Speed 
Rotorcraft Propulsion System Model 



Figure A.1. — Two-speed rotorcraft driveline comprehensive simulation — main block. 
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Figure A.2. — Two-speed tiltrotor driveline/engine/clutch/rotor/fuel control/flight dynamics block. 
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Figure A.3. — 2-DOF (altitude and forward speed) vehicle flight dynamics block. 



Rotor blade pitch 


Figure A.4. — Rigid rotor and blade pitch forward speed controller block. 


Here note the gas turbine engine model, engine fuel control and dual clutch transmission model 
simulation blocks are given in Sections 4.0 and 6.0, respectively. 
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